NASA CONTRACTOR 
REPORT 


CO 

CO 

04 



0 £ 


<X 


ANALYSIS OF THE FEASIBILITY 
OF AN EXPERIMENT TO MEASURE 
CARBON MONOXIDE IN THE ATMOSPHERE 


by M. H. Bortner, F. N. Alyea, R. N. Grenda, 
G. R. Liebling, and G. M. Levy 


Prepared by 

GENERAL ELECTRIC COMPANY 
Philadelphia, Pa. 19101 
for Langley Research Center 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • OCTOBER 1973 


1. Report No. 2. Government Accession No. 

NASA CR-2303 

3. Recipient's Catalog No. 

4. Title and Subtitle 

ANALYSIS OP THE FEASIBILITY OF AN EXPERIMENT TO MEASURE CARBON 
MONOXIDE IN THE ATMOSPHERE 

5. Report Date 

October 1973 

6. Performing Organization Code 

7. Author(s) 

M. H. Bortner, F. H. Alyea, R. N. Grenda, G. R. Llebllng and 
G. M. Levy (Barringer Research Limited) 

8. Performing Organization Report No. 

10. Work Unit No. 

630-52-00-01 

9. Performing Organization Name and Address 

General Electric Company 

Space Division - Space Sciences Laboratory 

P. 0. Box 8555 

Philadelphia, PA 19101 

11. Contract or Grant No. 

NAS 1-10139 

13. Type of Report and Period Covered 

Contractor Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


This is a topical report. 

16. Abstract 

An analysis has been carried out to determine the feasibility of measuring atmospheric carbon 
monoxide from a remote platform using the correlation Interferometry technique. The Carbon 
Monoxide Pollution Experiment (COPE) has been found to be a feasible method for obtaining a global 
carbon monoxide map and a vertical carbon monoxide profile. It has been determined that CO 
data can be obtained with an accuracy of 10 percent using this technique on the first overtone 
band of CO at 2.3^i . That band has been found to be much store suitable than the stronger 
fundamental band at 4.6^>-. Calculations for both wavelengths are presented which illustrate the 
effects of atsiospheric temperature profiles, inversion layers, ground temperature and emlsaivity, 
CO profile, reflectivity, and atmospheric pressure. The applicable radiative transfer theory on 


which these calculations are based is described together with the principles of the technique. 


17. Key Words (Suggested by Author(s)| 

Carbon monoxide 
Atmospheric measurement 
Air pollution 
Atmospheric chemistry 
Interferometer 
Correlation 


18. Distribution Statement 

Unclassified - Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 




Domestic, $6.25 

Unclassified 

Unclassified 

261 

Foreign. $8.75 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 


i 





















FOREWORD 

This report was prepared for NASA as part of contract NASI -10 139 
with Langley Research Center under the Advanced Applications Flight Ex- 
periments (AAFE) Program. The objective of this contract is the develop- 
ment of the Carbon Monoxide Pollution Experiment ("COPE"). This experi- 
ment is designed to obtain data for the investigation of mechanisms by which 
CO is removed from the earth's atmosphere. The approach uses an orbiting 
platform to remotely map global CO concentrations and determine vertical 
CO profiles using a correlation interferometer measurement technique being 
developed by Barringer Research Ltd. The instrument is to be capable of 
measuring CO over the range .02 to 20 atm. -cm. and of measuring other 
trace atmospheric constituents. 

The report covers one aspect of the study - the analysis of the feasi- 
bility of the experiment. A previous report was concerned with the CO prob- 
lem and a forthcoming report will be concerned with the breadboard instru- 
ment - its design, fabrication and testing. The technique is, of course, 
suitable for the measurement of other atmospheric trace species, and such 
applications are currently under study. 

The authors would like to express their appreciation to several co- 
workers who have contributed significantly to the' work described herein and 
to related work - specifically, to S. H. Chandra who did most of the work on 
scintillation effects described herein; to D, N. Vachon who provided meteo- 
rological information such as that on precipitable water and cloud cover; to 
J. C. Burns who made numerous helpful suggestions and reviewed much of 
the work; to S. L. Neste who assisted with the calculations; to H. W. Gold- 
stein who reviewed and evaluated much of the work; and R. H. Kummler who 
contributed many valuable suggestions in his consultation throughout the work. 
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ANALYSIS OF THE FEASIBILITY Qf Ai^ 

EXPERIMENT TO MEASURE 
CARBON MONOXIDE IN THE ATMOSPHERE 

By M, H. Bortner, F.N. Alyea, R.N. Grenda. 
and G. R. Eiebling, General Electric Cp. , 

Space Sciences Laboratory, and G. M. Levy, 

Barringer Research Ltd. 

1. INTRODUCTION 

One of the majqr problems in the field of remote sensing of atmospher- 
ic pollutant concentrations consists of determining what measurement? should 
be made and what information can be obtained from a given meas virement. It 
is, therefore, economically advisable to simulate (to the extent possible) mea- 
surements with the computer and utilize the results to determine the most 
useful instrument measurements. This report describes passive measure- 
ment of pollutants via their absorption of sunlight. The simulation consists 
of defining a representative model atmosphere and subsequently solving the 
radiative transfer equation to obtain the spectrum falling on the instrument 
and then calculating the action of the instrument to determine ^he signal pro- 
duced. Effects suc;h as reflected sunlight, earthshine, atmospheric absorp- 
tion and atmospheric emission must be considered in order to yield physically 
meaningful solutions tp the equation. The specific ppllutpr* py chemical spe- 
cies of interest v/ill de^ermipe the spectral region for which the atmospheric 
transmission will be computed. For the case in \yhicl> php InstrveTipnt is an 
interferometer, the Fourier transform of the spectrum is then calculated to 
produce an interferogram. Treatment of the interferograrn produces the 
measurement of the pollutant. 

The various portions of this type of calculation are described herein as 
applied to the calculation of the measurement of atmospheric carbon monox- 



ide by the correlation interferometer. The results of calculations are dis- 
cussed. 

The model has been used to investigate the feasibility of the use of the 
correlation interferometer to measure GO, to consider the usefulness of the 
measurement in looking for a CO sink* and to define the exact spectral band 
which should be used for most reliable and useful results. The measurements 
to be made include a' mapping of CO concentrations by looking downward to- 
ward the earth and a determination of the vertical CO profile by looking at the 
sun through the earth's limb, as described in Section 4. 

The model is the basis for a program which computes the upwelling 
radiation for both the mapping and the limb modes. It calculates the spec- 
trum incident on the instrument and the resulting interferogram. Another 
program then is used to invert the results to obtain CO densities. The model 
has been used to calculate transmission to determine the range of sensitivity, 
the effect of various atmospheric and source parameters, and the effect of 
interfering species. 


Bortner, Kummler, and Jaffe (1972) have reviewed the sources, sinks, and 
concentrations of carbon monoxide in the earth's atmosphere. 


2 



2 . 


ATMOSPHERIC MODELS 


2. 1 Composition Models 

The analytical work performed employed several atmospheric 
models involving several variations of the carbon monoxide, water, and tem- 
perature profiles and one carbon dioxide and one methane profile. 

Specifically, ten CO profiles were used. These are given in 
Table 2. 1. 1 and shown graphically in Figure 2. 1.1. Profile 1 is that for a 
constant mixing ratio of 0. 1 ppm; profile 2 is that for 10 ppm up to' 2 km and 
0. 1 ppm above that; profile 3 is that for a constant mixing ratio of 0. 01 ppm; 
profile 4 is that representing a sink in the 20-45 km range with the mixing 
ratio dropping from 0. 1 to 0.01 ppm over this altitude range; profile 5 is that 
representing a low altitude sink with an effect up to 9 km, having a mixing 
ratio of 0. 01 ppm at 0 km and 0. 1 ppm at 9 km; profile 6 represents a low 
altitude sink with an effect up to 1 km, having a mixing ratio of 0.01 ppm at 
0 km increasing to 0. 1 ppm at 1 km; profile 7 also represents a low altitude 
sink with an effect up to 3 km, having a mixing ratio of 0. 005 ppm at 0 km 
and 0. 1 ppm at 3 km; profile 8 represents a low altitude sink with an effect 
up to 1 km as in profile 6, but having a mixing ratio of 0. 01 ppm up to 1 km 
and 0. 1 ppm above. It is important for the mapping measurements to show 
the effects of these differences except for profile 4, the effect of which should 
be shown in the limb experiment. The total number densities from which the 
CO number densities were obtained with the above mixing ratios were taken 
from Bortner and Kummler (1968) in which they were derived mainly from 
CIRA (1965) representing mean conditions throughout the year for latitude 
near 30°. 

For purposes of making calculations (Section 7. 3) to compare 
with those of the Convair group (Ludwig, 1970) CO models 9 and 10, Table 2. 1. 1 


3 



- — - 

- — ■ 

- — - 



x-4 

x-4 

x«4 

x-4 

X-4 

.x—. 

x-4 

x-4 

x— 
















*“ ■ 1 






•-* 

PM 



— H 

o 

o 

x-4 


x-4 

X-4 

_ 

x-4 

x— . 

X4. 

X— * 

x-4 

x-4 x-4 

an. 

x— 

x-» 





>—4 

— t 


r— « 

-H 

— < 

M 


—4 

«— < 

O' 

O' 

O' 

O' 

00 

CO 

00 

r- 

r- 


r- 3 

M3 

in 

in 

'4-* 




4x* 

4 — • 

«w> 

4«X 

«»- 

4x" 

4x- 

4«X 

4X» 

4x- 















in 

00 

pm 

pM 

m 

Csl 

PM 

00 

PM 

PO 

m 


m 

00 

00 

o 

00 

PM 

in 

<M 

O' 

co 

m 

00 

Tf pm 

m 

o 

O 

m 

r- 

PM 

r- 

(M 

00 


o 

r- 



o 

nD 

o 

m 

*<* 

o 

O 

m 

O' 

in 

M3 

in 

PM 

o -• 

M3 

PO 

vO 

vO 

m 

m 


V 

rd 

PO 

PO 

PM 

rM 

no 


M* 

pm 

O' 


CM 

-I 

in 

PM 


oo 


tM 

•d M* 

-d 

sd 

PM 


cr| 


PM 

iM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

^4 

—i 

o 

© 

o 

O 












~ 4 

"4 


4 -* 

•— i 




— * 

<—» 





«-4 

>“H 


—* 

o 

00 

00 

00 

oo 

00 

00 

r-~ 

r- 

MD 

s£> 

4 - ' 

r" 

4 —' 

4— ' 

4«* 



4««x 

4-- 

4 — »■ 

4x*" 

4 — > 

4x- 

*»«» 

4x- 















m 

00 

PM 

PM 

m 

PM 

PM 

00 

PM 

PO 

in 

— < 

PO 

00 

00 

o 


PM 

m 

PM 

O' 

ao 

in 

00 


PM 

m 

O 

o 

PO 

r- 

PM 

h> 

PM 

00 


o 

r- 

M* 


O 

M3 

o 

in 


o 

O 

ro 

O' 

m 

vO 

m 

PM 

O 


vO 

<o 

sO 

sd 

in 

m 


xf 

PO 

fO 

PO 

PM 

PM 

PM 

-1 


PM 

O' 


PM 


m 

PM 


00 

M* 

PM 


yf 


sd 

PM 


00| 


cO 

i 


a 

o 

CO 

W 

t-H 

h 

O 

Dh 

o 

o 


(M 

w 

CQ 

< 

H 


•'M 


*1 


"M 



— 


PM 

PM 


PM 


PM 

PM 

PM 

PM 

—4 

—4 

x-4 

— i 

o 

o 

O 


X-^ 

X— 

XI. 

X4. 

X-4 

X-V 

X4 





XX4 

X4 

— ^ 

— < 


— < 

— H 


4-4 



^-4 

x-4 

— H 


—4 

—4 

—4 

—4 

X4 

H 

O' 

O' 

O' 

O' 

00 

00 

00 

r- 

r- 


sO 

M3 

M3 

— 






"«W 


•— X 

X— x 

*— x 

Xwx 

•w 

4.^* 

*w» 

■*— ' 

*— x 

4X 

4— X 














rr 



O 

O' 


o 


m 

r- 

m 

O' 

—4 

o 

m 

in 

m 

ro 

sO 

—4 

O' 


r- 

M) 

n- 

PM 

— i 

vO 

in 

O' 

PM 


m 

PO 


o 

00 


r~ 


m 

m 

PM 

o 

r- 

SO 

o 

00 

m 

CD 

r- 

ro 

o 

-* 


m 

M* 

00 

o 

-4 

M3 

in 

in 

O 

PM 

PM* 


<M 

— 




-4 


-X 

-* 

O' 

od 


- 

od 

PM 

- 

00 

M* 

(M 

— 

sd 

ro 

— 

o 


— 

M3 

PM 

-* 



fc 






























— < 

— < fM 


PM 

PM PM 


PM 


PM 

PM 

PM 

PM 


x^ 

x-4 

xH 

o 

O 

o 


X-4 

X-4 

X-4 

X-4 

X-4 

X-4 

x-4 

X-4 

X-4 

X4 

X-4 

X-4 

— H 

r—t 

■ — 1 


^ ■ — « 




x-< 

x-4 

—4 

—4 

—4 

x4 

x-4 

x-4 

x-4 

x4 

x-4 

O' 

O' 

O' 

O' 

00 

ao 

00 

r- 

r- 

r- 

M) 

M> 

NO 





"*4* 


•>^x 


■X— ' 

xw » 

«— x 

4_X 

x— x 

XX 

4— X 

x— x 

x— - 

X— —X 















r- 

PO PM 

o 

— t 

O O' 


O 


pn 

r- 

m 

O' 

—4 

o 

m 

in 

m 

ro 

vO 

X, 

O' 

M* 

r- 

M3 

n- 

PM 

—< 

M3 

in 

O' 

PM 

M* 

PM 

M 4 MD 

o 

PM 

sO 00 


n- 


in 

m 

PM 

o 

n- 

s£) 

o 

ao 

m 

OD 

n- 

ro 

© 


—• 

ro 


00 

© 


M3 

in 

n 

O 

— 

PM 

O' 


— 




— 

-d 



O' 

oo 

if 


od 

ro 

x4 

od 

M* 

PM 


m> 

ro 

— 

O' 



sd 

PM 

— 


































PM 


PM 

(M 


PM 


PM 

PM 

PM 

PM 

■— 4 

—4 

-4 

—4 

© 

O 

o 


X— 

X-4 

X-4 

X-4 

^4 

X-4 

x-4 

xs 

, 

X-4 

X-4 

X-4 

<—* 






~4 



•—4 


—4 

«— 4 

—4 

x-4 

x-4 

x-4 

—4 

4—4 

O' 

O' 

O' 

O' 

00 

00 

00 

r- 

r- 

r- 

M3 

M3 

M3 

'w ' 








^■x - 


X^XX 

x— x 

4— X 

x— x 

X— X 

X— X 

■^-x 

4*x 

4X 















<M 



O' 


o 


pn 

r- 

m 

O' 

—4 

o 

in 

in 

m 

ro 

sO 

-4 

O' 



M3 


PM 


M3 

m 

O' 

PM 

Tf 

m 

<M 


PO 

CO 


r- 


m 

m 

PM 

o 


sO 

o 

00 

m 

X 


ro 

© 

x-4 


ro 

M* 

CO 

o 

X4 

M3 

m 

in 

O 

PM 

— 


M 

— 


— 


— 

-* 


-* 

O' 

od 

Tf 


od 

ro 

-* 

od 

M* 

PM 

-* 

M3 

ro 

-> 

O' 

Tf 

-4 

M3 

PM 


__ _ 


__ _ 

. 

.- . 

J( _ L 













^_ L 















— < ^— < 

< 


— • — 1 

—4 

—4 

<—* 

«— 4 

—i 

—4 

1—4 

<— 4 

x^ 

x^ 


o 

O 

o 

X-4 

X-4. 

X-4 

X-4 

X-4 

X-4 

X-4 

X-4 

X-4 

X-4 

X— 4 

X-4 

X-4 


— < — • 

— H 

■— * 

— 1 -—l • 

•— * 

—I 

•—4 

—4 

—4 


—4 


x-4 

—4 

x-4 

x-4 

—4 

-4 

O' 

O' 

O' 

O' 

00 

00 

00 



r- 

M> 

M3 

M> 

•*— ' 

•*— ■* 



' — ■” 



s— X 

^>x 

X^X- 

x— x 

*— * 

x— ' 

■X— X 

x— X 

»- X 

4-X 

4— - 

4— X 

4— X 

4-X 

4— X 

4-X- 

4-X 

4-X 

4— ' 

4— ■ 



«— 

4— 

■W 

M* 

rr <M 

PM 

r- 

r- sc 

o 

o 

o 

00 

O' 

M* 

—4 

—4 

© 

m 

in 

m 

ro 

M3 

x4 ■ 

O' 

M* 


M> 

r- 

PM 


M3 

m 

O' 

PM 

M* 

m 

sO sO 

sO 

PM 

o m 

— * 

m 

CO 


in 

00 

00 

r- 

sO 

o 

CD 

m 

00 

r- 

ro 

o 

x4 


ro 

M* 

00 

O 


M> 

m 

m 

O 

(M 

CO Tf 

in 

nO 

r- r- 

oc 

od 

od 

O' 

O' 

o 

o 

O' 

ao 


-* 

ad 

CO 


od 


PM* 

— • 

M3 

rd 


o 

M" 


sd 

PM 

- 


_ 


___ 






























(M 

(M 


PM 

(M 


PM 


PM 

PM 

PM 

PM 


x-4 

—i 

X— ( 

© 

, 

X-4. 

X-— 

X-4 

x-4 

X-4 

X-4 

,-4 

X-4 

—4 

X-4 

X-4 

^4 

X-4 

X-4 

•—* 



— c 



—i 




r-t 

•—4 


x4 

x-4 

x4 

x4 

O' 

O' 

O' 

ao 

® 

OO 

r- 

1^ 


M3 

M> 

M3 

m 

m 

m 

s— ' 

s—' 


>««» 

^x” 





4-x 

4—” 

«X-X 

X— x 

X— X 

x 

x-x 

•w* 

4— X 

4-X 

4-X 

4»» 

4-X 

4-X 

4^ 


4-X 

4— 

4— 

4-X 1 

4-X 

4X 

4^ 




O' 

O' 


© 


m 


ro 

O' 

H 

o 

m 

m 

sO 

00 

ro 

M* 

O' 

M- 

n- 

M3 

f4 

PM 

-4 

sO 

m 

O' 

PM 

-«f 

m 

fO 


o 

00 


r- 


in 

cn 

PM 

o 

r- 

vO 

o 

CO 


m 

O' 

O 

o 



ro 

M* 

00 

o 

** 

M3 

m 

in 

o 

PM 

PM 


PM 

— 


-* 



-’ 



O' 

CO 

*r 

-* 


o 

PM 



PM 

-* 

M3 

rd 

-< 

O' 

Tf 


M3 

PM 





___ 

_ -_ 


_ . __ 


























— H 






— H 


*~4 

—4 

—4 

—4 

o 

© 

© 

o 



xs 

X-4 

X-4 

X-4 

X-4 

X-4 

X-4 

X-4 

XS 

X-4 

X-4 

X-4 

X-4 

X-4 

— H 

•-< 


—4 

■4 


—H 



—4 

<— 4 


—4 

x-4 

—4 

x4 

O' 

O' 

O' 

oo 

00 

CD 

00 

r- 


r* 

M3 

sO 

M3 

in 

m 

in 







»^X 


■XxX 



4— » 

4— X 

■—x 

•x— x 

4— X 


















— < 


O' 

O' 


o 


m 

r- 

ro 

O' 

■4 

o 

in 

in 

ro 

ro 

M3 

x-4 

O' 

M* 

r- 

M3 

r- 

PM 

-4 

M3 

m 

O' 

PM 


in 

ro 


o 

CO 


r- 


m 

m 

PM 

© 


vO 

© 

oo 

ro 

00 

f- 

ro 

© 



ro 


CD 

© 


M> 

in 

m 

o 

PM 

PM 


PM 

— 


-4 


— 



-I 

O' 

00 

*r 

-* 

ao 

ro 


oo 

M* 

PM 

— 

M3* 

rd 


O' 

* 


vd 

PM 

— 

___ 




— 














__ 














M* 

M* 


M* 

PM 


PM 


PM 

PM 

«M 

PM 



H 

■4 

o 

© 

© 

X— 

X-4 

X-4 

x-4 

X-4 

x-4 

X-4 

— 

— 

X-4 

X-4 

X-4 

X-4 


— H 


"N 



•~H 


■4 



•H 


—4 

—4 

x-4 

x4 

x4 

x^ 

O' 

O' 

O' 

o 

CO 

ao 

00 

r- 

r- 

r- 

si) 

M3 

M3 


s— ' 


s**-' 



■^-x 


>— X 

4— »■ 

-4^ 

x— X 

4«X 

X"^ 

«— X 

4-^ 

*— > r 

4-X 

4— X 

4-X» 

4—" 

4^* 

4-X 

4-X 

4-X 

—X- 

4-X 



4—X 

"4-X 

4-X 

M* 

*-« 


O' 

O' 


o 


m 


ro 

O' 

x-4 

© 

m 

m 

ro 

ro 

M> 

— 

O' 


r- 

M3 

n- 

PM 

4 

M3 

in 

O' 

PM 


in 

<o 


o 

00 


r- 


m 

m 

PM 

o 


sO 

© 

© 

ro 

00 

r- 

ro 

o 



PO 


CO 

© 


nO 

in 

in 

© 

PM 

PM 


PM 

— 


— 


-4 


-* 

-X 

O' 

od 

Tf 

— 

oo 

PO 


00 

M* 

PM 

— 

M> 

rd 

— 

O' 



M3 

PM 


































PM 

PM 


PM 

PM 


PM 


PM 

rM 

PM 

«M 

X4 

—4 

—4 

x4 

o 

O 

o 

X-4 

xX4 

X— 

X-4 

x-4 

X-4 

X-4 

x-4 

— 

X-4 

X^ 

X4 

xs 


— * 





<~4 




_4 


—4 

x-4 

—4 

x-4 

x4 

—4 


O' 

O' 

O' 

O' 

00 

00 

00 


r- 

r- 

o 

sO 

M) 

S"' 






^»x 


X 

4—X 


*4X 

X—» 

»— X 

»— X 

4< 

4^ 

•^x 

4— X 

4— X 













M* 



O' 

O' 


© 


m 

r*- 

ro 

O' 

—4 

o 

m 

m 

ro 

PO 

M3 

—4 

O' 


n- 

© 

1-4 

PM 

-4 

M3 

m 

O' 

PM 

rf 

m 

rO 


o 

00 


n- 


m 

m 

<M 

o 


© 

o 

ao 

rO 

00 

n- 

ro 

© 


-4 

ro 

M* 

00 

© 


vO 

in 

m 

O 

PM 

PM 


PM 



-d 


-I 

x-4 # 


-4* 

O' 

ao 

V 


oo 

ro 

- 

oo 


(M 

-—4 

M3 

rd 

x"4 

O' 



M3 

PM 



« 

■o .j 

3 c 

" C 

£ jc 


OO — -< N N ^ ^ 


tfisor~aoc'©iri©in©ir>OkA©irioirto-,w... — 

>o*>o r~ r~ co oo O' o-- o 


tn o it o >n o 


(Sj 

rt ' 

o E 


4 


.13(18) 5.83(19) 2.13(17) 2. 12 (18) 1.35(18) 2.01 (18) 1.84(18) 









and Figure 2. 1. 1, representing 0.025 and 0.25 ppm CO, were also treated. 

The optical thickness (total number of CO molecules per cm^ 
column of sight) is given for the mapping experiment in Figure 2. 1 . 2 and for 
the limb experiment in Figure 2. 1.3. 

The CO^ profile used is given in Table 2.1.2 and Figure 2. 1.4. 
The CO^ is based on a constant mixing ratio of 320 ppm'. The CH^ profile, 
taken from Cadle and Powers (1966) is shown in Table 2. 1.2 and Figure 
2. 1.4, and corresponds to a constant mixing ratio of 2 ppm. 

The three water profiles, corresponding to a dry^, normal 
(3) 

and a wet atmosphere, were taken from Gutnick, 1962; Oppel, 1963; and 

A * 

Linquist, 1965; respectively. These are given in Table 2.1.2 and Figure 
2.1.5. These are equivalent to approximately 0.2, 1.5, and 3 percipitable 
cm H^O. 

The planned satellite experiment may be best carried out in a 
polar orbit. Without noting reasons, advantages, arid disadvantages, this is 
excellent for the mapping experiment, but restricts the limb measurements 
to the polar regions. This is not a problem, however, since above the tropo- 
pause, it is to be expected that there is no' appreciable latitude effect on the 
concentrations. Keneshea (1971) has calculated the OH concentration at dif- 
ferent latitudes and found little variation. Since the chief CO removal mech- 
anism at these altitudes is, in all probability, CO + OH reaction, the rate of 
CO removal should be about constant with latitude. The small temperature 
effect on the rate constant should not cause significant variation of CO with 
latitude. 

2.2 Temperature and Pressure Models 

Four specific temperature profiles were used. These are 
shown in Table 2. 1.3 and Figure 2‘. 1.6. Profile 1 corresponds to a cold at- 
mosphere (AFCRL, 1966);. profile 2 corresponds to an average atmosphere 
(AFCRL, 1966); profile 3 corresponds to a hot atmosphere (AFCRL, 1966); 
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Figure 2.1.2 Optical Thickness as a Function of 
Angle from Vertical 




0 10 20 30 40 50 60 70 

ALTITUDE (KM) 


Figure 2. 1.3 Optical Thickness for Limb Transmission 

Experiment as a Function of Grazing Altitude 



TABLE 2.1.2 CONCENTRATION PROFILES" • 


Altitude 

C °2 

CH. 

(km) 

c 

4 

0 

7. 8 (15) 

5. 0 (13) 

1 

6. 8 (15) 

4. 1 (13) 

2 

5. 9 (15) 

3.7 (13) 

3 

5. 0 (15) 

3. 1 (13) 

4 

4.3 (15) 

2. 7 (13) 

5 

3.9 (15) 

2. 3 (13) 

6 

3. 3 (15) 

2. 1 (13) 

7 

3. 0 (15) 

1.8 (13) 

8 

2. 5 (15) 

1.4 (13) 

9 . 

2.2 (15) 

1.3 (13) 

10 

2.0 (15) 

1.2 (13) 

15 

1.0 (15) 

5. 8 (12) 

20 

5. 0 (14) 

2. 9 (12) 

25 

2. 4 (14) 

1. 3 (12) 

30 

1.2 (14) 

7. 0 (11) 

35 

6. 3 (13) 

4. 4 (11) 

40 

3.2 (13) 

1.8 (11) 

45 

1.7 (13) 

8. 5 (10) 

50 

8. 2 (12) 

4. 2 (10) 

55 

4. 1 (12) 

2. 0 (10) 

60 

2. 1 (12) 

1.0 (10) 

65 

1.2 (12) 

5.0 (9) 

70 

5.2 (11) 

2.4 (9) 

75 

1.2 (11) 

1,0 (9) 

80 

1.0 (11) 

5.0 (8) 

85 

3. 7 (10) 

2. 1 (8) 

90 

. 2. 0 (10) 

1.0 (8) 

95 

7.5 (9) 

4.0 (7) 

100 

3. 1 (9) 

1.0 (7) 


■«2° 

(D 

' (2) 

(3b 

3. 0 (16) 

2. 3 (17) 

4. 6 (17) 

1.2 (16) 

9.9 (16) 

1.9 (17) 

5. 1 (15) 

3.8 (16) 

7. 4 (16) 

1.2 (15) 

1.5 (16) 

2. 7 (16) 

4. 1 (14) 

3.9 (15) 

8. 2(15) 

9. 1 (13) 

5. 1 (14) 

1.7 (15) 

1.? (13) 

5.9 (13) 

:l. 5 (13) 

3.4 (12) 

8.4 (13) 

2. 8 (12) 
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NUMBER DENSITY CM 




TABLE 2. 1. 3 TEMPERATURE PROFILES (°K) 


Altitude 


(km) 


2 

3 

J 

0 

257. 3 

288. 2 

302. 6 

1 

259. 3 

281.7 

2.95.9 

2 

256. 1 

275. 2 

289. .4 

3 

252. 8 

268. 6 

284. 3 

4 

247. 8 

262. 2 

’ 277. 4 

5 

241. 0 

255. 7 

270. 6 

6 

234.2 

249. 2 

, 263. 8 

7 

227.4 

242. 7 

257. 1 

8 

220. 6 

236. 2 

250.4 

9 

217. 2 

229.7 . 

243. 7 

10 

217. 2 

223. 3 

, 237.0 

15 

217. 2 

216 . 6 

203.7 

20 

214,2 

216. 6 

206.7 

25 

211.2 - 

221.6 . 

221. 4 

30 

216.0 

226. 5 ' 

232. 3 

35 

222. 3 

236. 5 

243. 2 

40 

234. 6 

250. 4 

254. 0 

45 

247. 0 

264.2 

264. 8 

50 

259- 3 

270.6 

270. 2 

55 

259. 1 

265. 6 

263.4 

60 

250. 9 

255. 8 

253. 1 

65 

248. 4 

239. 3 

236.0 

70 

245. 4 

219. 7 

218.9 

75 

234. 7 

200. 2 

201. 8 

80 

223. 9 

180. 6 

184. 8 

85 

213. 1 

180. 6 

177. 1 

90 

202. 3 

180. 6 

177. 0 

95 

211.0 

195. 5 

184. 3 

100 

218. 5 

210. 0 

190.7 




profile 4 involves a low altitude inversion layer (up to 2 km), based on aver- 
age values for Vandenberg AFB, June 1970 (Vachon, 1971). It is important 
that these variations should have no large effect on the CO measurements . 
Any such effects require detailed accurate supplemental measurements and 
make the interpretation of the data much more complex and, for practical 
purposes, means that many fewer data can be interpreted. Analysis of such 
effects is covered in subsequent sections of this report. 

For calculations, described later, made to compare results of 
these studies with those of the Convair group (Ludwig, 1970), four other tem- 
perature models were used. These are given in Tables 2. 1.4 and Figure 
2.1.7 as models 5 through 11. 

2. 3 Other Atmospheric Characteristics 

In the development of the COPE experiment a number of atmo- 
spheric characteristics must be considered. Those of greatest interest are 
cloud cover, precipitable water, pressure, temperature, and inversion lay- 
ers. A report (Vachon, 1971) was prepared which presented a base of de- 
tailed information on these characteristics which can be extended or refined. 

Significant cloud cover denies the possibility of making the de- 
sired observation since radiation (of either the 2. 3 or 4. 6 p band of CO) does 
not penetrate the clouds. The amount of cloud cover has been found to vary 
mainly between 0. 3 and 0. 8 around the world. There is no part of the globe 
where there is not an appreciable time when the atmosphere is cloud free. 

An empirical formulation was developed to relate the probability of having a 
cloud-free line-of-sight at various look angles as a function of the frequency 
of occurrence of scattered, broken, and overcast cloud cover. - 

Data on precipitable water in the atmosphere shows variation 
between about 0. 5 and 5 cm. Although there is considerable variation there 
are predictable trends for any region and season. The knowledge of these 
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TABLE 2.1.4 TEMPERATURE MODELS FOR COMPARATIVE CALCULATIONS 
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trends now appears to be unnecessary as can be seen in later discussions 
(Section 7. 4). 

One effect which must be considered in remote sensing from a 

fast moving vehicle is the possibility of a Doppler shift. If there is a large 

velocity component in the direction of the observed radiation, there is a 
.... v s 

Doppler shift which is determined by Ay = V — where v is the velocity com- 

c s 

V 

ponent in the line-of-sight, or Ay =' y— cos 0 where 0 is the angle between 

the line-of-sight and the vehicle movement. In the limb experiment, the 

satellite may be traveling very nearly direct toward or away from the sun so 

5 -1 v s Ay 

that cos 0 « 1 , and v is of the order of 7 x 10 cm sec . Thus — is 

s c y 

5 -1 

of the order of 2. 3 x 10 . Thus at 2. 3 p, Ay = 0. 1 cm , and at 4. 6 p, Ay = 
0. 05 cm These shifts are of the order of a Doppler line-width, which is 
the width of the line at the upper altitudes, i.e. those of interest in the limb 
experiment. If the technique is spectral and employs a reference cell of any 
sort on board this is extremely important and probably means that the tech- 
nique is not suitable for the limb measurement. For the correlation inter- 
ferometric technique a Doppler shift essentially appears as a phase shift 
which can be readily taken into account. 


2.4 The Effect of Scintillations 

2. 4. 1 Introduction: - Scintillations in the light intensity of a beam 
are introduced by atmospheric turbulence. When the light from a small 
light source of constant intensity passes through a long path in the atmos- 
phere, the energy collected through a small aperture at the other end is not 
constant: rather it fluctuates with time. The effects of such fluctuations 
(usually called scintillations) on the performance of a ground-looking scan- 
ning Michelson interferometer are discussed here. The output of a scanning 
Michelson interferometer is normally a time varying signal. Therefore, any 
scintillations in the light beam, before it is collected by the Michelson inter- 
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ferometer, can introduce spurious fluctuations in the wave form of the inter- 
ferometer output. However, the effects of scintillations will be negligible 
if the amplitude of these scintillations turns out to be smaller than the error 
of digitalization in the interferometer output wave form. First, we shall 
qualitatively explain why the scintillation effects will not be significant with 
the earth-looking interferometer. Next, we shall quantitatively determine 
the extent of the effect of scintillation on the interferometer data. 

2.4.2 Qualitative Explanation: - It is generally said that to the eye, 
stars, which have an angular subtense of a hundredth arc sec or less, ap- 

J . , 

pear to twinkle. On the other hand, objects of larger angular subtense - 
such as planets, which have angular subtense from several arc secs to 30- 
40 arc sec - do not. Consider Figure 2.4. 1 and suppose that with a receiver 
(a telescope or the pupil of the eye) of diameter, d, on the ground we observe 
an object outside the earth's atmosphere, subtending an angle 0. For the 
purpose of scintillations, we take the height of the "top" of the atmosphere 
as 100 km. (The exact value of the height of the "top" is not important for 
the present discussion. ) For an eye-pupil, AD = d =» 6 mm so that, if the 
object is a star (9 =* 0. 01 arc sec), BC 10 mm. On the other hand, if the 
object is a planet (0 =“ 20 arc sec), BC ^ 10 , 000 mm (10 meter). Thus , in 
the case of the planet, it is the increased number of atmospheric turbulence 
eddies perpendicular to the line-of- sight and within the light cone ABCD 
that are responsible for reducing the scintillations and stabilizing the inten- 
sity. 

If we view a star (9 = 0.01 arc sec) with a 100 cm aperture 
telescope, then AD “ 100 cm and BC =** 100 cm, and again the number of 
turbulence eddies perpendicular to the line-of-sight will have increased 
and, consequently, scintillation will be decreased. This is, of course, an 
observed fact: scintillations become negligible when photo-electric measure- 
ments of the brightness of a star are made with a large aperture telescope. 
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The situation of an earth-looking scanning Michelson interfer- 
ometer is shown in Figure 2.4. l b. Consider an interferometer receiver of 

diameter, AD, 3-inches. Then for a wavelength' of 2 jj, and resolution, &v , 

- 1 2 - 3 

of 0.5 cm , the throughput is 0.006. (cm ster), which gives 0 =- 6 x 10 

radians (= 1200 arc sec) and therefore A / D , ’= 600 meters and BC = 1200 
meters. Even if the height ofthe satellite or that of the "top" of the atmos- 
phere is different from the above values, it is quite evident that, due to the 
large value of 0, dimensions A / D / and BC of Figure 2. 4. 1 b will be very much 
larger than the corresponding dimensions for the case of a planet or a large 
aperture telescope considered above. Consequently there will be a very 
large number of turbulent eddies in the direction perpendicular to the line- 
of-sight and the effects of the scintillation will be negligible. 

2.4.3 Criterion for "Significant" Amount of Scintillations: - The in- 
strument used for this purpose will be a Michelson interferometer in which 

* i . 

the beam of light falling on the instrument is split with a beam splitter into 
two mutually coherent beams which are then made to interfere with each 
other to produce Haidihger fringes. Light from the central region of this 
fringe pattern is collected by a photo detector and the output of the instru- 
ment is an electrical signal. The spectral nature of the radiations is ana- 
lyzed by varying the optical phase difference between the two beams and re- 
cording the fluctuating electrical output; that is, the interferogram . The 
fluctuations in this recorded signal contain the information on the spectral 
properties of the radiations. However, the light beam falling on the instru- 
ment fluctuates in intensity. The observed interferogram will have spurious 
fluctuations resulting from the modulation of the interferogram by light in- 
tensity fluctuations. In Section 2.4. 5, therefore, we shall investigate the 
influence of incident light beam intensity fluctuations upon the energy spec- 
trum obtained with the interferometer. ; • 
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2.4.4 Additive and Multiplicative Noise: - Before we make quantita- 
tive calculations of scintillations, we must point out the distinction between 
multiplicative and additive noise. 

2. 4. 4. 1 Additive Noise: - Let the amplitude of the noise be 
independent of the signal. This will be the case with digitization noise. 

For a six-bit word, the digitization noise amplitude will be: 1/2^ = .016 =- 

1.6%. Under such a situation, the greater the bandwidth -- (i.e. the greater 
the phase difference for scanning)-- the spikier the reconstructed spectrum. 
Therefore, to reduce digitization noise the path difference for scanning 
should not be greater than that dictated by the needs of spectral resolution. 
Conversely, even a small additive error can produce large errors in the re- 
constructed spectrum if the bandwidth of the scan is large. Of course the 
greater the amplitude of the additive noise (e. g. digitization noise) the 
greater will be the amplitude of the spurious spikes in the reconstructed 
spectrum. 

2. 4. 4. 2 Multiplicative Noise : - The noise is multiplicative 
when it is dependent upon the signal level, i.e. when the amplitude of the 
noise is a certain fraction of the instantaneous signal. For example, let the 
true interferogram be F^uj) giving a true spectrum, f Q (\). 

Then, _ 

f(\) =7 [f(w)] 

where represents the usual linear operation of Fourier transform of a 
real signal. Let the multiplicative noise be a fraction, €(uj), of the signal, 
F(x), where the value of fraction €(m) is a random variable. Therefore the 
observed interferogram can be represented as 

F (a) - F (u>) + € ( uu ) F ( uj ) = f 1 + € ( u) ) ~| F (cju) 
obs o o L j o 
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Therefore the observed reconstructed spectrum will be 

£ obs< X) = 7[ F „bs<»>] =7[ f 0 (»)] *?[«(•>]*•„<»>] 


= f 


> (X) +^[€(0)) F o (w)~ 


The error in the reconstructed spectrum therefore will be 

A f (x.) = f Qhs (V - £ C (X) = ?[€(«!)) F Q (U))]. 

But, since €(uu) is a random variable, 


?[€(.) F 0 t.)] - ! € max l'f[ F 0 (“)] = 


= l€ •£ (X). 
max o 


where € is the maximum likely fractional error in the interferogram. 
Consequently, the maximum likely error in the reconstruction spectrum 
will be 


Af(X) 

f o (x) < 


€ 


max 


The effect of light intensity scintillations will be a 
multiplicative error on the observed interferogram and the reconstructed 
spectrum. Therefore, for example, if the scintillation amplitude is 1%, the 
error in the observed intensities in the reconstructed spectrum will not be 
more than 1%. 

2. 4. 5 Calculation of the Scintillation Effect: - The effect of an atmo- 
spheric refractive index inhomogenity on the light propagating through tur- 
bulent atmosphere and the consequent scintillation effects have been investi- 
gated by various authors, as noted in the following discussion. We shall 
make use of these investigations in calculating the magnitude of the scintilla- 
tion expected from an earth-looking satellite. 
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. An exact mathematical description of atmospheric turbulence 
is not possible. However, with some reasonable assumptions , the problem 
becomes,, tractable. . In the work of earlier investigators the atmospheric 
turbulence is assumed to be homogeneous and isotropic and the spectral en- 
ergy distribution in the turbulence is assumed to be given by the Kolmogoroff 
theory. Let C (p) be the convariance of the log-amplitude. between two mea:- 
surement points separated by a distance, p,. Then, scintillations in a laser 
beam traveling .through a uniform atmosphere have been calculated by Fried 
and Seidman (1967) showing that for p = o, C (O) can be written as a product 
of two factors: 


C (O) = C (O)-f (Q) 

\a \a O 


Eq. 2. 1 


■ S, 


where the first factor, C (O), the value of C (O) associated with a point 
source/ is given by: 


C- S (0) = : constant- k ?/ 6 z I1/6 C 2 
t hi 


Eq. 2 . 2 


2 7T 2 

(k =-^— ') V = wavelength of light, z = path length, and C is the refractive 

index structure constant). The second factor is a function of the parameter 
2 . 

a (where fl = k« /z, and a is the linear source size). 1 The value of f(A) is 
o 

given -by Fried and Seidman (1967, Equation 4. 12). • 

• For . the case of an earth-looking satellite, the atmosphere in 

2 , 
the line-of-sight is not uniform and therefore C" is not a constant. Based 

N 

on the data presented' by Hufnagel ( 1 964) it has been : shown (Fried and Cloud, 

2 

1966) that the value of at an altitude, h, is given by: 


C. 2 = 4.2 x 10 ^ h exp(-h/h ) 
N o 3 


Eq; 2. 3 


where the scale height, h , is 3. 2 x 10 meters. Fried (1967b) has made 

° 2 
an evaluation of C (O) and f (A) for the case of altitude dependent given 

l n N 

by Eq. 2. 3. Following him we can derive: 
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C (O) = 0. 122 sec 
l 


il/l6' rt r2_ i7/6 
X _ 


0 


where the wavelength, X , is in microns. For \ = 2.3 microns and 0 = 0; 

S ■ 

C (O) = 0. 122, i.e. for a point source on the ground C (O) = 0. 122. 

As As 

To calculate the magnitude of scintillations we first calculate 
the fractional intensity variance C^(O) from the relation (Fried, 1967a). 


Cj(O) = exp^C^O)] - 1 


Eq. 2/4 


thus for C^(O) = 0. 122 we get: 

C (O) = 0.629. 

Noting that the above value is for fractional intensity and comparing Eq. 4 
of Fried (1957a) with Eq. 13-13 of Tartskii (1961), we obtain the RMS devi- 
ation of the scintillation = /0. 629 - . 79 °* 80%. 

The above RMS value of scintillation has been ’obtained for a 
point source on the ground. However, the Michelson interferometer on the 

satellite will have a field-of- view of, say, four degrees (which gives 3. 8 x 

-3 

10 steradians for a circular field-of- view) . This increase in the field-of- 
view will decrease the magnitude of scintillations. 

S 

Fried has evaluated the function f (Q) 1 = C„(0)/C„ (O) that 

n -l l 

gives the reduction in C (O) with increasing value of the source size param- 

^ S 

eter, fj. This dependence of C AO) I C (O) on O is shown in Figure 2.4.2. 

As As 

For our case of an earth- viewing instrument in a satellite we put fl in a 
more useful form: 


o = 


uuh sec 0 

Th" ‘ 


where 
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Ci(0)/Ci(0) Upon the Normalized So\ 
Observations from Space of a Ground 










<JU 


solid angle of the field-of- view in steradians, 

h = height of the satellite orbit about the ground, 

3 

h^ = scale height of turbulence (3.2 x 10 meters), 

0 = zenith angle of the instrument with respect to the 

source point on the ground, and 

X = wavelength. 


To calculate the upper bound on the magnitude of scintillation we take: 



0) = 

10 k steradians (for two arc minute field-of- view) 


h 

5 

10 meters 


h 

o 

3 

3. 2 x 10 meters 


0 

0 


X 

2.3x10 ^ meters 

Thus, j 

?. =- 10 6 . 



This is a very large value for and from Figure 2. 4. 2 it is quite clear that 

S 6 - 6 

the value of C (0)/C (O) for Q 10 is essentially zero (less than 10 ). 

c ^ ^ 

Since C (O) = 0. 122 as calculated above, C (O) < 0. 122 x 10 therefore 

'C / v 

C^O) = exp|^4C^(0)] - 1 <0.5xl0" 6 , 

and the RMS value of scintillation amplitude will be less than 0. 07%. There- 
fore, we conclude that, for the above described technique of measuring car- 
bon monoxide with a Michelson interferometer on an orbiting satellite, the 
effects of atmospheric scintillations will be less than 0. 1% in the recon- 
structed spectrum. 

It should be noted that for calculating the upper bound for scin- 
tillation, we have already chosen the greatly exaggerated values for the 
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parameters. For example, vye have taken uj = 10 radians while the actual 

-3 

value will be 4 x 1 0 radians or greater. Similarly the height of the satel- 

5 

lite will be. much greater than 10 meters chosen here. For the situation in 
question, therefore,, the conclusion is that the effect of scintillation will be 
negligible. 


28 



3. 


SPECTRAL CONSIDERATIONS 


3. 1 The CO Infrared Spectrum 

The infrared spectrum has been discussed in detail elsewhere 
(Bortner and Kummler, 1971). The two bands under consideration herein 
are the first overtone band centered about 2. 35 p, and the fundamental band 
centered about 4. 67 p. The positions and the energies of the various rota- 
tional lines of these bands is given in Table 3. 1. 1 for the first overtone band 
and in Table 3. 1.2 for the fundamental band. The line strengths of the lines 
of each band are given in Table 3. 1. 3. It can be seen that the fundamental 
band is stronger than the overtone by a factor of the order of 100. The fun- 
damental is less interfered with by absorption due to other atmospheric con- 
stituents. The overtone band is capable of much simpler interpretation, be- 
ing much less affected by atmospheric temperatures, ground temperatures, 
and ground emissivity. As seen in Figure 3.2. 1 the radiation seen in this 
experiment in the overtone band region would consist almost exclusively of 
reflected sunlight while that radiation for the fundamental band region would 
be mainly earthshine with some effect of reflected sunlight and serious atmo- 
spheric emission effects. These factors will be discussed later. 

3. 2 Intensity of Reflected Solar Radiation and Earthshine 

3.2. 1 Solar Reflection : - The solar radiation reaching the earth and 

being reflected can be calculated for any wavelength as follows. From black- 

body radiation tables (Bowen, 1963) the radiant intensity, J, for a 5900 K 

1 -2 -1 -1 

blackbody is calculated to be 9. 3 x 10 watt cm sr |j, at 2. 3 ^ and 8. 0 x 
0 

10 at 4. 6 (j,. The average irradiance, H, incident on the earth's atmosphere 
(and incident on the earth's surface in the absence of any atmospheric atten- 
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TABLE 3. 1. 1 ENERGY & WAVELENGTH OF VIBRATION- 
ROTATION LINES OF CO (2-0) 


R BRANCH 2-0 P BRANCH 

' // /■//'.. 


J 

CM 

-1 

X(|i) 

f r r 

J -J CM 

-1 


X(li) 

1 - 0 

42 63. 

831 

2. 34531 

0 - i 

42 56 

21 1 

2. 

34951 

2- I 

42 67. 

53 6 

2.34327 

1- 2 

42 52 

29 6 

2. 

351 67 

3.7- a 

42 7 1 . 

1 70 

2.34128 

2- 3 

4248 

31 1 

2. 

3 5388 

4-3 

4274. 

734 

2. 33933 

3- 4 

4244 

257 

2. 

35612 

5- 4 

4278* 

227 

2.33742 

4- 5 

4240 

133 

2. 

3 5842 

6- 5 

4281 • 

649 

2.33555 

5- 6 

4235 

940 

2. 

3 6075 

7- 6 

4285* 

000 

2.33372 

6- 7 

4231 

678 

2. 

3631 3 

8- 7 

4288. 

283 

2.33194 

7- 8 

422 7 

347 

2. 

36555 

9- 8 

4291 . 

488 

2.33019 

8- 9 

4222 

947 

2. 

36801 

10- 9 

4294. 

625 

2.32849 

9-10 

4218 

479 

2. 

3 70 52 

11-10 

4297. 

690 

2.32683 

13-11 

4213 

942 

2. 

3730 7 

12-11 

4330. 

684 

2.32 521 

11-12 

4209 

337 

2. 

37567 

1 3-12 

4333. 

60 5 

2. 32363 

12-13 

4204 

664 

2. 

37831 

14-13 

433 6. 

4 54 

2.32210 

13-14 

4199 

923 

2. 

38100 

15-14 

4339. 

231 

2.320 60 

14-15 

419 5 

1 1 4 

2. 

38373 

16-15 

4311. 

935 

2.31914 

15-16 

4190 

237 

2. 

38650 

17-16 

431 4. 

567 

2.31 773 

1 6- 17 

4185 

293 

2. 

38932 

18-17 

431 7. 

125 

2« 31 636 

17-18 

41 80 

282 

2. 

39218 

19-18 

4319. 

61 1 

2.31 502 

18-19 

41 75 

203 

2. 

39 509 

23-19 

4322. 

324 

2.31373 

1 9- 20 

41 70 

0 58 

2. 

3980 5 

21-23 

4324. 

3 63 

2.31248 

20-21 

41 64 

845 

2. 

40105 

22-21 

4326. 

628 

2. 31 127 

21-22 

41 59 

567 

2. 

40410 

23-22 

4328- 

821 

2.31010 

22-23 

41 54 

221 

2. 

40719 

24-23 

4333. 

939 

2. 30897 

23-24 

4148 

810 

2. 

41033 

2 5-24 

4332. 

983 

2. 30788 

24-25 

4143 

332 

2. 

41 352 

2 6-25 

4334. 

953 

2. 30 68 3 

25-26 

41 37 

788 

2. 

41 675 

27-26 

4336. 

8 49 

2.30 582 

26-27 

41 32 

179 

2. 

42003 

28-27 

4338. 

670 

2. 33 48 5 

27-28 

4126 

504 

2. 

42336 

29-28 

4340. 

41 7 

2. 30 39 3 

28-29 

4120 

763 

2. 

42674 

33-29 

4342. 

089 

2. 30334 

29-30 

4114 

957 

2. 

4301 6 

31-30 

4343. 

686 

2. 33219 

33-31 

4109 

086 

2. 

43363 

32-31 

4345. 

208 

2. 30139 

31-32 

4103 

1 50 

2. 

4371 5 

33-32 

4346. 

654 

2* 30362 

32-33 

4097 

1 49 

2. 

440 72 

34-33 

4348. 

326 

2.29989 

33-34 

4091 

083 

2. 

44434 

3 5-34 

4349. 

321 

2.29921 

34-35 

4084 

9 53 

2. 

44801 

36-35 

43 50. 

541 

2.29856 

3 5-36 

4078 

7 59 

2. 

451 73 

37-36 

4351 . 

685 

2.29796 

36-37 

40 72 

501 

2. 

45549 

38-37 

4352. 

7 53 

2.29740 

3 7-38 

4066 

1 78 

2. 

45931 

39-38 

43 53. 

745 

2.29687 

38-39 

40 59 

792 

2. 

46318 

43-39 

4354. 

660 

2.29639 

39-40 

40 53 

342 

2. 

46710 



TABLE 3. 1. 2 ENERGY & WAVELENGTH OF VIBRATION- 
. ROTATION LINES OF CO (1-0) 



R BRANCH 


I 

- 

P BRANCH 



/ // 




/ // 

1 



J -J 

cm' 1 . 


Ml*) - 

J -J 

CM 

■ I 


Mm.) . 

1-0 

21 47.083 

4. 

65748 

0- 1 • 

2139. 

428 

4. 

6741 5 

2- 1 

21 50.8 58 

4. 

64931 

1-2 

2135. 

548 

4. 

68264 

3- 2 

21 54. 598 

4. 

64124 

2- 3 

2131. 

634 

4. 

69124 

4- 3 1 

21 58.301 ’ 

4» 

63327 

3- 4 

2127. 

684 

4. 

69994 

5-4 

21 61 .9 70 

4. 

62 541 

4- 5 

2123. 

70 1 : 

4. 

70876 

6-5 

2165. 632 

4. 

61 765 

5- 6 

2119. 

683 

4. 

71 7 69 

7- 6 

21 69. 1 98 " 

4. 

61300 

6- 7 

2 1 l 5* 

631 

4. 

72672 

8-' 7 

2 1 72 • 7 S8 

4. 

60244 

7- 8 

2111. 

545 

4. 

73587 

9- 8 

21 76.282 

4. 

59499 

8-9 

2107. 

42 6 

4. 

74513 

10- 9 

21 79. 770 

4. 

58764 

9-10 

2103. 

273 

4. 

7 5449 

11-1,3 

2183.220 

4. 

580 39 

13-11 

2099. 

086 

4. 

76398 

12-11 

2186. 63 4 

4. 

57324 

11-12 

2094. 

867 

4. 

77357 

13-12 

2190.01 1 

4. 

56619 

12- 13 

2090. 

61 4 

4. 

78328 

1 4-1.3 

21 93.3 51 

4* 

55923 

13-14 

2086. 

328 

4« 

7931 1 

1 5-14 

2196. 653 

4- 

55238 

14-15 

2082. 

010 

'4* 

8030 5 

1 6-1 5 

2199.918 

4. 

54562 

15-16 

20 77. 

659 

4 . 

81311 

17-16 

2233-1 45 . 

4. 

53897 

16-17 

2073. 

275 

4. 

82329 

18-17 

2206.334 

4. 

53243 

17-18 

20 68. 

8 60 

4. 

83358 

19-18 

2239.486 

4. 

52 59 4 

18-19 

20 64. 

412 . 

4* 

84399 

20-19 

2212. 599 

4. 

51 957 

19-20 

2059. 

932 

4. 

854 53 

21-23 

221 5-674 

4 « 

51 333 

23-21 

23 55. 

421 

4 . 

86518 

22-21 

221 8. 713 

4. 

53 712 

21-22 

23 50. 

878 " 

4 . 

87596 

23-22 

2221 . 738 

4. 

53104 

22-23 

2046* 

30 3 

4 . 

88 68 6 

24-23 

2224- 667 

4. 

49 53 5 

23-24 

2041 . 

69 7 

4 . 

89789 

2 5-24 

2227.587 

4. 

4891 6 

24-25 

2337. 

063 

4 . 

93903 

2 6-25 

2233.469 

4 . 

4833 6 

25-26 

2032. 

393 

4 . 

920 31 

2 7-26 

2233*310 

4 . 

47766 

2 6-27 

2327. 

69 4. 

4 . 

931 71 

28-27. 

2236. 1 1 3 

4. 

4723 5 

27-28 

2322. 

965 

4. 

94324 

29-28 

2238.875 r 

4. 

46653 

28-29 

2318. 

20 5 

4. 

95493 

3'?-29 

2241.598 

4. 

461 10 

29-30 

2013. 

41 5 

4. 

9 6669 

3 1 - .33 

2244.281 

4. 

4 5 577 

33-31 

2308. 

59 5 

4 . 

97863 

32-31 

2246.924 

4* 

4 53 53 

31-32 

2333. 

745 

4. 

990 65 

33-32 

22 49. 52 7 

4. 

44538 

32-33 

1 998. 

865 

5. 

00284 

34-33 

2252.389 

4. 

440 32 

33-34 

1 993. 

956 

5. 

01 51 6 

35-34 

22 54. 61 1 

4. 

43 53 6 

34-35 

1 989. 

017 

5. 

02761 

3 6-35 

2257.092 

4. 

43048 

3 5-36 

1 984. 

049 

5. 

0 4320 

37-36 

2259. 532 

4 . 

42 573 

36-37 

1 979. 

051 

5. 

3 5293 

38-37 

2261 .931 

4 . 

42100 

37-38 

1974. 

025 

5* 

36579 

39-38 

2264.289 

4. 

41 643 

38-39 

1968* 

9 73 

5. 

0 7880 

43-39 

2266. 635 

4. 

41 1 88 

39-40 

1963. 

886 

5. 

39194 



TABLE 3. 1. 3 CO LINE STRENGTHS (cm' 2 atm -1 ) 


FUNDAMENTAL OVERTONE 

(Benedict, 1962) (Kostkowski, 1961) 



R 

P 

R 

P 

0 

1.970 


. 0163 


1 

3. 878 

1.927 

. 0321 

. 0158 

2 

5. 643 

3. 707 

. 0467 

. 0303 

3 

7. TOO 

5. 251 

. 0592 

. 0428 

4 

8. 261 

6. 491 

. 0691 

. 0527 

5 

9.059 

7. 386 

. 0760 

. 0598 

6 

9. 482 

7. 920 

. 0797 

. 0638 

7 

9.544 

8. 105 

. 0805 

. 0651 

8 

9. 286 

7.980 

. 0785 

. 0638 

9 

8.760 

7. 591 

. 0742 

. 0604 

10 

8. 032 

7. 002 

. 0682 

. 0555 

1 1 

7. 170 

6. 278 

. 0610 

. 0495 

12 

6. 240 

5. 480 

. 0532 

. 0430 

13 

5. 301 

4. 665 

. 0453 

. 0365 

14 

4.400 

3. 876 

. 0376 

. 0301 

15 

3. 568 

3. 145 

. 0306 

. 0243 

16 

2.830 

2. 494 

. 0243 

. 0192 

17 

2. 196 

1. 934 

. 0189 

. 0148 

18 

1. 668 

1. 486 

. 0144 

. 0112 

19 

1.241 

1. 091 

. 0107 

. 00828 

20 

0. 904 

0. 793 

. 00779 

. 00597 

21 

0.645 

0. 566 

.00557 

. 00423 

22 

0. 451 

0. 395 

. 00390 

. 00294 

23 

0. 309 

0. 290 

. 00268 

. 00200 

24 

0. 208 

0. 181 

. 00180 

. 00133 

25 

0. 137 

0. 119 

.00119 

. 000872 

26 

0.0884 

0. 0766 

. 000765 

. 000556 

27 

0. 0589 

0. 0483 

. 000486 

. 000350 

28 

0. 0349 

0. 0300 

. 000302 

. 000215 

29 

0. 0212 

0. 0182 

. 000184 

. 000130 

30 

0. 0127 

0. 0108 

. 000110 

. 0000770 

31 

0. 00743 

0. 00634 

. 0000642 

. 0000446 

32 

0. 00427 

0. 00364 

. 0000371 

. 0000254 

33 

0. 00241 

0. 00205 

. 0000209 

. 0000142 

34 

0. 00134 

0. 00113 

.0000116 

. 0000078 
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-3 -Z -1 

uation is (Thekaekara, 1971) 6. 8 x 10 watts cm y at 2. 3 y and 5. 5 x 

-4 -2-1 

10 watts cm y at 4. 6 y. 

The reflected solar radiation, N, from the earth's surface is 

the 


-2 -1 - 1 dH 

N (watts cm ster y ) = 

V 

-3 -2 -l -i 

The values of N are thus 2,16x10 p watt cm ster y at 2. 3 y and 
1. 68 x 10 p watt cm ster y at 4. 6 y. The intensities (in watt cm 
sr y are shown in Figure 3.2. 1 for various reflectivities. 


3.2.2 Earthshine: - From blackbody radiation tables (Bowen, 1963) 
the radiant intensity is given for a 300 °K blackbody and by the ratio 


R 1. 439 x 10 4 X x 300 

I 6 — 1 

R 300 1. 439 x 10 4 x T , 

e - 1 

-2 - 1 - 1 

the intensities (watts cm sr y ) for 273.2 and 288 K blackbodies were 
calculated giving the following, as shown in Figure 3.2. 1. 


T 

g 

2. 3 y 

4. 6 y 

273. 2 

2.0 x 10" 8 

6. 15 x 10 5 

288 

7. 54 x 10" 8 

-4 

1.11x10 

300 

1 . 60 x 1 0 

-4 

1.71x10 


3.2. 3 Ratio of Reflected Solar Radiation and Earthshine Intensities: - 
The ratios of reflected sunlight to earthshine at 2. 3 y are thus 
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and rl 


p 

T : 
g 

273. 2K 

288 K 

300 K 

0.02 


2 x 10 3 

5. 4 x 10 2 

2. 5 x 10 2 

0 . 1 


4 

1 x 10 

3 

2.7x10 

1.2 x 10 3 

J. 0 


1 x 10 5 ! 

4 

2. 7 x 10 

1.2 x 10 4 

3 of earthshine to reflected 

sunlight at 4. 

6 ^ are 

R 

T : 
g 

273. 2K 

288 K 

300 K 

0.02 


1. 78 x 10 1 

3.20 x IQ 1 

4. 95 x IQ 3 

0.1 


3. 56 x 10° 

6. 40 x 10° 

9. 90 x 10° 

* 

1.0 


3.56 x 10' 1 

-1 

6.40 x 10 

-1 

9.90 x 10 
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4. RADIATIVE THEORY 


4. 1 Theory 


To study the feasibility of the CO density measurements and the 
effect of various atmospheric parameters on such measurements, calcula- 
tions have been made using the radiative transfer theory described in this 
i 

section. First, the physical model and solution thereof are presented and 
then the physical parameters and their scaling relationship with altitude are 
defined. 

The model includes absorption and emission of photons along a 
ray from the source to the detector. The two geometries considered are di- 
rect solar observation - limb experiment, and absorption of the earth’s al- 
bedo - mapping experiment. Provisions are made for adding earthshine to 
reflection of sunlight from the earth’s surface. However, phenomena such 
as scattering resonance fluorescence and radiation trapping have been ne- 
glected. Boltzmann populations as described by the local atmospheric tem- 
perature have been assumed for the upper and lower states of all transitions 
and the line profiles are given by the Voigt function at all altitudes. Doppler 
shifts as caused by atmospheric winds and by relative motion between the 
atmosphere and the satellite have been neglected. 

The monochromatic radiative transfer equation is given 
(Chandrasekhar, I960) by: 


Where: 


dl 

v_ 

dS 


€ 

V 


I a 


v v 


Eq. 4.1 


I = spectral intensity (watts /cm -ster-cm ) 
S = distance along ray (cm) 
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€ = spectral emission coefficient (watts /cm -ster-cm ), 

a function of S 

CL ^ = spectral absorption coefficient (cm ^), a function of S 

Formally, 

1^ = C exp(-Ja dS 7 ) + exp(- Ja^dS 7 ) J exp£ Ja ^dS ,7 J dS 7 Eq. 4.2 

i 

For the limb boundary condition: 


I (z) = I exp(- r Z a dS") + exp(- P a dS u ) [* <■ exp(P S a. dS") dS 7 
V O , V J V ^ V J V V 


or 


Eq. 4. 3 


I (z) = I exp(- f a dS ) + F c exp(- f a dS ) dS 
V o, v J V Ju >’ / P 

o o S 


Eq. 4.4 


The geometry for this situation is shown in Figure 4.1,1. 

For the mapping experiment, two conditions are required. The 
solar energy is attenuated as it penetrates to the earth's surface. At this 
point, it is reflected and the total intensity increased by earthshine. The 
resulting flux is then further attenuated by the atmosphere until the satellite 
is reached. It is pointed out that atmospheric emission which is reflected 
by the earth into the field-of- view of the detector is negligible and therefore 
not included in the model. The geometry of the mapping experiment is given 
in Figure 4. 1,2 and is described mathematically by Equation 4. 5. 

z z z 

l U = % [ I o,„ exp( -J' ladS " } + J ?***&] Eq. 4.5 

o o S 


z 

+ I e exp(- J c^dx 77 ) 
o 


+ 


z z 

f ex P(-J a dx") dx 7 
1/ /l/ 
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Where 


2 - 1 

I = earths hine (watts / cm -ster-cm ) 

e 

p = earth reflection coefficient 
V 


Equations 4.4 and 4. 5 have been numerically integrated along the rays indi- 
cated in Figures 4. 1. 1 and 4.1.2 using Simpson's rule. The properties 

and a are functions of altitude which must be evaluated prior to integration. 
V 

This is considered below. 

On a microscopic basis, for a single rotational line (Kulander, 

1964), 


n hey 
u 

e - cep = — ; A cp 

v v 4tt ul T i/ 


Eq. 4. 6 


a = (a. - a „) cp = (n B, - n B „ )• ~ 7 ~ cp Eq. 4.7 

y A S* ^y l -tu u ul 4u ^y 


where 


9 

y 

n , n 
u L 

hey 

*A 


nt 


B 


-tu 


Voigt profile (1/cm - ^) 

3 

upper and lower state number densities (particles / cm ) 

energy per transition (watt sec/transition) 

absorption coefficient (cm */cm) 

stimulated emission coefficient (cm */cm) 

Einstein coefficient for spontaneous emission 

T ransitions 
Sec. Particle 

Einstein coefficient for absorption 

Transitions Cm-Ster 

Particle Watt Sec. 
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B . = Einstein coefficient for stimulated emission 

ut 

( T ransitions Cm-Ster \ 

Particle Watt Sec. 1 


It is noted that the definition of the Einstein coefficients is arbitrary but that 
they are related by Kirchoff's Law for equilibrium at temperature, T. 


or 


B (T) a. 

V V 


Eq. 4.8 


Thus , 


2 h 2 3 2 1 

= 'hct>/kT (watts / cm -ster-.cm ) 

e -1 


n A = (n B - n B .) B (T) 
u ut t tu u ut ir 


Eq. 4.9 


A /B 
ut/ ut 


2 V, 2 3 

2 he v 


( \ B -tu/ n u B u^ ’ ex P(hct//kT)-l 


Eq. 4. 10 


But, at equilibrium the Boltzmann distribution is: 


n 


l 


n 

u 


g ; 

exp(hcp/kT) 

g„ 


and by comparison 


g B = g B 
l -tu u u-t 


A - B * 2hc 2 y 3 
u-t ut 


Eq. 4. 1 1 


Eq. 4. 12 
Eq. 4. 13 


The Voigt profile, tp^, is defined in terms of the Doppler and Lorentz half 
widths, i.e. half the line width at half the peak intensity (Armstrong, 1967). 
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1 / ln2 


v a D l 17 


K(x,y) 


Eq. 4. 14 


K{x, y) = f 


Y / exp(-t ) 

IT 2 X / 

J y + ( x - 


L n ?x 1/2 

y = — ( ln2 ) 

"d 


■'o' „ ,.1/2 

— — (1H2) 

D- 


a = 
D 


Doppler half width (cm ) 


2R(ln2) T 1/2 , r „ ,„-7 T 1/2 -1 

= 3.58x10 i/ — cm 

c M o M 


a L 


= Lorentz half width (cm ) 


= r/2irc 


1,4 PCT 12 N o(kTRt) 


, 2 
2n c 


gas collision frequency (sec ) 

2 

pressure (dynes/cm ) 


optical collision cross-section (cm ) 


reduced mass of the colliding species (gm/gm mol) 
gas constant (erg/gm mol °K) 


temperature ( K) 
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N q = Avagadro's number (particles /gm mol) 
c = velocity of light (cm/sec) 

The Doppler and Lorentz widths are functions of temperature and pressure 
and thus are easily obtained at any point along the ray from the atmospheric 
properties. The integral, K(x,y), is then calculated using a subroutine de- 
veloped by Armstrong (1967). It should be noted that 

J K(x, y) dx = tt 1/2 


The emission coefficient number density product necessary for Equations 
4.6 and 4. 7 is most easily obtained by scaling experimentally derived line 
strengths which are defined as (Penner, 1959, p.20)S 


Where 



1 _ 

„ 2 

8ncu 





1 - exp 


hcv 

kT 



Eq. 4.20 


S = line strength (cm -atm ) 

V = wavenumber (cm *) 

P = pressure (atm) 

N„ = lower state population for conditions where S is 

l 3 

measured (particles / cm ) 
g , g = statistical weights (2 J + 1) 

li. 

A = Einstein coefficient 

uf. 

Transitions 
Sec. Particle 
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Using equation 4. 11 and rearranging 


n 

u 



„ 2 - 

8 ircy P s 

[exp(hci//kT)-l] la 


Where 


Eq. .4.21 


n = upper state population for conditions where S is measured, 
u ^ £u 

(particle /cm ) 


But 


_ ft > Q exp (-hcp/kT) Eq. 4.22 

n n Q exp (-hcp/kT) 

u t 

total number density 
partition function 
upper state energy 

Thus, combining Equations 4. 19 and 4.20 


Where 


Q 

E 


n 

u 




Eq. 4.23 


If several lines contribute to the absorption at a specific wavenumber, the 

total absorption coefficient is obtained by summation. The computational 

algorithm then uses S. as a basis, n A from Equation 4.21 gives 
6 -tu u u£ 

(n B - n B ) from Equation 4.9. This, in conjunction with cp . leads 

' l -tu u nl V 

to a via Equation 4.7. Equation 4.5 can then be integrated over altitude. 



4.2 Synthetic Spectrum - Fourier Transform Program 

4.2.1 Introduction : - Early in the development of this program it be- 
came necessary to devise a means for simulating the types of measurements 
to be performed. Because of the complexity of the atmospheric spectra, a 
digital computer simulation appeared most suitable. 

A computer program was required which would compute the 
transmission spectrum of the atmosphere over some predetermined infrared 
region for multiple gas species with overlapping spectra, taking into account 
that for certain strong transitions the atmosphere is not optically thin. The 
program must be able to include the effects of a bandpass filter to simulate 
the bandpass of the measurement instrument. Since the instrument is an 
interferometer, rather than a spectrometer, the Fourier transform of the 
spectrum must be obtained to simulate the instrument output. 

A computer program was written for use on the Langley Re- 
search Center - CDC computer system which computes theoretical absorp- 
tion-emission spectra of the earth's atmosphere as would be observed by an 
exo-atmospheric detector, and the Fourier transform of the spectrum. The 
program permits computation of a spectrum for three possible geometric 
configurations. Case 2 is for use when the instrument looks directly at the 
s\an through the earth's atmosphere. It calculates transmission, including 
absorption and emission of solar radiation through the atmospheric limb. 
Cases 1 and 3 are for use when the instrument looks downward. Case 3 cal- 
culates absorption of solar radiation reflected off the surface of the earth; 
while case 1 calculates transmission, including absorption and emission, of 
blackbody radiation emitted by the surface of the earth. Geometry of the 
cases considered are shown in Figure 4. 2.1. Geometry case 1 is used when 
blackbody radiation of the earth's surface dominates solar radiation. This 
occurs at wavelengths greater than about 4 y,. Since temperatures in the at- 
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Case 1. Blackbody Emission from Earth's Surface. 

Includes Absorption arid Emission in Intervening Atmosphere. 



Case 2. Limb Transmission of Solar Radiation. 

Includes Absorption and Emission in Intervening Atmosphere. 



Case 3. Solar Radiation Reflected from Surface of Earth, 
Includes Only Atmospheric Absorption. 


Figure 4.2.1 Geometries of Atmospheric Transmission Calculations 



mosphere are comparable to the surface temperature, the effects of emis- 
sion within the atmosphere itself must be included in the computation of the 
spectrum in addition to absorption. Conversely, geometry case 3 is used 
when solar radiation dominates blackbody radiation from the earth's surface, 
which occurs below about 4q,. In this case, the source temperature is far 
greater than the atmospheric temperature and hence only absorption in the 
atmosphere need be considered. 

Spectra are computed for atmospheres comprised of up to three 
gaseous components, having a total of 150 discrete absorption lines. Inputs 
to the program consist of a list of molecular absorption lines to be included 
in the synthetic spectrum, some thermodynamic properties of the constitu- 
ent gases, the physical properties of the model atmosphere to be employed, 
and the limits of the spectrum to be synthesized. Since the object is to sim- 
ulate the spectrum arriving at an exo-atmospheric instrument, it is also 
possible to superimpose a bandpass filter on the spectrum. Outputs from 
the program are tabulated spectra, plotted spectra and Fourier transforms, 
and the Fourier transform punched on cards for additional processing. 

4. 2. 2 Theoretical Problem Formulation : - A detailed formulation of 
the problem is given in Section 4.1. The problem consists basically of com- 
puting infrared atmospheric transmission spectra as determined by absorp- 
tion and emission from molecular rotational- vibrational transitions of the 
atmospheric constituents. 

As shown in Section 4. 1, the basic equation which is to be 
solved (Chandrasekhar, I960) is Equation 4.4. The first term of this ex- 
pression gives the absorption of incident radiation between the source and 
the top of the atmosphere. The second term gives the additional radiation 
contributed by emission in the atmosphere, as diminished by absorption by 
the part of atmosphere lying between the emitting part and the top of the 
atmosphere. This expression gives the net intensity at the top of the atmos- 
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phere at a single frequency. The program described herein solves this ex- 
pression for enough frequency values to define a spectrum over some de- 
sired frequency range. 

4.2.3 Program Structure : - The program solves the radiative trans- 
fer equations set forth in Section 4. 1 for an input set of molecular transi- 
tions and an input model atmosphere. 

Theoretical spectra may be calculated which cover frequency 
ranges of several hundred wavenumbers. Since at each frequency point the 
entire set of radiative transfer equations must be solved and integrated 
through the atmosphere, the problem program is clearly CPU bound. In 
order to keep processor time within reasonable limits, use is made of data 
tables, calculated initially and then combined in proper sequence to produce 
absorption and emission coefficients. Although this approach conserves 


execution time, it does so at the expense of core storage. To conserve 
storage, an overlay structure is used, resulting in the following storage 

quirements (numbers given are thousands of words. 

octal) : 

Main segment, including I/O buffers 

10 

Overlayed segments 

4 

Common block 

44 

System, including plot routines 

_1_5 

Total 

75 


This is equivalent to about 32,000 decimal locations. To fur- 
ther conserve storage, extensive use is made of the EQUIVALENCE state- 
ment to reuse storage locations of data tables. Some arrays are used for 
storage of four different tables at different points in the program. 

The program is written entirely in FORTRAN IV. An overall 


49 



view of general program flow is given below, followed by a detailed descrip- 
tion of individual components. A listing is given in Appendix A, 

4. 2. 3.1 Program Flow : - As noted above, the overall con- 
struction of the program involves the generation of data tables which, are 
combined to produce the desired absorption and. emission coefficients,,- which 
are then used to compute the overall transmission. The following is a gen- 
eral outline of the sequence of operations in the program. These steps fol- 
low to a great extent the development of the theoretical model given in Sec- 
tion 4. 1. „ , 

Step 1 . - Calculate Voigt Profiles - tables of constants are 
'calculated which are used in determining the Voigt function. Tables 
are calculated for varying values of the Voigt function parameters: 

X, the displacement between frequency of interest and the line cen- 
ter frequency; and Y, the ratio of Lorentz to Doppler broadening. 

' Step 2 . - Calculate Atmospheric Tables - model atmos- 

phere parameters are read in, including atmospheric tempera- 
ture, pressure and species concentrations as a function of altitude. 
Data are read in defining the desired geometry case; and atmo- 
spheric properties are computed for points along the radiation 
path. Tables of constants are calculated which will be used in 
computing the Lorentz and Doppler widths as functions of altitude 
and species type. 

Step 3 . - Calculate Line Tables - line information for all 
molecular lines to be included in the spectrum 'is read in and 
sorted by frequency. For each line a frequency band is defined 
which is the frequency region in which each particular line must 
be considered. 

Step 4.- Calculate Absorption and Emission Coefficient 
Factors - for each line at each altitude an absorption and (if re- 
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quired for the desired geometry case) an emission coefficient is 
calculated. The coefficients are Calculated for the line center 

i 

frequency. 

Step 5 . - Calculate Spectrum - for each frequency within 
the desired spectral region the frequency bands defined in Step 3 
are searched to determine which lines must be considered at that 
- particular frequency. For each line at each altitude a Voigt pro- 
file is calculated from the tables of Step 1, and multiplied by the 
absorption and emission coefficients of Step 4. Coefficients for 
overlapping lines are added, and the integration of equation (4.4) 

. is performed. After completion of one frequency point the fre- 
quency is incremented by a variable amount based on the proximity 
of line centers. If. a line center is close, then the transmission 
will be varying rapidly and small frequency increments are taken. 
Tabulations of the spectrum are printed. 

Step 6 . - Apply Filter and Take Transform - the spectrum 
calculated in Step 5 is multiplied by an input bandpass filter func- 
tion. The resultant spectrum is then Fourier transformed to give 
the interferogram. 

Step 7 . - Plot - the spectrum and interferogram are plotted 
on a Calcomp drum plotter. 

4. 2.3.2 Program Routines : - 

Main Segment . - Resident in core at all times is a short main 
program. This program, designated segment zero in the CDC segmentation 
scheme, performs no calculations, but simply calls the subroutines that per- 
form the calculations in the proper sequence. The main segment also re- 
serves space in blank common for the extensive arrays which will be used. 

A flow chart for the main segment is given in Figure 4.2.2. 
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Subroutine VOIGT . - The Voigt profile is a molecular line shape 
function which results from a superposition of independent Lorentz and Dop- 
pler line broadening (Armstrong, 1967). Lorentz broadening in the atmos- 
phere is a function of atmospheric pressure and is dominant at sea level, 
while Doppler broadening is temperature dependent and is dominant at high al- 
titudes. At altitudes between these two extremes there is a continuous degree 
of superposition of the two effects. 

The Voigt profile at a frequency, v , for a line centered at v is 
a function of two variables, X and Y, given by Equations 4. 16 and 4. 17. For 
the lines which will be included in the computation of atmospheric spectra, 
the range of values for X and Y which can reasonably be expected to be en- 
countered are 

O £ X £ io 4 
O £ Y £ 24 . 

As noted above, the general approach in this program is to com- 
pute data tables which will be combined to produce the desired spectra. Since 
it is known that Voigt functions will be required later with X and Y values fall- 
ing in the ranges defined above, tables of these functions are derived as indi- 
cated in Figure 4.2.3. Tables of 27 values of X and 25 values of Y are de- 
fined in a DATA statement. The Voigt function, $ (X, Y) is computed for all 
675 combinations of X and Y using a subroutine developed by Armstrong 
(1967)* It is also jcnown that interpolation of Voigt functions will be required. 
For this reason, the actual values of Voigt function are not stored. Instead, 
for a constant value of Y and three consecutive values of X, a quadratic func- 
tion is fitted through the three corresponding values of <J>, the Voigt function. 
The resulting quadratic equation is 

§(X.,Y.) = A..(X.) 2 + B. . (X.) + C.. Eq. 4.24 

11 ij 1 ij 1 ij 

The tables stored for later use are the 25 x 25 tables of A, B, and C. 
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Figure 4.2.3 Program SPECTRA First Segment 
(Subroutine VOIGT) 





Subroutine HEIGHT . - The primary function of subroutine 
HEIGHT is to set up tables of atmospheric properties along the radiation path. 
The routine reads data which defines the molecular species to be included in 
the computation of the spectrum. The program will accept from one to three 
species. For each species a data card is required which includes the symbol 
of the species, e. g. , CO, H^O, N^O, etc., the molecular weight of the spe- 
cies, the number of molecular lines which will be furnished for each species, 
and the partition function of the species. The partition function will be re^ 
quired in the program as a function of temperature. To facilitate this, the 
partition function, Q, is calculated for several temperatures from data in the 
JANAF Tables. It is found that over the temperature range of interest to this 
program, Q can be represented as a quadratic, i. e. 

Q (T) = AT 2 + BT ■+. C Eq. 4.25, 

The coefficients A, B, and C are evaluated manually and used as inputs to the 
program in place of the actual partition function. 

The next set of input data defines the geometric case to be con- 
sidered, as shown in Figure 4.2.4. If limb transmission (case 2) is indi- 
cated, then the grazing altitude, h (KM) , is defined. In case 1 the zenith an- 
gle of the observer a (degrees) is specified, while in case 3 the zenith angles 
of the source, a., and observer, /3 , are specified. The routine assigns a 
case number 4 to the special case of geometry case 3 where G! = /J. This re- 
sults in a shortening of computation time by a factor of two. 

Data is then read which defines the model atmosphere to be em- 
ployed. These data consist of a table of altitudes (KM) and tables which 
specify the following parameters at corresponding altitudes: temperature 
(degree K) , pressure (millibar), and concentrations of species specified 

earlier (particles /cc) . The set of altitudes for which these properties are 

• ' . / 

specified is arbitrary, but should cover the range 0 - 100 KM and must be 
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Figure 4.2.4 Program SPECTRA Second Segment 
(Subroutines HEIGHT, INTERP) 
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monotonic. These data are then combined with the geometric case data to 
compute tables of atmospheric properties along the actual line of observation. 
This procedure is similar for cases 1 and 3, but different for case 2, and 
will therefore be described separately. 

For the limb transmission case, case 2, the top of the atmos- 
phere is defined to be 100 KM. Using this condition, and the grazing altitude, 
h, the total path length through the atmosphere is calculated. This path is 
then divided into 48 equal segments, and the altitude at the end of each seg- 
ment is calculated. An interpolation is performed on the input atmospheric 
model tables for this altitude to obtain the atmospheric properties at this al- 
titude. A linear interpolation is performed for the temperature parameter, 
while a logarithmic interpolation is performed for pressure and species con- 
centrations. Counting the endpoints, this results in tables of atmospheric 
properties for 49 points uniformly distributed along the radiation path. 

For geometry case 1, a similar approach is taken in that tables 
of parameters are computed for points along the radiation path. In these 
cases the top of the atmosphere is' defined to be 96 KM. The atmosphere is 
divided into 48 altitude steps as follows; starting at 0 KM the first 16 steps 
are incremented by 1, KM each, the next 16 are incremented by 2 KM each 
and the third 16 by 3 KM each. The total radiation path length is computed, 
using the zenith angle, and the path lengths corresponding to the predefined 
altitude steps are calculated. Using the same interpolations as for case 2, 
the atmospheric properties at these altitude steps are computed. Again 
counting the end points, a table of atmospheric properties at 49 points along 
the radiation path is generated. For geometry case 3, this calculation is 
performed twice, for the entering and reflected radiation, using the two dif- 
ferent zenith angles. If the zenith angles are equal, then only one table is 
computed, but in effect it will later be used twice. 

The Doppler halfwidth is given by Equation 4. 18. Since will 
be repeatedly calculated later, considerable time can be saved if as much as 
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possible of this expression can be precalculated. • Therefore, a table of Dop- 
pler width factors is computed for .each .path position, which determines T, 
and each species, which determines M. Hence, •; . 


DOPP (I, Jj 


1_ 

C 


2RT.ln2 

J 

M. 

, . l- . 


1/2 


Eq. 4. 26 


It is now only necessary to multiply the approximate DOPP value by v , the 
line frequency to obtain for an transition at any point in the atmosphere. 

Similarly, the Lorentz halfwidth, a , for any line is assumed 

_ 1 1 ... 

to be 0.06 cm at STP. Tables of a. are therefore calculated for all points, 

L 

j, along the. radiation. path by 


whe re 


a£(j) = 0.06 


P. 



s 



P 

s 

T 

s 


standard pressure = 1000 mbar 

standard temperature = 273 °K 


Eq. 4. 27 


Subroutine LINER . - Subroutine LINER reads data concerning 
the individual lines to be included in the computed spectrum, sorts the lines 
by wavenumber, and computes the frequency band over which each line con- 
tributes to the spectrum. 

In Subroutine HEIGHT, up to three chemical species were de- 
fined as contributing to the spectrum. Subroutine LINER reads the set of 
transition lines to be considered for each species. Up to 150 lines are per- 
mitted, distributed in any way among the three species. For each line it is 

-1 -2 

necessary to specify the center frequency (cm ), the line strength (cm 
atm 1) and the lower state energy (cm *). 
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After reading this line data, LINER sorts the lines and corre- 
sponding line data in order of ascending center frequency. A search is per- 
formed to determine the minimum line strength. For each line a frequency 
band is determined over which the line must be considered in computing the 
spectrum. The line inclusion band is determined by fitting a Gaussian line 
shape over each line with a nominal 0.0 6 cm ^ halfwidth. The line inclusion 
band is the frequency band over which the line has a strength which exceeds a 
previously defined minimum. This minimum is presently taken as 1% of the 
strength of the weakest line. If the line inclusion band is computed to be 
greater than ± 5 cm it is truncated to ± 5 cm 

Subroutine COEFFS . - Subroutine COEFFS computes, for each 
line at its center frequency and for each altitude in the model atmosphere, 
the absorption coefficients and, if required by the geometry case being con- 
sidered, the emission coefficients. Using the definitions of Einstein coeffi- 
cients, the absorption and emission coefficients for line SL of species i at al- 
titude h can be shown to be: 
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where 




concentration of species i at altitude h 


/ particle ^ 

cm 


s chematic /numbe r 
lower state energy of line t (cm *) 
second Planck constant 
standard temperature = 300°K 


temperature at altitude h (°K) 

center frequency of line £ (cm *) 

-2 -] 

strength of line t (cm atm ) 

partition coefficient of species i at the temperature 
at altitude h 


These factors do not represent the complete absorption and 
emission coefficients since they are lacking the Voight line profile function. 
This last term will be included in the next subroutine. 

Subroutine SPECT . - Subroutine SPECT performs the actual 
spectrum evaluation. For each frequency point in the spectrum, the precal- 
culated data tables are accessed, and equation i is evaluated, giving the ob- 
served intensity. 

In order to keep the computer running time within reasonable 
limits, it is necessary to restrict the total number of frequency points at 
which the spectrum is computed. On the other hand, it is necessary to in- 
clude a number of points sufficient to maintain adequate spectral resolution, 
particularly in the rapidly varying regions of line centers. To satisfy these 
requirements, a frequency increment system was devised whereby large in- 
crements, up to 0.512 cm 1 are used in those parts of the spectrum distant 
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from line centers, while small increments, as small as 0.001 cm are used 
in the vincinity of line centers. The frequency increment scheme operates as 
follows. 

The first frequency point is defined by input data. After the in- 
tensity is computed, for this frequency value, as will be described below, a 
search is made of the table of line center frequencies to determine the closest 
line center, at either higher or lower frequency. The absolute value of the 
frequency difference between the nearest line center and the current frequency 
point is used in an algorithm to compute the next frequency increment. If this 
difference is zero, an increment of 0.001 cm * is used. At every second 
point from the line center the increment is doubled until a maxim um of 0.512 
cm * is reached at 1.0 cm * from the line center. The computed increment 
is added to the current frequency value, the spectral intensity is computed for 
the new frequency, and the process is repeated until a final frequency, de- 
fined by input data, is reached. 

Each group of 100 spectral points is printed and also written on 

J 

a magnetic scratch tape. This permits the computation of spectra with large 
numbers of points without extensive storage requirements. 

At each frequency value, the spectral intensity is computed as 
follows. The table of line inclusion bands, determined by subroutine LINER, 
is searched to determine which lines must be included in the computation of 
the spectral intensity at this frequency. A table is made of lines to be includ- 
ed. For a given altitude step in the model atmosphere, the Lorentz and Dop- 
pier line widths are obtained for each line in the line inclusion table. The 
Lorentz width was computed in subroutine HEIGHT. The Doppler width is 
computed by multiplying the appropriate entry in the DOPP table, also com- 
puted in subroutine HEIGHT, by the corresponding line center frequency. 

The Doppler and Lorentz line widths are used to compute the X and Y param- 
eters of the Voight function, as given in Equations 4. 16 and 4. 17. 
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These X and Y values are used in a table lookup of the Yoigt 
function tables generated in Subroutine VOIGT, and, using a three point La- . 
grange interpolation, the Voigt function is evaluated for a given line at a par- 
ticular frequency and position in the model atmosphere. The Voigt function is 
multiplied by the appropriate absorption, and if required, by the geometry 
case in use, the emission coefficient factors generated in Subroutine COEFFS, 
giving the final absorption and emission coefficients for a given line at a spec-, 
ified frequency and altitude. If the search of line inclusion bands indicated 
that more than one line must be included at this frequency, i. e. , the lines 
overlap, this process of determining absorption and emission coefficients is 
repeated for each overlapping line. The absorption and emission coefficients 
are summed for all contributing lines to give a single net absorption or emis- 
sion coefficient. 

This process of computing absorption and, if required, emission 
coefficients is repeated for all altitude steps in the input model atmosphere, 
as defined in Subroutine HEIGHT. A table of absorption and emission coeffi- 
cients is generated for a single frequency covering all altitude steps and in- 
cluding the effects of overlapping lines. The radiative transfer equation 
(equation 4.10) is evaluated using the tabulated coefficients, employing Simp- 
son's rule integrations. The incident radiation for geometry case 1 is black- 
body radiation from the ground. The ground temperature is taken as 288.2 
°K, with an albedo of 1.0; however, both parameters can be varied by input 
data. For geometry cases 2, 3, 4, the incident radiation is solar radiation 
based on a solar temperature of 5800 °K, with an emissivity of 1.0. For ge- 
ometry case 3, where incident solar radiation is reflected from the ground at 
an angle different from the incident angle, the intensity computed in the man- 
ner described above gives the attenuated solar radiation at ground level. The 
entire process is repeated for the reflected beam to give the radiation at an 
exoatmospheric observing platform. For geometry case 4, where solar radi- 
ation is reflected from the ground at the same angle as the incident angle. 
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the computed absorption is simply doubled to give the total absorption from 
the double pass. In both geometry cases 3 and 4, a ground reflectivity may 
be specified by input data. The intensity at ground level is further reduced by 
a factor of ff to account for lambertion diffusion by the ground. 

By this method the spectral intensity at a single frequency is 
computed. The subroutine then computes the next frequency point by the line 
center proximity method described above, and the entire process is repeated 
for the new frequency point. This procedure is continued until a stopping fre- 
quency, specified by input data, is reached. 

Segment Six . - The sixth segment of this program contains three 
subroutines. Subroutine TRANF prepares the spectrum computed above for 
the Fourier transform, applies a bandpass filter function as defined in Sub- 
routine FLTRS, and calls Subroutine FRXFM to perform the Fourier trans- 
form. Subroutine TRANF then prepares both the spectrum and interferogram 
for plotting. 

In the description of Subroutine SPECT it was noted that the 

spacing between points on the frequency scale of the computed spectrum is 

variable father than uniform. In order to perform the Fourier transform, 

however, it is necessary that the frequency spacing between points be uniform. 

The first task of Subroutine TRANF, therefore is to perform conversion of 

the spectrum to a uniform scale. The interval between the spectral limits is 

divided into a number of uniform divisions. The number of divisions must be 

12 

a power of two, and is taken as 2 unless otherwise specified by input data. 
The spectrum is read from the intermediate scratch tape in blocks of 100 
points, as written by Subroutine SPECT. The conversion to the uniform scale 
between the limits of each uniform interval, i. e. : 
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where 
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Following conversion of the spectrum to a uniform frequency 
scale, a filter function is applied to simulate the bandpass of a measurement 
instrument. Several types of filter functions are available, the choice being 
specified by input data. A Lorentz or Gaussian shaped filter function may be 
specified, in which case the input data consists of the center frequency, the 
peak transmittance and the half height - half width. A problem has been noted, 
however, in using these two functions in conjunction with the Fourier trans- 
form. Because the Lorentz and Gaussian functions asymptotically approach 
zero, but do not become zero in the defined spectral range, there is introduced 
into the Fourier transform high frequency components. These components can 
be eliminated by forcing the filter function to zero while still within the spec- 
tral range. This is accomplished by using a third filter function, referred to 
as a power function, having the form: 


Spectral intensity of point in spectrum 
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where the center frequency, p , the width, w, and the exponent, n, are all 
specified by input data. This function has the property of becoming zero with- 
in the frequency range of the computed spectrum and thereby preventing high 
frequency components in the Fourier transform. 
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In the event that a filter bandpass is required which does not fit 
any of these three analytic forms, it is also possible to provide input filter 
data in the form of the calibration curve, i. e. , a table of frequencies and cor- 
responding transmittances may be specified. Up to fifty calibration points 
may be specified. 

It is also possible to superimpose any number of filters, using 
any mixture of the above filter types. 

Following the application of the filters to the computed spectrum, 
the spectrum is written on magnetic tape for subsequent plotting. The Fourier 
transform of the spectrum is then computed. The Fourier transform subrou- 
tine, FRXFM, is a standard library routine based on the Cooley-Tukey algo- 
rithm, and returns the real and imaginary components of the Fourier trans- 
form. These are converted to amplitude and phase components of an interfer- 
ogram by the expressions 



where 

R, I = real, imaginary components of Fourier transform 

S = plotting scale factor. 

The computed interferogram looks much like an amplitude modu- 
lated radio signal in that it consists of a carrier wave modulated by the rela- 
tively lower frequency envelope. In the interferogram the carrier contains in- 
formation concerning the spectral region for which the interferogram is com- 
puted. Since this information is already known, it is advantageous to eliminate 
the carrier and retain only the envelope, since this relaxes the required ac- 
curacy of the sampling point. 
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The elimination of the carrier is performed by "synchronously 
detecting" it with an input reference frequency which lies approximately in the 

middle of the spectrum being computed. 

The delay scale for the interferogram is based on the frequency 
range of the initial spectrum. The delay for the first interferogram point is 
zero, and the delay increment for each subsequent point is 


where 


Ad 


1 




V , = initial, final frequencies of the computed spectrum. 

The amplitude, phase and delay tables of the interferogram are 
written on magnetic tape for subsequent plotting. In addition, the real and 
imaginary Fourier transform components are punched on output data cards for 
subsequent use in deriving interferogram correlation functions. 


Subroutine SPPLOT . - Subroutine SPPLOT reads the spectrum 

n 

and interferogram from the magnetic tape written by previous subroutines, and 
generates plots of them. Plots are made on the twelve inch CALCOMP drum 
plotter at Langey Research Center, using the standard system library of plot- 
ting routines. 

Following generation of the plotted output, program control is re- 
turned to the main segment where the entire spectrum computation process is 
repeated if desired. ' 3 1 
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5. CORRELATION INTERFEROMETRY 

5. I Principles of Interferometry 

The use of spectral techniques for the remote measurement of 
concentrations of trace atmospheric species is dependent on separating out 
the effects of the species being measured from those of all other species 
present in the optical path. In techniques where a part of the spectrum is 
measured, the separation must be obtained by spectral resolution. Thus, 
some separable part of the spectrum must show an appreciable effect of the 
species being measured and any significant effects of other species must be 
such that they can be eliminated. If a different technique is used effects 
must still be separable but the separation criteria is no longer spectral reso- 
lution but a type of resolution peculiar to that technique. Such a technique is 
interferometry. In this technique the separation is a'ccorriplished by a reso- 
lution of path differences. 

.The instrument being developed for this program is a correla- 
tion interferometer. The following discussion will describe the basic theory 
of its operation. More detail on this and on the specific instrument itself 
will be given elsewhere (Bortner, et. al. , 1973). , 

An interferometer is a fairly simple device (Figure 5.1. 1). The 
essential elements are a beam splitter and two mirrors, plus a detector to 
measure the radiation output. Light from the source is incident on the beam 
splitter, B. At the beam splitter it is divided into two paths; one portion of 
the light goes to one mirror, the other portion of the light goes to the 

other mirror, The two portions recombine at the beam splitter and the 

intensity of the light once they recombine is registered by the detector, D. 

The intensity of the radiation received by the detector will depend on the dif- 
ference between the lengths of the paths traveled by the beams in the two 
arms. The length of the path F-B-M^-B-D can be different from that of the 
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path.F-B-M^-B-D. If the two optical paths are exactly the same, the path 
difference (delay' 1 ') is zero, and there is a peak in intensity. If monochro- 
matic radiation enters the instrument, and if the path difference is increased 
by one -half the wavelength, the intensity reaching the detector goes dov/n es • 
sentially to zero. Then as it increases again towards one wavelength path 
difference, another peak occurs. This sinusoidal oscillation about a mean 
level repeats at intervals of one wavelength if the light is monochromatic. 

The instrument actually does a Fourier transformation of the spectrum of the 
radiation entering. The way in which the delay is scanned, the way in which 
the path difference is changed is in most instruments a matter of shifting one 
of the mirrors. This certainly does the job, but it causes some unnecessary 
alignment problems. One of the features of the COPE instrument is that the 
problem of having to scan and maintain the position of the mirror accurately 
is avoided by not scanning the end mirror, but, instead, scanning a plate of 
refractive material in one arm of the interferometer. Many interferometers 
have such a refractive plate, but generally it is left in a constant position. If 
this plate is rotated, the path length in that arm of the interferometer will 
vary. This will accomplish the same effect as moving one of the mirrors 
back and forth without the alignment problems. This is a specific advantage 
of the technique of this instrument. Other advantages, some of which are 
advantages of interferometry in general, include a large light throughput, 
multiplexing of spectral elements, a compact yet flexible instrument, and a 
handy output. 

5. 2 Relationship Between Spectra and Interferograms 

In Figure 5.2. 1 a few particular cases of the relationship between 
spectra and interferograms are shown. As shown on the previous figure, 
monochromatic spectrum produces an interferogram which is essentially a 

* The term "delay" refers to the temporal variation of the interferometer sig- 
nal. This is related to the variation of the path difference by the scan velocity. 
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sinusoidal variation as the delay is scanned. If instead of a single monochro- 
matic line there is a pair of monochromatic lines, then the two sinusoidals 
due to lines will beat together. Again a characteristic maximum is obtained 
at zero delay. As the oscillations beat together, they get out of phase and go 
through a minimum, then maximize again at a point characteristic of the line 
separation. If there are a number of regularly spaced lines, the beat pat- 
terns will all reinforce at zero delay again, but they very. rapidly get out of 
step with each other and decay to very small values. Then there appears a 
point where they all reinforce again at a delay which is inversely proportional 
to the spacing between the lines. This is almost exactly the effect seen in the 
case of carbon monoxide as shown in Figure 5.Z.2, which gives the actual 
spectrum' 1 ' and interferogram. The spectrum displays fairly regularly spaced 
lines. There is certainly a noticeable change in the spacing from one end to 
the other, but it is certainly not a random spacing. As seen in the previous 
figure, there is reinforcement for zero delay which dies out rather rapidly, 
and then, at a delay which is characteristic of the spacing between the lines 
(about 3 cm ^), the interferogram amplitude peaks up again. In general the 
relationship between the spectrum of the radiation and the interferogram of 

t . " . i 

the radiation is given by a Fourier transformation, the cosine Fourier trans- 
formation. That is, the interferogram, apart from the constant term, is the 
cosine Fourier transform of the spectrum (Figure 5.2.3). The interferogram 
signal, as ' a function of path difference is given by: 

I (X) = 2 J 0 n. (s.(a)) cos (2 ttcxX) da 

where 

CT = frequency (cm *) , 

2 -1 -1 ” 1 -1 
(a) = spectral input at a (ergs-cm -ster (cm ) -sec ) 

X = path difference (cm) 

*This spectrum is that given by Plyler (1952) with the mercury emission 
lines omitted. 
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SPECTRUM 



Figure 5. Z. 3 Heterodyning of Interferogram 



Several relationships exist between the interferogram and the 
spectrum from which it is produced. The interferogram has a carrier fre- 
quency directly related to the mean frequency of the spectrum producing it. 

It will be amplitude and phase modulated, the amplitude and phase modulation 
being described by the band envelope shown. From an interferogram up to a 
given delay, the spectrum could be reconstructed with a resolution given in- 
versely by the delay out to which the interferogram is obtained. So, from an 
interferogram from zero to 1 centimeter, inverse Fourier transformation - 
reproduces a spectrum to resolution of one wave number (cm ?). Conversely, 
if the information characteristic of one wave number resolution is not needed, 
if the information need only be characterized by a 3 cm ^ resolution, then it 
is only necessary to scan the interferogram out to 1/3 cm (3 millimeters) to 
obtain information equivalent to a resolution of three wave numbers. It is 
very important to do this. Significant improvements in signal to noise ratio 
are obtained by limiting ourselves to scanning only the portion: of the interfer- 
ogram which gives the best signal to noise, ratio. If essentially all the infor- 
mation on any given species (all the effect of that species on the interfero- 
gram) occurs over a small part of the interferogram, only that range of path 
difference need be scanned. The operation of the correlation interferometer 
involves the treatment of the interferogram data directly to obtain data on 
species densities rather than the use of the spectrum obtained by the Fourier 
transform of the interferogram. With such an operation the concept of spec- 
tral resolution loses meaning. Although the spectrum which can be obtained 
from an interferogram has a resolution which is the reciprocal of' the path 
difference scanned, the ability to obtain data on one species in the presence 
of others by direct examination of the interferogram, is not dependent on the 
length of path-difference scan. For example, if a scan of 0 to Z cm were 
used to produce an interferogram but a given species of interest has an effect 
only between 1.0 and 1.2 cm, as much information on that species can be ob- 
tained by scanning from 1.0 to 1.2 cm as- from 0.0 to 2.0 cm even though the 
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longer . scan could produce a spectrum of 0.5 cm resolution whereas the 
shorter scan could not. Rather than spectral. resolution, the separation of ef- 
fects of various species is accomplished by a spatial resolution of path differ- 
ence which is dependent on such factors as smoothness of path-difference scan. 

Another relationship between the spectrum and the interferogram 
is that the spectrum of the incoming radiation can be severely band-limited by 
optical filters. If it is band-limited, and if an interferogram that would re- 
produce the spectrum only to a limited resolution is taken, then the spectrum 
could be characterized by a fairly small number of points. On the other hand, 
to characterize the interferogram in this form with its very rapid oscillations 
(the oscillations are essentially proportional to the mean frequency and not to 
the width) with any degree of accuracy, relatively large number of points on 
the order of at least one per cycle of oscillation would be needed. This would 
be many more points than are necessary to characterize the spectrum to a 
corresponding resolution. There is a lot of redundancy in the information be- 
cause of the band limitation. The method used in the correlation interferom- 
eter to eliminate the redundancy is simply to take the amplitude and phase 
modulated sinusoid and heterodyne it down with a local oscillator, where the 
local oscillator is the interferogram of radiation somewhere around the mean 
frequency of the spectrum. By heterodyning the interferogram down with a 
cosine or a sinusoidal variation, the interferogram is reduced to its essential 
variations. This is illustrated in the bottom part of Figure 5.2.3. All the in- 
formation necessary to characterize the interferogram can be retrieved by 
sampling a relatively small number of points, a few points for each of the 
much longer cycles shown in the bottom figure. There is a slight difference 
whether the interferogram is beat with a sinusoid at a given wavelength or a 
cosine; slightly different beat-down interferograms are obtained. Actually 
both of them are used but, through most of the remaining discussion this fact 
will be ignored. In generating the local oscillator for carrying out the hetero- 
dyning, radiation which passes through the same interferometer as does the 
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signal radiation is used. This insures that the local oscillator always has the 
correct phase relationship with the signal interferogram. It also relaxes 
some of the accuracy needed in the knowledge of the plate drive (the scan 
drive). If this were not done, the scan drive would have to be known accu- 
rately to within a very small fraction of the wavelength. By having the local 
oscillator going through the same interferometer, this problem is alleviated. 
Probably the main objection to interferometers is that most people think in 
terms of using a spectrum and in order to obtain a spectrum from the inter- 
ferometer, the interferogram must be transformed. This problem is avoided 
completely by not looking at the spectrum at all. It is not really necessary to 
use the spectrum. The measurement can be made quite adequately on the in- 
terferogram itself. This is a major advantage of the correlation interferom- 
eter. 


5. 3 The Measurement in the Presence of Interferents 

Figure 5. 3. 1 outlines very briefly the basic principles by which 
the measurement on the interferogram is made. There are essentially two 
problems involved in a measurement of this type. We are trying to measure 
the CO burden on the basis of the radiation received from a satellite. The 
problems in making the measurement are: (1) can the measurement be made 
with sufficient accuracy in view of the noise limitations, and (2) can the mea- 
surement be made under conditions where the radiation received is affected 
not only by the gas that we are trying to measure, in this case carbon mon- 
oxide, but is also severely affected by other gases. In fact, the radiation is 
much more affected by such things as water vapor and methane. Ignoring for 
the moment this problem of the interference species, consider how a mea- 
surement would be made. Consider an interferogram due to a target signal, 
say CO, such as shown. There are different ways by which one could mea- 
sure the CO burden which produced this signal level. One might sit at a con- 
stant delay and measure the signal level at that point. If the CO burden were 
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doubled, the effect of CG on the intensity we would measure would be doubled 
(for the' moment assuming variations are linear with the gas burden); the in- 
tensity would essentially be- a measure of the carbon monoxide burden. If 
there were a constant noise level associated with the measurement of any 
point then the best point to make a measurement would be at the point where 
the; signal' level is maximum. . If a finite range of delay were scanned, then 
the measurements made at all points would' be combined to get some sort of 
average measurement of the carbon iriOnoxide burden.' In general, the' opti- 
mum ^measurement that can be made in ; such circumstances is given by com- 
bining all the measurements, all the intensities- of the various points ,- in a 
manner related directly. to the intensity of the signal shape. That is a weight- 
ing function'or correlation function is generated and this is multiplied togeth- 
er with^the signal. The measurement rnaide is the integral of the correlation 
function times the signal, the integral- over the delay range which is scanned. 
It is integrated and the measurement obtained is directly proportional to, in 
this case, the carbon monoxide 'burden.' The' signal-to'- noise ratio in such a 
measurement can be shown to be optimum when the correlation function looks 
exactly in shape like the target signal.. ..That .result holds when there are no 
interferents. Now. consider another gas species which is affecting the signal 
received, for example water vapor. If a point measurement of the intensity 
were used, the signal received would not be simply the signal due to carbon 
monoxide but the , sum of the signals due to the carbon monoxide and to the 
water vapor. Even if the. carbon monoxide level were to remain the same, 
drastically varying measurements due- to variations, in water vapor might be 
obtained so that a very poor, measurement of CO would be obtained by sitting 
at that point. A region where a signal is unaffected by all interferents may 
not be available or if it is, the signal-tomoise ratio achievable at. the point 
might not be adequate for the measurement. This problem is handled as fol- 
lows, The correlation function used is. not matched exactly to the target gas 
signal shape. The correlation function (W) is adjusted so that when it is 
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cross-correlated with the interferents and the result is integrated over the 
range scanned, all the positive correlation regions are balanced exactly by 
the negative correlation regions, so that the total area under the curve, the 
total result of the measurement, comes out to zero as illustrated, still 
maintaining the correlation function as close as possible to that constraint to 
the signal interferogram (target interferogram) in order to. still get as large 
a positive correlation between the correlation function and the target gas as 
possible. The result of the measurement is still proportional. to the carbon 
monoxide burden. In principle, if there are a number of interferents of this 
sort, rejection of these interferents can still be achieved so long as we have 
at least as many points to describe our correlation function as we have gases 
which are affecting the radiation. The measurement of the target gas can be 
made in real time. The first gas, in this case CO, and any other gases 
which significantly affect the interferogram can be measured. 

5. 4 An Example of the Use of Interferograms 

As an example of the interferograms resulting from overlapping 
spectra (using three imaginary species), Figure 5.4.1 shows spectra (drawn 
with relative frequency scale in cm and the major features of a portion of 
the interferogram envelopes for species A, B, and C, and for all combina- 
tions of them. Thus, from the spectrum of A, the lines at relative positions 
of 1 and 3 cm * give a peak in the interferogram envelope at 0. 5 cm; those 

at 1 and 6 cm * give one at 0. 2 cm * ; those at 3 and 6 cm * give one at 0. 33 

cm Similarly, B gives peaks at 0.28 cm and 0.55 cm. However, in the 

interferograms of AB an additional effect occurs between . 25 and . 36 cm. 

This is due to the combination of the line at the relative frequency of 1 cm * 
in the A spectrum and of that at 4. 8 cm * in the B spectrum. Other effects 
are also present in the AB interferogram due to other combinations of lines. 
If it is desired to determine the concentration of species C in the presence of 
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A and B, the spectrum or the interferogram of the A-B-C mixture must be 

different from that of the A-B mixture. The spectra show an appreciable 

1 -1 -1 
difference only around v> , = 5 cm where A-B-C extends about 0.2 cm 

rel 

past A-B. Thus, a .2 cm'* spectral resolution would be required to make 
the measurement. In the interferograms , the region around .25 cm path dif- 
ference shows a variance.' 1 The region from . 24 to .26 could be used to make 

a measurement C. Actually a larger region would be used and the effects of 

£ 

A and B separated out. 

The spectra of A, B, and C are related to each other somewhat 
as those of H O, CH , and CO in a real case, except that the latter are more 
complex. In that case a spectral region much affected by CO is chosen and 
separated out by an optical filter. 

In the correlation interferometer, the path difference scan is 
limited principally to that portion of the interferogram which is most affected 
by the gas whose density is being measured (excluding that portion near zero 
path difference). In the case of CO that portion of the interferogram obtained 
with path differences between 2.70 and 3. 95 mm has been selected. Since the 
lines in the first overtone band of CO are separated by about 3 cm ^ on the 
average this shows up in that part of the interferogram centered at about 1/3 
cm * = . 33 cm = 3. 3 mm. Thus, the CO spectrum is as shown in Figure 
5.2.2 and the interferogram (using a spectral filter centered at 4278 cm *) is 
also shown in Figure 5. 2.2. In this the major effect of the CO is seen to be 
in the region of 0. 3 cm. Taking this portion of the interferogram between 
0. 27 and 0.395 cm gives essentially all of the information on CO that can be 
obtained from the interferogram. The determination of CO is not, however, 
the simple measurement of this peak in the interferogram, since other gases, 
particularly CH and H O also have some effect in this region. Thus, this 

Tt £ 

peak is divided up into a number of sections, 32 in the case of our instrument. 
The effect of each gas on each of these 32 sections is determined in effect, 
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by calibration.* 

Thus, in the simplest case, a, portion of the interferogram af- 
fected only by HO is chosen to determine the amount of H O** and the re- 
mainder of the interferogram is, in effect, corrected using these data and 
calibration information. CH^ is dealt with similarly. Thirty-two different 
effects could be dealt with this way. Since in a real case portions of the in- 
terferogram affected by a single gas are not always available, the set of 32 
simultaneous equations is solved to determine the amounts of individual gases 
causing the combined effects. 

The output of the interferometer will be voltages, each related 
to the concentration (column density) of one species. The calibrations can 
be thought in terms of these voltages. Thus, to be able to separate out the 
effects of individual species it is necessary to resolve parts of the interfero- 
gram which show the effects of these species. In the A-B-C example this res 
olution must be of the order of 0.2 mm which is readily attainable. The spec- 
tral resolution required for a similar separation of effects would be 0.2 cm 
which is not attainable with the spectrum which would be produced from the 
interferogram obtained with this instrument. Thus, by directly using the in- 
terferogram the effects of specific species and hence its density is measure- 
able whereas it would not be if the interferogram were converted to a spec- 
trum and that used for the measurement. This comparatively simple mea- 
surement can be made very much faster than can a measurement of the same 

^Actually, the effects of specific gasbs, other than that being measured 
(such as CO), need not be specifically known by calibration. The calibration 
can consist of 31 different atmospheric conditions with one amount of CO and 
a different amount of CO with one of these conditions. The measurement of 
CO is then made by determining the effect on each section of the interfero- 
gram of other gases of the atmosphere ad determined by a best fit of the 32 
points. 

;,;5,5 Since the H£0 spectrum and hence its interferogram is affected by temper- 
ature, two or more portions of the interferogram must be used to take this 
effect into account. 
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species by a high- resolution spectromatic technique, whether by a spectrom- 
eter technique or' by an interferometric technique which employs the Fourier 
transform of a interferogram obtained over a large range of path difference. 
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6. CORRELATION TREATMENT 

6. 1 Basic Principles 

6. 1. 1 Instrument Output: - The interferometer receives radiation, 

, — . — <n — j — !— 

the spectral characteristic^ of which are affected by a number of parameters. 
We wish to measure one of these parameters (CO burden), with as high ac- 
curacy as possible and as little interference as possible from variation in the 
other parameters. Within the instrument, optical filters limit the range of 
the radiation to a fairly narrow spectral region where CO absorption has the 
greatest relative effect on the radiation. 

The correlation interferometer operates by yarying the de)ay 
between the two beams of a Michelson interferometer over a defined range. 

As the delay is swept, a sinusoidal type signal is generated at the detector 
output. This signal is both amplitude and phase modulated. The band center 
and width are directly related to the center and width of the optical spectral 
filter, while the modulations are characteristic of the more detailed spectral 
information. 

6. 1.2 Preprocessing : - An important function of the interferometer 
electronics hardware is to reduce the high sampling rate implied by the rel- 
atively high center frequency of the interferogram signal, by synchronously 
detecting the signal with a reference signal, (similar to the local oscillator 
in a homodyne system). This reference is generated from the interferqgram 
of radiation having stable spectral characteristics. The reference radiation 
passes through the same interferometer system as the signal radiqtiop, so 
that to a very large extent the reference derived from its interferogram is 
fixed in phase relative to the interferogram of any given signal radiation, 
regardless of any irregularities in scan, or variations in the interferometer 
arm lengths. 
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When this synchronous detection is carried out, the result is a 
signal' containing only the information of the sidebands if the original signal, 
beat down to zero frequency mean. (Actually two such signals are produced, 
from references both in phase and in quadrature phase with the reference in- 
terferograrh. ) T ! he signal's information is still bind limited, by virtue of the 
original spectral band limitations. The signal is therefore integrated and 
digitized at an interval sufficient to retrieve this information, producing a 
preprocessed, digitized, iriterfe /ogram. * 

6. 1.3 Final Processing : - Final processing consists of applying a 
linear digital filter ("weighting function" or "weights") to the preprocessed 
interferogram. The filtering operation consists of taking the weights, one 
for each point of the interferogram, multiplying them together with the cor- 
responding interferogram points 1 , and summing the results. The weights are 
chosen, in a manner to be described below, so as to give a final result which 
is insensitive to variations in all parameters except the desired one. (The 
weights are also chosen to give a result which is optimized against a mixture 
of noises. ) The result may be directly converted using a defined zero point 
and scale factor, to units of target gas burden (atm- cm* ppm-m); etc. 

6. 2 Determination of Weights 

6.2. 1 Basic Philosophy and Theory : - The determination of the 
weights is based on the possibility of representing any observed interfero- 
gram as a linear combination of component interferograms. The number of 
such component interferograms required to represent a range of actual ob- 
served interferograms to some accuracy, will depend on a number of factors 
depending on the range of conditions under which the interferograms are ob- 
served. The number will depend mainly on the number of parameters which 
vary significantly over the range of conditions, and to a lesser extent on the 
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degree of non-linearity in the variations of the interferograms with these 
parameters. (Parameters will characteristically represent factors such as 
water vapour burden, methane burden, and temperature profile, as well as , 
the target CO burden). 

Consider then any interferogram I. within the defined range. It 
may be written as a linear combination of N constituent interferograms. 
Using the index K to define the point on the interferogram, and J to define 
the particular constituent interferogram we may write 

N 

I(k) = T qj Ij(k) 

■ ; . ’■ j = 1 : 

The q 's are the strengths of the constituent interferograms in the particular 
•J . 1 

interferogram observed. Let us assume one of these, q , is the target gas 

J X 

burden. . 

Now the final processing consists of the application of a weight- 
ing function, H(k), to the interferogram. The resulting measurement is 

M = £ ,H(k) I(k) 

k 

= £ q, E H(k) I (k) . 

J J J 

From this, it can be seen that if we can find weights H(k) that give 
• 1 Nk 

t H(k) I (k) = , 1 J = J1 . ‘ 

k = 1 ^ 

, . « , 0 J * J1 

then the measurement will be M = q , the target burden. If J varies from 

J i 

1 to N, then this is a set of N equations in as many unknowns H(k) as there 
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are interferogram points (NK). In order for these equations to be satisfiable, 
there must be at least as many interferogram points as there are constituents 
(Nk ^ N). In the case that Nk = N, there are N linear equations in N un- 
knowns, and the solution for the H(k) is straightforward. 

In the case that there are more interferogram points than there 
are constituent interferograms , there are many possible weighting functions 
which will satisfy the N equations. It is however possible to choose a unique 
set of weights by considering the noises in making the measurement, and 
minimizing their effect on the measurement. 

6.2.2 Noises in Measurement: - The noises operate in the following 
manner. The total integration time for measurement of a given point may be 
written D(k) (its duration, assumed unity up till now). The contribution to 
the final measurement for that point is actually AM(k). 

AM(k) = j^H(k) D(k)] I(k) 

The noise contribution for that point is significant in how it can vary this 
from its "correct" value. The noises may or may not decrease with inte- 
gration time: that is, they may be either random or synchronous with the 
scan. They may or may not depend on the expected level of the interfero- 
gram itself <I> , (an RMS value for the constituent interferograms at a 
given point): that is, they may be either additive or multiplicative. 

A discussion of various physical origins for these noises will 
be presented elsewhere (Section 6.2.4). However, it is sufficient at this 
point to be able to characterize them by four numbers, a^ .... a^, which 
represent the severities of the four possible combinations of noise type: 
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Random Additive 


(RA) 


- Random Multiplicative (RM) 

a^ - Synchronous Additive (SA) 

a^ - Synchronous Multiplicative (SM) 

The noises will combine in an RMS manner to give the mean square noise 
2 

contribution AN (k) to the measurement at a particular interferogram point, 
and we may write 

AN 2 (k) = [H(k) D(k)] 2 x [( aj + a., <I>^)/D(k)' 

+ (a 3 + a 4 

= [H(k) D(k)] 2 G(k) . 

G(k) is an effective mean square interferogram error for point k (although 
it may originate through error in either H, D, or I), 

6.2.3 Derivation of Weights : - 

The total measurement is thus 

M = E AM(k) = E |~HD ll 
k k L J 

and the total squared noise is 

N 2 = E AN 2 (k) = eThdI 2 G. 
k k*- 


It is possible to choose, given the various Ij(k) and the duration D(k), those 
values for which minimize the above noise term; subject to the con- 
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straint that M = q , the target quantity, regardless of what values the 

J I 

other quantities may have (to the extent that the assumed linear representa- 
tion is valid).' The solution, derived elsewhere (Grenda, et al, 1971)., is 

I HD ]k = a ;!ji IjW/GCM. 

J 

where A / is the inverse of the matrix. 

A = S V k) I L (k)/G(k). 

k 

(Knowing D(k), H(k) = HD k /D(k). 

A couple of points should be noted. First, the weights pro- 
duced in this manner may be multiplied by an arbitrary scale factor, to give 
a new set of weights; however, the signal to noise ratio will remain un- 
changed. Secondly, if G(k), (the squared noise at a given point) is multiplied 
by a scale factor; (i. e. if all noise terms a^ are increased by the same fac- 
tor), then there is no change to the weights produced. In fact, G(k) may be 
written as: 

r 2 

G(k) = const x 1 + R <I> 

where R = £a 2 /D(k) + a 4 J y/^aj/D(k) + a 3 J 

represents the relative importance of the multiplicative terms with respect 
to the additive terms. If we assume D(k) is constant, then this one param- 
eter will determine the weights produced; whether they optimize against 
multiplicative noises (R>>1), additive noises (R<<1), or some intermediate 
mixture. Analysis of the way R varies with the a’s for reasonable values of 
the noises should indicate which terms will have most affect on the shape of 
the weights produced when varied. 
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Similarly, .it is possible to examine the relative importance of 
the two multiplicative terms, (in the numerator). In fact, the ratio of their 
contribution to the total final noise. is given by a x D/a_. The same is true 
for the two additive terms. On the other hand, to examine the total contri- 
bution of both multiplicative terms in comparison to the contribution of the 
additive terms, it is necessary to first calculate the weights. Once this has 
been done the absolute importance of the four noises may then be calculated, 
in the same units as those in which the target burden is measured: 

N 2 = ai 2;(HD) 2 /D 

N 2 = a 2 E (HD) 2 

N 3 = a 3 E (HD <I> 2 /D ■ 

. N 2 = a 4 S(HD<I>) 2 . . ..... 

6.2.4 Physical Origins for. the Noise Terms : - The various noise 
terms may each arise due to one or more sources. Some possible sources 
are listed below; 

Detector noise. 

Photon noise. 

Other electrical noises. .... . ..... 

Digitization of integrated outputs. 

Scintillation or rapid variations in target, 
albedo or illumination. 

Random errors in derivation of reference signal. 
Random variations in scan waveform. 


RA: - 


RM: - 
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SA: - Presence of spectral signatures which were 

not represented in sample used to derive 
weights. 

Crosstalk from reference into signal channel. 

SM: - Inaccuracies in carrying out multiplications 

(truncation of weights). 

Change in scan waveform from the shape 
when the weights were derived. ' 

6. 3 Software 

6. 3. 1 Description : - As described in detail above, processing of inter- 
ferograms, either theoretical or experimental, consists of applying a linear 
digital filter ("correlation function" or "weights") to the interferogram. 

These weights are chosen so as to give a final result which is indicative of the 
quality of one of the atmospheric components contributing to the interferogram, 
and which is insensitive to variations in all other parameters. The filtering 
operation consists of taking'the weights, one for each point of the interfero- 
gram, multiplying them together with the corresponding interferogram points, 
and summing the results. (The weights are also chosen to give a result which 
is optimized against a mixture of noises.) . The results are also directly con- 
verted using a defined zero point and scale factor, to units of target gas bur - 
den (atm -cm, ppm-m), etc. 

A software package has been developed for use with a large com- 
puter system to input interferograms deriving the weights for a given noise 
mixture, evaluate the noise severities, and to simulate measurement on. inter- 
ferograms using the weights. 

The procedure for the determination of the weights is to* accept 

S_f 

a set of interferograms and define the delay region over which the weights are 
to be applied. One interferogram is designated as "nominal", and another 
designated as "target". Weights are derived, based on the assumption that 
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all interferograms except "target" represent a single value of the desired 
parameter (CO burden), while other parameters are varied. The interfero- 
gram designated as "nominal" is assumed to represent the same values for 
the other parameters as does "target"; i. e. "nominal" and "target" are as- 
sumed to differ only in the size of the target parameter. 

The weights determined will then give measurements which have 
a zero for the target parameter value used in "nominal", and a scale of 1 for 
the target parameter variation between "nominal" and "target". 

(In order to reduce the sizes of numbers to be handled, the nom- 
inal interferogram is subtracted from all others, so that difference interfer- 
ograms are actually used. Because of assumptions of linearity, this makes 
no difference to the algorithm for determining the weights.) 

(It may also be noted that the program actually deals with in- 
terferograms containing both inphase and quadrature phase components. 

They are just handled as if there were twice as many data points I(k); al- 
though an extra index is needed to reference the two components.) 

Values are input to the program representing the noise levels 
(RA, RM, SA, SM), and the total time, T, available for measurement. 

(The measurement time should include any duty cycle factor which may make 
it shorter than elapsed time. Values for RA, etc. will be discussed later.) 

Based on these, the a. are calculated. 

1 

The interferograms and noise parameters are then input to a 
subprogram which carries out the weight generating algorithm, returning the 
weights. 

Provision is included in the program for putting out the noise 
components for given a^ and T, and to also put out the expression relating to 
the size of the noises to the parameters. 

An AGC function is also carried out by the program. In the 
case of computer generated theoretical interferograms, the normalization is 
done on the basis of the interferogram value at the first point, (assumed to 
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be the zero delay point, proportional to the total radiation received). For in- 
terferograms generated by the actual instrument, a good deal of the AGC 
function is carried out by the hardware. However, a term proportional to the 
duration of the measurement remains (number of scans). This factor is re- 
moved for tapes generated by the NOVA by dividing all interferogram points 
by the total number of fringes duration of the measurement at each point, as 
the values are put into core. A flow chart of the program is given in Figure 
6. 3. 1. A listing is given in Appendix B. 

6.3.2 Values for Noise Inputs : - Values for the noise terms may be 
determined by either knowledge or estimation of system parameters, or else 
empirically by the examination of instrument generated interferogram s. 
Conversions from the values input to the a., also depend on whether the inter- 
ferogram set was computer or instrument generated. Certain terms depend 
on whether the program is running with instrument generated (INS) or com- 
puter generated (COM) interferograms. , 

The inputs are RMS noise levels for the following conditions: 

RA - Noise for f unit (1 fringe) integration time (INS) -> , interferogram 
1 unit band width (COM) J units 

RM - Noise for unit integration time, fractional interferogram units 

SA - Noise in limit of long integration normalized to 1 fringe 

integration units 
(INS) 

as a fraction of 
DC level x 10 
(COM) 

SM - Noise in limit of long integration, fractional units 

6. 3. 3 Noise Values Based on Instrument Measurements : - If empir- 
ical observations are available from the instrument the appropriate varia- 
tions, AV, of measured values, V, taken under conditions of N fringes inte- 


94 
















gration per point, the necessary conversions are listed below: (Time, T, 

used with these inputs should be the total duration of the measurement to be 

; - ‘ ( Jt «1 v * *■ r> 

optimized for, in fringes.) 

RA = <AV>//N 

.. RM = < AV/V> x /N 

SA = < A V / V > 

‘ 6. 3. 4 RA and SA f or Computer Interferograms : - In the case of com- 
puter generated interferograms, for the RA term, a system NEP for 1 Hz 

2 

bandwidth may be known, along with the system throughput, T (cm sterad), 

including instrument inefficiencies). The conversion to RA will depend on 

2 

the the inte rferogram units used. These may be either watts/cm sterad, 

or equivalent width of the instrument plus atmosphere transfer function. 

2 

(These correspond to units in the original spectrum of watts/cm sterad 
cm 1 , and fractional transmission, respectively.) The conversions are: 


Units for Original Spectrum 

Detector Limited 

Photon Limited 

2 - 1 

W/cm sterad cm 

NEP/ (T x DNU) 

/ZttN hy/T/DNU 

V 

Fractional Transmission 

NEP/ (T x N^ x DNU) 

/2irhl//N T/DNU 

V 


is the background radiance, in W/cm sterad cm , and DNU is the spec- 
tral interval at which the original spectrum was sampled when the Fourier 
transform was taken. (Provision is included in the program for transferring 
DNU along with the interferograms. If this has been done, it may be set to 1 
in the above conversion.) 
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A rough estimate for the SA term may be made for use with the 

computer interferograms as follows. The mean optical depth, x, may be 

estimated for those absorbing constituents which have not been included in 

the computer modelling, . those constituents which might not be represented 

in the range of observation used to define the weights). Assuming that the 

spectrum of their absorption is random, then the value to be used for SA 

would be SA = X (giving a^ = X times the interferogram DC level). In 

practice, the spectrum is never completely white. If we assume a 1/f type 

of variation, then RA should be decreased according to the distance of the 

interferogram working region from zero, a factor of 10 (implying a distance 
2 

of 10 fringes from the origin) has been built into the program. Thus, for a 
system working N fringes from zero delay, 

SA ~ X/(N/100) 1/2 

This procedure does not work for non-random interfering spectra. If a di- 
rect estimate can be made of the interferogram amplitude due to an unac- 
counted interferent, as a fraction X of the DC level, then the number to use 
is 

SA ~ 10 X. 

(It may be recalled that it is possible to decide on the basis of and a^ 
whether or not the SA noise is significant in comparision to RA noise. ) 

In terms of the inputs, the quantities to consider are listed below. 



Interferogram 1 

Instrument 

Computer 

QA 

RA/SA /DT 

4 RA/SA x WO x DNU x/DT 

QM 

RM/SM /DT 

RM/SM /DT 
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For the computer QA term, W0 is the interferogram DC level, DNU the 
spectrum step if transferred with the interferogram. DT is the duration 
for the measurement of a single point, (fringes for the instrument, seconds 
for the computer interferograms). 

When a Q is much greater than 1, the associated random term 
dominates the synchronous term, and vice versa. 

6. 3. 5 Data Output : - The output of this program is a list of the target 
species burden in both test units of target gas burden. One test unit is the 
difference in species burden between the target and nominal interferograms. 
Since the actual burden is stored along with the interferogram, usually in 
atm-cm, it is possible for the program to convert the test units into burden 
units. 

In addition to applying the correlation function and listing the 
target species burden for the test interferograms, as a check the correlation 
function is also applied to the basic set of interferograms from which the 
correlation function was generated. 
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7. 


RESULTS 

Calculations have been made with several programs. Initial calcula- 
tions were made with a single-line model. For each of several lines of the 
overtone and fundamental bands, calculations were made of absorption across 
the line, the line shape as a function of altitude, net absorption as a function 
of altitude for various CO density profiles (including sinks and various tem- 
perature profiles ,_ integrated net absorption as a function of ground temper- 
ature and emissivity, and of other factors which influence the net intensity. 
With a second program, calculations were made, for the limb mode, of the 
effect of instrument error on the inversion of measured total CO densities 
in the path to obtain a CO profile. With a multi-line program which com- 
putes the spectrum incident on the instrument, generates the corresponding 
interferogram, and, with another program, inverts this to CO densities in 
the path, calculations were made .to determine, the sensitivity of the technique 
and the effects of various atmospheric parameters. - 

7. 1 Atmospheric Transmission 

7.1.1 Calculation Model: - Absorption by the CO overtone (2. 3 p,) and 
fundamental (4. 6 p,) bands has been studied using the computational technique 
described above. The geometry selected for the calculation considers a sat- 
ellite observing the earth from the zenith. Consequently, lines of the funda- 
mental band are strongly absorbed whereas those of the overtone band are 
relatively weak. The results then are based on a comparison between these 
two cases with emphasis placed on the effect of assumed atmospheric pro- 
files on the observed signal. 

Two programs for the calculation of atmospheric spectral trans- 
mission were written. One is a single line model used to calculate atmo- 
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spheric transmission for a variety of conditions appropriate to the CO prob- 
lem. Specifically, it was used to compare the first overtone and fundamen- 
tal bands, to determine the effect of ground temperature variations, to de- 
termine the effect of variation and uncertainty ground emissivities (at spe- 
cific wavelengths), to test the effect of different Lorentz half-widths, and to 
test the effect of low-altitude CO sinks on the transmitted signal. 

The second program is a multi-line model which computes at- 
mospheric transmission for multiple overlapping spectral lines. It has been 
used to calculate the spectrum incident on the instrument, and the output in- 
terferogram. A separate program establishes weighting functions and uses 
them to compute CO density. 

The program can handle a total of 150 different lines, distrib- 
uted in any way between three chemical species (CO plus any two others). 
Three different cases (Figures 4. 1. 1, 4. 1.2) are incorporated in the pro- 
gram: 

1. Ground blackbody radiation (4.6p,) 

r 

2. Limb transmission (2. 3 or 4. 6 |j,) 

3. Reflected sunlight (2.3p,) 

Parameters which can be varied by means of input data include 
the temperature and species' concentration profiles of the atmosphere, ground 
temperature and emissivity, reflectivity, viewing angles, incident radiation 
angle, and grazing height. 

The line overlap system operates basically as follows. For 
each spectral line, a "line inclusion band" is computed. This is the frequen- 
cy band over which each line must be considered. The limits of each band 
are the distance one must go from each line center to reach a threshold inten- 
sity. This threshold intensity is defined as a preselected percentage (usually 
1%) of the weakest line. For the purposes of determining these bands, a 
Lorentz line shape of nominal width 0.06 cm * is assumed. 
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The spectrum is not computed at uniformly spaced points. The 
spacing varies from as low as 0.001 cm * near the line centers to as much 
as 0. 5 cm * between lines. The actual computation of the net transmission 
is done by solving the basic radiative transfer equation in much the same 
way as previously described for the single-line calculations. The overlap of 
lines is represented by addition of their absorption and emission coefficients. 

The complexity of the computation has resulted in a program 
too large to fit a 32 K word computer memory. It was therefore necessary to 
use an overlay system, whereby each section of the program is stored on a 
disc and brought into memory as needed. Using this system, the program 
occupies a maximum of about 24 K words. Preliminary runs using geometry 
case 1 and 13 CO lines had a running time of about two minutes. 

7.1.2 Single- and Multi-line Models and Line Shapes: - The net 
change in intensity for the R7 line of the fundamental is given in Figure 7. 1. 1. 
This shows by the solid curve the results obtained by the single-line model 
and by the.x's the results obtained by the multi-line model. These agree 
within 1% and it can be concluded that the comparison of the fundamental and 
overtone bands can be carried out using the single-line model. 

A variety of atmospheric models were used for the computations 
(AFCRL, 1966). The pressure was taken from the 1966 U. S. Standard At- 
mosphere. Temperature and CO concentration profiles were varied as shown 
in Section 2. For all calculations discussed in the following pages, the con- 
ditions used, if not otherwise stated, are: 

CO Profile 1 

Temperature Model 2 

Ground Temperature = 0 km Atmospheric Temperature 

Emissivity = 1.0 
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Reflectivity = 0. 1 

Bandpass = 20 cm 

Lorentz Half-width (R7) = 0.06 cm * 

Calculations have been made for six CO profile models (Table 2.1.1, Figure 

2 . 1 . 1 ) 

1. Standard (0. 1 ppm) 

2. Low Altitude Urban Source 

3. Low CO (0. 01 ppm) 

4. High Altitude Sink 

5. Low Altitude Sink (Effect to 9 km) 

6. Low Altitude Sink (Effect to 3 km) 

These models were used to evaluate the emission and absorption coefficient 
and compute the CO absorption spectrum which would be observed at an alti- 
tude of 80 km. 

The absorption and emission coefficients are functions of pres- 
sure, temperature, and CO number density. This is illustrated by Figure 

7. 1.2 which shows the spectral absorption coefficient of the 4.6^, P8 line as 
a function of altitude for the standard temperature and CO models. At low 
altitudes the lines are Lorentz broadened and then as altitude increases, 
Doppler broadening dominates. When these absorption coefficients are inte- 
grated over the atmosphere, the line profile variation with altitude is mani- 
fested by an emission peak at the center of the absorption line. This behav- 
ior is illustrated in Figure 7. 1. 3. Physically, earthshine is absorbed up to 
an altitude of about 30 km. At this point, the atmospheric temperature rises 
and the light flux is increased by emission. However, because the line pro- 
files are very sharp at this altitude, the emission is observed only near the 
line center. 

A second example of the importance of atmospheric emission 
is illustrated in Figure 7. 1.4. Here the net change in absorbed intensity 
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Figure 7.1.4' Variation of Rate Change in Absorbed Intensity 
Fundaments (4. 6 (j, ) (R7 Line) 



with altitude as integrated over the P8 line of the fundamental is shown as a 
function of altitude for the four temperature models and the standard CO 
model. Two temperature profiles have inversion layers which cause a net 
increase in the light flux (negative absorption). When the temperature drops 
there is less photon emission and the net absorption increases. For the case 
of the P8 line of the overtone, this behavior is not observed as shown in Fig- 
ure 7. 1. 5. The reason for this is source temperature. In the case of the 
overtone, reflected sunshine dominates the intensity at the earth's surface, 
whereas at 4. 6 p, earthshine is the dominant contributor to the ground inten- 
sity. The atmospheric temperature can be higher than the ground tempera- 
ture and hence atmospheric emission becomes important. 

7. 1.3 CO Profile Effects: - The effect of CO number density profile 
on the intensity has been examined using the standard temperature profile 
and the six CO models. The results for the change of intensity with altitude 
for the P8 line of the fundamental are shown as a function of altitude in Fig- 
ure 7. 1.6. In all cases, except the urban atmosphere, the change in inten- 
sity peaks between 6 to 10 km whereas it should peak at ground level except 
for the low-altitude sink models. This peak is caused by a balance between 
the changing CO concentration and the emission. In the case of the urban at- 
mosphere, the CO concentration profile dominates. This balance should be 
emphasized in that the change in absorption at 4. 6 microns is not directly 
proportional to the CO concentration. On the other hand, similar calcula- 
tions for the P8 line at 2. 3 p, are a direct measure of the CO number density. 
This is shown in Figure 7. 1.7. The difference between the two wavelengths 
is again caused by the effective source temperature. 

These data show that absorption spectra as observed from 80 
km depend on the atmospheric profiles of the absorbing species and on tem- 
perature. The temperature dependence of the spectra is controlled by the 
temperature of the absorbed source. For the case where earthshine is im- 
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Figure 7. 1. 6 Variation of Rate Change in Absorbed Intensity Fundamental (4. 6 |j,) 



portant the absorption profiles are strongly effected by the atmospheric tem- 
perature profile. Also in this case, the change in absorption with altitude is 
not a direct measure of the absorbing species. On the other hand, when the 
absorbed source is the sun, the absorbed signal is almost independent of the 
temperature profile and the change in absorption with respect to altitude is 
directly proportional to the concentration of the absorbing species. 

The results are shown in Table 7. 1. 1 for the six models and in 
Table 7. 1. 2 for models in which concentrations of the CO were multiples 
(.1, .2, .4, .8, 1., 1.6, 3.2, and 6. 4) of the standard atmosphere at all 
altitudes. 

The effect of low altitude sinks on measured CO densities is 
seen in Tables 7. 1.3 and 7. 1.4 using temperature models 2 and 4, respec- 
tively. The numbers presented are those of the percent difference of the ac- 
tual CO density in the model and that calculated using the computed net ab- 
sorption for that model and converting that by Paeans of the data of Table 
7.1.2. It can be seen that at 4. 6 p the errors are larger than the differences 
in the amounts of CO in the standard and the sink models. To improve these 
results, an atmospheric CO profile would have to be known for measured 
CO densities; this cannot be assumed to be obtainable; The errors for the 
temperature model 2 are less drastic than those of model 4 which has an in- 
version layer. The latter errors are extremely large for the 4. 6 p band. 

7. 1.4 Temperature Profile Effects: - Calculations have been made 
for temperature models 1, 2, 3, and 4 (Table 7. 1 . 3 and Figure 7. 1 . 6) using 
the 0. 1 ppm CO model. The calculations shown in Figures 7. 1.4 and 7. 1.5 
were made for the R7, PI, P8, and P26 lines of the fundamental band and the 
R7, PI, P8, and P26 lines of the overtone band, using CO model 1. The 
results for this are shown in Table 7. 1. 5. The tabulated results show a great 
dependence on the atmospheric temperature model for the 4. 6 p band and little 
dependency for the 2. 3 p band. The graphical data of Figures ,7. 1.4 and 7.1.5 
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TABLE 7. I. 1 EFFECT OF CO PROFILE ON ABSORPTION 


FUNDAMENTAL (4. 6p.) 


CO MODEL 

n co |cm ’ 2) 

R7 

PI 

P8 

P26 

STANDARD 

2. 15 (18) 

. 2002 

. 0810 

. 1700 

. 00104 

LOW ALTITUDE SOURCE 

5.21 (19) 

. 6620 

. 2974 

. 5695 

.0167 

LOW CO 

2. 15 (17) 

. 0439 

. 0119 

. 0350 

. 000107 

HIGH ALTITUDE SINK 

2.09 (18) 

. 2006 

. 0807 

. 1703 

. 00102 

LOW ALTITUDE SINK 
(9 km) 

1. 37 (18) 

. 1580 

. 0642 

. 1335 

. 000614 

LOW ALTITUDE SINK 
(3 km) 

1.93 (18) .1 963 

OVERTONE (2. 3 (j. ) 

.0795 

.1 666 

. 000965 

CO MODEL 

"co (cm ’ 2) 

R7 

PI 

P8 

P26 

STANDARD 

2. 15 (18) 

.00647 

. 00144 

. 00462 

. 0000098 

LOW ALTITUDE SOURCE 

5.21 (19) 

t 

. 1226 

. 0299 

. 0953 

. 000431 

LOW CO 

2. 15 (17) 

. 000662 

. 000146 

. 000471 

. 0000011 

HIGH ALTITUDE SINK 

2. 09 (18) 

. 00631 

.00140 

. 00451 

. 0000097 

LOW ALTITUDE SINK 
(9 km) 

1. 37 (18) 

. 00418 

. 00095 

. 00292 

. 0000037 

LOW ALTITUDE SINK 
(3 km) 

1. 93 (18) 

. 00582 

. 00130 

.00413 

. 0000068 
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TABLE 7.1.2 EFFECT OF CO COLUMN DENSITY 
ON FRACTIONAL NET ABSORPTION 


FUNDAMENTAL (4.6 m,) 


n co ,cm ’ i ) 

R7 

PI 

P8 

P26 

2. 15 (17) 

. 0439 

. 0119 

. 0350 

. 000107 

4.3 (17) 

. 0741 

. 0223 

. 0604 

. 00020 

8.6 (17) 

. 1174 

.0405 

. 0978 

. 000410 

1.72 (18) 

. 1769 

. 0692 

. 1496 

. 000820 

2. 15 (18) 

. 2002 

. 0810 

. 1710 

. 00104 

3. 44 (18) 

. 2585 

. 1 108 

.2205 

. 00163 

6. 88 (18) 

.3714 

. 1683 

. 3183 

. 00322 

1.37 (19) 

. 5292 

. 2470 

. 4545 

. 00629 


OVERTONE (2. 

3 M ) 


, -2 





n co (cm • 

R7 

PI 

P8 

P26 

2. 15 (17) 

. 000662 

. 000146 

. 000471 

. 0000011 

4.3 (17) 

. 00131 

. 000286 

. 000926 

. 0000015 

8.6 (17) 

. 00260 

. 000569 

. 00185 

. 0000032 

1.72 (18) 

. 00518 

.00115 

. 00369 

. 0000071 

2.15 (18) 

. 00647 

. 00144 

. 00462 

. 0000099 

3.44 (18) 

. 0102 

. 00228 

. 00729 

. 000014 

6. 88 (18) 

. 0197 

. 00455 

. 0142 

. 000030 

1.37 (19) 

. 0370 

. 00895 

. 0272 

. 000059 



TABLE 7. 1. 3 f EFFECT OF CO PROFILES ON DERIVED 
CO CONCENTRATIONS (% DEVIATION) 


FUNDAMENTAL (4. 6 u ) 


CO MODEL 

n CO (cm ’ 1 

R7 

PI 

P8 

P26 

HIGH ALTITUDE SINK 

2.09 (18) 

0.2 

0 . 0 

8. 3 

2.5 

LOW ALTITUDE SINK (9 km) 

1. 37 (18) 

14.6 

10. 5 

18. 0 

4.9 

LOW ALTITUDE SINK (3 km) 

1.93 (18) 

8. 1 

2.7 

12.7 

1.0 

• ! ■ • ' ■ l 

CO MODEL 

OVERTONE 

(2.3n) 




n co (cm ’ 2) 

R7 

PI 

P8 


HIGH ALTITUDE SINK 

2.09 (18) 

3.5 

4. 5 

3. 5 


LOW ALTITUDE SINK (9 km) 

1.37 (18) 

6. 1 

6. 1 

6. 1 


LOW ALTITUDE SINK (3 km) 

1.93 (18) 

2.7 

2.7 

2.7 
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TABLE 7. T. 4 EFFECT OF PROFILES ON DERIVED 
• CO CONCENTRATIONS - TEMPERA- 
TURE INVERSION (% DEVIATION) 


• FUNDAMENTAL (4. 6 n) 


' CO MODEL l '"' : 

n co (cm ’ 2) 

R7 

PI 

' P8 • 

P26 

LOW ALTITUDE SINK (9 km) 

• I. 37 (18) 

13. 3 

10.4 

17.0 . 

0.04 

LOW ALTITUDE SINK (3 km) 

1.93 (18) ' 

0.9 

29. 5 

32.6 

35.2 


OVERTONE 

(2. 3 p.) 




CO MODEL 

”co ,cm ' 2) 

R7 

PI 

P8 


LOW ALTITUDE SINK (9 km) 

1.37 (18) 

7.4 

7.4 

7.4 


LOW ALTITUDE SINK (3 km) 

1.93 (18) 

3.6 

4. 7 

z: 6 




TABLE 7.1.5 EFFECT OF TEMPERATURE 
PROFILE ON ABSORPTION 


4 . 6 (j, 

R7 

PI 

P8 

P26 

LOW TEMPERATURE 

. 1613 

. 0677 

. 1328 

. 000489 

NORMAL TEMPERATURE 

.2002 

. 0810 

. 1700 

00104 

HIGH TEMPERATURE 

. 1942 

. 0772 

. 1662 

.00134 

inversion Layer 

. 127b 

. 0521 

. 1063 

. 000433 


2 . 3 |jl 


LOW TEMPERATURE 

. 00660 

. 00152 

. 0b462 

.0000062 

NORMAL TEMPERATURE 

.00647 

. 00144 

. 00462 

. 0000098 

HIGH TEMPERATURE 

. 00637 

. 00140 

. 00461 

.000014 

INVERSION LAYER 

. 00639 

. 00140 

. 00461 

. 000013 
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show the contribution of various altitudes to the net absorption. The 2. 3 p, 
band shows absorption as a function of altitude to be as would be expected 
whereas the 4. 6 |i band shows much less absorption below about 5 km (and, 
in fact, a significant net emission in some regions) which is an important 
region for CO sink studies. 


7.1.5 Ground Temperature Effects : - The calculations described 
above included a ground tejnperature of each model at 0 km. Calculations 
have also been made for a different ground temperature. The following cases 


were carried out. 


T 


Ground 


Atm. T. Model 


T 


Model Atm 




288.2 2 

288.2 1 

288.2 3 

288.2 . 4 

273.2 2 

284.2 2 

287.2 2 

289.2 2 

291.2 2 

308.2 2 


288. 2 
257. 3 
302. 6 
283. 2 
288.2 
288. 2 
288. 2 
288. 2 
288. 2 
288. 2 


The results of these calculations are shown in Tables 7. 1. 6 and 7. 1. 7 and in 
Figure 7. 1.8. Drastic effects are seen for the fundamental band <4.6 (jl) and 
negligible effects are seen for 2. 3 pu . From Figure 7.1.8 for the R7 line of 
the fundamental band it can be seen that an underestimate of the ground emisr 
sion by 14% (as would be caused by an underestimate of the ground tempera- 
ture by 4°) causes a 12% overestimate in the CO density while a 12% over- 
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TABLE 7.1. 6 •, EFFECT. OF. GROUND TEMPERATURE » 



•• ON ABSORPTION 

■' ' <! v . • V . 




FUNDAMENTAL (4. 

6 p, ) 





R7 

PI 

P8 

P26 


257.3 

1613 

. 0677 

. 1328 

. 000489 

LOW TEMPERATURE 

288.2 

2500 

.1030 

.2117 

. 00104 


302.6 

1942 

. 0772 

. 1662 

. 00134 

HIGH TEMPERATURE 

288.2 

1342 

.0550 

. 1128 

. 00415 


283.2 

1270 

.0521 

. 1063 

. 000433 

INVERSION LAYER 







288.2 . 

1540 

. 0622 

.1302 

. 000793 


OVERTONE 

<2. 3ji 

) 





R7 

PI 

P8 

P26 


257.3 

00660 

. 00152 

. 00462 

. 0000061 

LOW TEMPERATURE 

288.2 

00657 

. 00150 

. 00459 

. 0000053 


302.6 

00637 

. 00140 

. 00461 

. 0000137 

HIGH TEMPERATURE 

288.2 

00634 

. 00138 

. 00457 

. 0000122 


283.2 

00639 

. 00140 

. 00461 

. 0000128 

INVERSION LAYER 

288.2 

00635 

. 00139 

. 00458 

. 0000144 
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TABLE 7.1.7 EFFECT OF GROUND TEMPERATURE 
ON APPARENT ABSORPTION 
OF SOURCE RADIATION 


FUNDAMENTAL 

<4.6n) 


R7 

273. 2 

. 1337 

284.2 

. 1866 

287. 2 

. 1971 

288. 2 

. 2002 

289.2 

.2033 

291.2 

1 

.2091 


308.2 


2428 
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Figure 7. 1.8 Effect of Change in Earth Emission Intensity- 
on Predicted CO Density (4. 6 p - R7 Line) 






estimate of the earth emission (a three degree overestimate of the ground 
temperature) would cause about an 8% underestimate in CO density. These 
calculations all used the same temperature profile (2) and the same CO pro- 
file (1), varying the earth emission. When searching for sinks which cause 
changes in CO density of the order of 10% such errors are larger than can be 
tolerated. Other lines of the fundamental show similar effects. The curve 
for the R7 and other lines of the overtone would lie on the abscissa showing 
that there is no ground intensity effect on this band. 

7. 1.6 Emissivity Effects : - In calculations described above, an emis- 
sivity of 1 was used. Since certain ground areas have appreciably different 
emissivities , calculations have been made using emissivity values of 0.95, 
0.9, 0.8, and 0.7 for the 4.6 p band and a value of 0.7 for the 2. 3 p band, 
with each of the four temperature models. 

Drastic effects in the fractional net apparent absorption are 
seen by the data in Table 7. 1 . 8 at 4. 6 p, while those at 2. 3 g, are negligible. 
Figure 7.1.8 shows that a change (AI^) in ground emissivity from 1 to . 9 to 
.8 gives changes in measured CO density of 8% and 18%, respectively. CO 
density is obtained from intensity changes by use of data in Table 7. 1.2. 

It can be seen that atmosphere and ground temperatures and 
ground emissivities would have to be very accurately measured to be able to 
interpret any CO data obtained by remote measurements at 4.6y,. It appears 
that the best expected temperature data will have an accuracy of ± 2 K 
(Houghton and Smith, 1970; Abel, et al, 1970a; Abel, et al, 1970b). These 
data would not be expected to be this accurate if an inversion layer exists. 
Even ± 2K would present substantial effects at 4. 6 (j, . Further, the ground 
temperature and emissivity would have to be measured accurately and would 
have to be known to be that of 4.6y,, since the emissivity may vary with 
wavelength. 
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TABLE 7. 1. 8 EFFECTS OF EMISSIVITY ON 

FRACTIONAL NET ABSORPTION 
FOR R7 LINES 


TEMPERATURE MODEL • 


LOW 

TEMPERATURE 


€ 4. 6 jo, 2. 3 (j, 

1 .161 .00660 

0. 95 .155 

0.9 . 149 ' 

0. 8 .133 

0.7 .113 .00657 


1 .200 

NORMAL 0.95 • 196 

TEMPERATURE °- 9 * 191 

0. 8 . 180 

.. 0. 7 .165 


.00647 


. 00644 


NORMAL 1 .200 

TEMPERATURE 0.9 .191 

Tg =. 288. 2 0.7 .. 165 


HIGH 

TEMPERATURE • 


1 

. 194 

. 00637 

0. 95 

.190 : 


0.9 

. 184 


0. 8 

. 172 


0. 7 

. 156 

. 00634 


INVERSION 

LAYER 


1 . 127 . 00639 

0. 95 . 119 

0.9 .110 

0. 8 . 0886 

0.7 .0710 


. 00635 



7. 1.7 Effects of Other Parameters :. - Calculations made to test the 
effect of reflectivity ; show this effect to be small for the 4. 6 g, band and negli- 
gible for the 2. 3 g, band. These data are given in Table 7. 1.9. The effect is 
due to a change in source intensity with no accompanying change in atmo- 
spheric emission. 

Calculations made to show the effect of Lorentz half-width on 

absorption. Significant effects were found for some conditions for the 4. 6 g, 

band while the effects on the 2. 3 g, band were small. Results are shown in 

2 

Tables 7. 1. 10 and 7. 1. 11. This shows that the value of a in Equation 4.21 

' ■' » 14 

must be known to accuracies of the order of 10% for the fundamental, while 
for the overtone an accuracy of 50% is sufficient. 

Calculations to test the effect of bandpass (the distance from the 
center of the line over which the absorption is integrated) were'made, As 
expected this has an appreciable effect if too small a bandpass is taken. 

Data for four lines of the 4. 6 and 2. 3 g, bands are given in Tables 7. 1. 12 and 
7. 1. 13 for bandpasses (v = V q ) of 20, 2, and 0. 1 cm * . 

7.2. Expansion of Overtone and Fundamental Bands 

Using the data given in the preceeding tables, a summary may 
be made of the two general wavelength regions, the 2. 3 p, band of the first 
overtone of CO and the 4. 6 g, band of" the fundamental of CO. All other CO 
bands have been found to be impractical (Bortner and Kummler, 1971). 

The absorption of the 2. 3 g, radiation has been found to be suffi- 
cient to produce the required sensitivities and ranges using the correlation 
interferometer technique for optical thicknesses appropriate to both the map- 
ping and the limb experiments, considering nominal ambient atmospheric con- 
centrations arid reasonable fractions thereof. 

The absorption of the 4. 6 g, radiation is, of course, sufficient for 
the experiments. One difficulty arises at the higher optical thicknesses ap- 
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TABLE 7. 1. 9 EFFECT OF REFLECTIVITY 
ON FRACTIONAL NET 
APPARENT ABSORPTION 



P 

R7 


0.4 

.2058 

FUNDAMENTAL 
(4. 6 u) 

0 . 1 

.2002 


0.02 

. 1986 


0.4 

.00651 

OVERTONE 

(2.3U) 

0.1 

.00647 


0.02 

.00655 



TABLE 7. 1. 10 EFFECT OF LORENTZ HALF - WIDTH 
ON APPARENT ABSORPTION OF 
FUNDAMENTAL (4. 6 p,) 




R7 

PI 

P8 

P26 

STANDARD CO 

.06 

. 2002 

.0810 

.1700 

.0010 


.09 

.2387 

.0905 

.2012 

.0010 

LOW ALTITUDE 

.06 

.6620 

.2974 

.5695 

.0167 

SOURCE 

.09 

.8052 

. 3555 

.6921 

.0178 

LOW CO 

.06 

.0439 

.0119 

.0350 

.0001 


.09 

.0470 

.0121 

.0370 

.0001 

HIGH ALTITUDE 

.06 

.2006 

.0807 

. 1703 

.0010 

SINK 

.09 

.2390 

.0898 

.2014 

.0010 

TABLE 7.1.11 

EFFECT OF LORENTZ HALF-WIDTH 


ON 

ABSORPTION OVERTONE (2. 

3u) 



R7 

PI 

P8 

P26 

STANDARD CO 

.06 

.00647 

.00144 

.00462 

.0000098 


.09 

.00647 

.00142 

.00461 

.0000085 

LOW ALTITUDE 

.06 

. 1226 

.0299 

.0953 

.000431 

SOURCE 

.09 

.1312 

.0303 

.1005 

.000420 

LOW CO 

.06 

.000662 

.000146 

.000471 

.0000011 


.09 

.000649 

.000140 

.000457 

.0000006 

HIGH ALTITUDE 

.06 

.00631 

.00140 

.00451 

.0000097 

SINK 

.09 

.00631 

.00139 

.00450 

.0000083 



TAELE 7.1. 12 


EFFECT OF BANDPASS ON FRACTIONAL 
NET APPARENT ABSORPTION 


FUNDAMENTAL (4.6^) 




R7 

PI 

P8 

P26 


2° 

. 1613 

. 0677 

.1328 

. 00049 

LOW 

TEMPERATURE 

2 

. .1593 

.0670 

. 1314 

. 00047 


0.1 

. 1086 

. 0142 

.0951 

. 00042 


20 

. 2002 

. 0810 

. 1700 

. 00104 

NORMAL 

TEMPERATURE 

2 

. 1979 

•• 0804 

. 1682 

. 00102 


0. 1 

. 1286 

.0635 

. 1155 

. 00093 


20 

. 1942 

. 07J2 

.1662 

. 00134 

HIGH 

TEMPERATURE 

2 

. 1921 

. 0766 

. 1646 

. Q0 1 33 


0. 1 

. 1264 

. 0610 

.1138 

. 00120 


20 

. 1270 

. 0521 

. 1063 

. 00043 

INVERSION 

LAYER 

2 

. . 1260 

. 0518 

. 1056 

. 00041 


0. 1 

. 0977 

. 0450 

. 0852 

. 00040 



TABLE 7. 1 . 1 3 EFFECT OF BANDPASS ON FRACTIONAL 
'NET' APPARENT ABSORPTION 




OVERTONE 

(2.3 m, ) 





R7 

PI 

P8 

P26 

■ ■ ’ 

20 

. 00660 

00152 

. 00462 

. 0000062 

LOW 

TEMPERATURE 

2 

.00656 

.00150 

. 00459 

.0000053 


0. r 

.00544 

. 00125 

. 00383 

. 0000053 

1 'i . * 

20 i 

.f . 00647 

.00144 

. 00462 

. 0000098 

NORMAL 

TEMPERATURE 

2 

,. 00643 

.00143 

. 00459 

.0000089 

■ i . 1 

0. 1 

; . . 00531 

. 00118 

. 00381 

.0000087 

r :■ ... , . • • 

20 

.00637 

. 00140 

. 00461 

. 000014 

HIGH 

TEMPERATURE 

2 

.00633 

.00138 

. 00457 

. 000012 


0. 1 

.00521 . 

. 00114 

. 00378 

. 000012 


20 

. 00639 

. 00140 

. 00461 

. 000013 

INVERSION 
LAYER .... 

2 

. 00635 

. 00139 

. 00458 

. 000012 


0. 1 

. 00523 

. 00115 

. 00378 

. 000011 
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propriate to the lower altitude limb measurements. Under such conditions 
the absorption is so strong that it is well off the linear portion of the curve- 
of-growth with ambient atmospheric CO concentrations resulting in a loss of 
sensitivity in the CO measurements. 

The data show that there is considerably more error in mea- 
surement at 4. 6 |jl of the CO column density in the important case of a low al- 
titude sink than for the same measurement at 2. 3 p. Since a major objective 
of the experiment is to search for a potential low altitude sink, this indicates 
a definite advantage for 2. 3 p. In the case of a temperature inversion drastic 
errors are found in the case of a low altitude sink. Since temperature in- 
versions are common the 2. 3 jjl has another distinct advantage. In general it 
can be said that the 4. 6 p is greatly affected by the temperature profile. 

This is as expected due to contribution of atmospheric emission. The 2. 3 p 
radiation is essentially unaffected, as is desired. Similarly the ground tem- 
perature and the ground emissivity, have larger effects (not desired) on the 
CO absorption signal at 4. 6p but not at 2. 3 p. 

The calculations show that, assuming measurement of CO ab- 
sorption can be made on the 2. 3 p band, that band is much to be preferred 
over the 4. 6 p band for data interpretation. The advantages are: 

1. Measurements of temperature of the atmos- 
phere, ground temperature, and emissivity 
are not needed for 2. 3 but accurate values 
are required for 4. 6 p. 

. t 

2. The signal at 2. 3 JJ. is affected significantly, 
as desired by the CO concentration in the low- 
est few kilometers in the atmosphere, but not 
significantly at 4. 6 u and in a manner very 
difficult to interpret. 
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3 . 


The measurements at 2. 3 g are readily in- 
terpreted directly in terms of CO densities 
whereas those at 4. 6 g must be interpreted 
by use of an atmospheric model' calculation 
at each point, assuming the required atmo- 
spheric data are available for the latter. 

7. 3 Comparative Calculations 

Calculations given above of predicted 4. 6 g CO band transmis- 
sion showed large effects caused by inversion layers, ground temperatures 
and emissivities , and little effect of low-altitude sinks. It was suggested to 
us that other calculations (Ludwig, 1970) did not show effects of such magni- 
tude. Further calculations were made to make a reasonably direct compar- 
ison of the computational methods. Two specific comparisons were carried 
out. One calculation was that of the transmission of a single line, the P19 
line. The second was for several lines to be compared with the results of a 
band model calculation. 

The single-line calculation was carried out for four tempera- 
ture models, shown in Table 2. 1.4 as models 5, 6, 7, and 8, and a 0. 2 ppm 
CO model. It is not known what exact temperature model was used in the 
calculation to which the present ones are being compared, but it is believed 
that it is in the range covered by these models; the authors reporting (Lud- 
wig, 1970, p22) those calculations merely state that they used a "represen- 
tative temperature model". The results of the calculations are shown in 
Figure 7. 3. 1, which plots the transmission over the line profile (The cir- 
cles showing the results of Ludwig (1970), the remaining symbols being our 
results for the four models.). It can be seen that there is excellent agree- 
ment between the two calculations. The only visible difference is at the line 
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center where the present calculations yield a 19 to 20% transmission while 
the calculations being compared show a 21% transmission. This small dif- 
ference may be due to differences in temperature models. From . 002 cm * 
from the center outward, the two calculations show no appreciable .difference. 
The integrated absorption over the line would certainly show no meaningful 
variation. 

The other comparison which is informative is one showing the 
effects, of a temperature inversion and of ground temperature variation. 

Data (Ludwig, 1970, p9) are given as models 9. 10, and 1 1 in Table 7. 3. 1. 
Model 10 shows the effect of a ground temperature different from the lowest 
atmospheric temperature, the results (Ludwig, 1970) showing a calculated 
signal change variation over the band of 50 and 68% for . 25 and . 025 ppm CO 

* / 

atmospheres, respectively. These can be compared with previously pre- 
sented results, Table 7. 1.7, of our calculations which show a 33% decrease 
: in absorption for a single line (R7) of the same band. Model 11 shows the 
effect of a temperature inversion layer giving variations of signal change 
over the band of 42 and 55% for the . 25 and . 025 ppm CO atmospheres, re- 
spectively (Ludwig, 1970). These may be compared w,ith our results, Table 
7, 1.5, showing changes of 37, 36, 37, and 76% for the R7, PI, P8, and P26 
lines of that band, respectively. 

For one case of interest, a low altitude sink model. Table 
7. 1.4. , no data were given on the other calculations (Ludwig, 1970) to com- 
pare with ours . 

It can also be noted that a graphical representation of the effect 
of ground temperature differences on CO calculations, Table 7. 1.8. , can be 
roughly compared to those of Ludwig (1970, pl6) and seen to give reasonable 
agreement. 

, It can be concluded that the two programs give the same results, 

at least as nearly as we can determine with available information. 
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TABLE 7. 3. 1 


TEMPERATURE MODELS AND RESULTS 
FOR COMPARISON BAND CALCULATIONS 


Model 9 ■ 10 11 Reference 


Signal Change 

.25 ppm CO 7. 14 3.59 4. 16 , a 

.025 ppm CO 1.35 0.43 0.61 a 

Relative Signal Change 

.25 ppm CO 1 0.50 0.58 a 

. 1 ppm CO 1 0,33 0,46 b 

.025 ppm CO 1 0.32 0.45 a 


a) Ludwig (1970) 

b) This work (Tables 7. 1. 7 and 7. 1. 8) 



. 4 


Limb Inversion Analysis 


The limb experiment has been analyzed for the effects of in- 
strumental error on the inversion of the instrument values of total CO in the 
path to yield CO concentration vs altitude. The calculations have been per- 
formed for the three CO sink models to which the limb experiment is rele- 
vant. These are the constant mixing ratio model with 0. 1 ppm at ground 
level), the upper atmosphere sink, and the low-level sink with maximum con- 
centration at 9 km. The error sources, which have been introduced, include 
a random error with standard deviation of 2%, 5%, or 10% alnd systematic 
errors of ± 2%, ± 5%, or ± 10%. 


For these calculations the atmosphere below an altitude, h^, 
was divided into n spherical shells (j = l,2,---n), a distance, Ah, apart 
(Figure 1 : 4. 1). The atmospheric properties within a shell are assumed to 
be constant. The altitude of the first shell is given by 


h j = n Ah, 

and the altitude of the jth shell is 

h. = hj - (j-1) Ah 

The ith ray from the sun to the instrument passes horizontally 
through the lower edge of the ith shell (j = i). The height of this ray is 


h.. = h. - Ah 
n i 

The rays passing through the shells can be divided into lengths, 
a^j, over which the atmospheric properties are assumed to be constant. The 
mgths a., are given by: 
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1/2 


{(r + h.f - 

<r * »..) } 

if i = j 

{(r + h., 2 - 

2l' /2 

<- ♦ v } 


k = i 

- ZZ 

k = j + 1 

a ik 

if i ± j 


A square matrix (n x n) of the elements a., can be formed. 

In a shell j, the average concentration of GO is giveh by: 

c . = c(h. - 1/2 Ah) 

J J 


where 


The integrated amount of gas in a horizontal path is given by: 


M = 2 { c jl [ a ij] 



column matrix of integrated amount of CO 

column matrix of the average concentration 
of CO 

square matrix of path length elements 
within the shells; 


The calculations were carried out with a 40 shell atmospheric 
model with a 2 km distance between shells. The matrix elements' u. were 

l 

calculated for the three CO models in Figure 7*4.2 (Models 1, 4, and 5 of 
Figure 2. 1. 1 and Table 2. 1. 1); The correlation interferometer gives an in- 
strument output which is proportional to the total CO ihi the path, or u^» In 
this case the known parameters will bfe u as a function' of altitude, and the 
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a elements will be known. A set of equations can then be set up to solve 

ij 

for the concentrations c.. 

1 

{ C J = 1/2 { u jl [ a ij]'' 

-1 

where £a..J is the inverted matrix £a_J . 

To test for the effects of random oi* systematic instrument 
errors, the values of u were computed for each of the three concentration 
models from the relation j^u.j- = 2-|c.j- These values are shown in 

Figure 7.4. 3. The computed values were then changed by fixed amounts 

a) random error with 2%, 5%, or 10% standard 
deviation. 

b) systematic error of ± 2%, ± 5%, or ± 10%. 

For each of these types of errors, the concentration profile was calculated 
.and compared with the exact profile used as the original input. The results 
of the calculations are shown in Figures 7.4.4 and 7.4.5 and in Tables 
7.4. 1, 7.4. 2, and 7. 4. 3. ' 

Figure 7.4.4 shows the result of imposing random errors with 
standard deviations of 5% and 10% on the measurements of u. With measure- 
ment errors of this magnitude the computations show that it is possible to 
distinguish the carbon monoxide constant mixing ratio concentration profile 
from the concentration profile which assumes a high altitude sink. The com- 
putations with a 2% random error show the same result with a much smaller 
error. These data were not plotted to avoid cluttering the graph with too 
many points. The data are tabulated in Tables 7. 4. 1 , 7.4.2, and 7.4.3 for 
the 2% random error calculations. 

Figure 7.4.5 shows the result of performing the data inversion 
to obtain the concentration profile, when the instrument measurement of the 
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Figure 7.4.4 Computed Atmospheric CO Concentrations with Simulated Instrument Error 



Figure 7. 4. 5 Computed Atmospheric CO Concentrations with Simulated Instrument Error 



TABLE 7.4. 1 LIMB TRANSMISSION ANALYSIS FOR 

STANDARD CO CONCENTRATION MODEL 


Altitude 

(km) 

Shell 

No. 

78.3 

43*3 

7 6.0 

39.0 

74.0 

38.0 

72.0 

37.3 

70.0 

36.0 

68.0 

3 5. '3 

66.0 

34.0 

64.0 

33.0 

62 • 3 

32.0 

60.0 

31 .0 

58.0 

30.3 

56.0 

29.0 

54. 0 

28.0 

52*0 

27.0 

50.0 

26.0 

48.0 

25.0 

4 6*0 

2 4.0 

4 4.0 

23.3 

42.0 

22.0 

40.0 

21 .0 

38.0 

20.0 

3 6.0 

19.0 

34-0 

18-0 

32.0 

17.0 

30.0 

1 6-0 

28*0 

15*3 

2 6.3 

14.0 

24.3 

13*0 

22.0 

12.0 

20.0 

11.0 

1 8.0 

1.0.0 

1 6.0 

9.0 

14.3 

8.0 

12.0 

7.3 

10.0 

6.0 

8.0 

5-3 

6.3 

4.3 

4.3 

3.0 

2.0 

2.0 

3. 

1.3 


No Error 
Mean CO Cone. 
(Part, / cm 3) 

4. 8S5E+0 7 

6. 673E+37 
9 • 0 1 4 E+ 3 7 

1 • 23 3 E+3 3 . 

I • 593 £>3 8 

2.079E+3S 

2. 695E+33 
3.46SE+3S 
4. 420 E+3 8 
5- 620 E+38 

7. 1 9 5E+33 
9. 1 72E+0B 
1 . 1 66E+09 
1 . 477E+39 

1 • 887E+ 39 

2.41 7 E> 0 9 
3.111 E+09 

4.0 89 E+39 
5.432E+09 
7. 189 E+39 

9. 620 E+09 
1 .296E+10 
1 . 763 E+ 1 0 

2 • 40 5E+ 1 0 
3.283E+10 

4. 4 66E+ 10 
6.094E+ 13 . 
8* 335F+ 1 0 

1 . 1 44E+1 1 
1 . 5 7 4 E+ 1 1 

2. 1 63E+1 1 
2.9 58 E+ 1 1 
4*0 50 E+ 1 1 

5. 542 E+ 1 1 
7. 5S6E+ 1 1 

9.711 E+ 1 1 

1 • 22 7 E+ 1 2 
1 . 531 E+12 
1 . 8 91 E+ 12 

2 • 3 1 2 E+ 1 2 


2% Random 
Mean GO Cone, 
(Part. / cm^) 

4* 88 6 £+13 7 

6. 60 6E+0 7 
9. 878E+07 
1 . 10 5E+3 8 
1 . 62 6 E+3 8 

2.31 7 E+38 
2 . 655E+38 
3* 423E+08 

4. 2 53 E+3 8 

6.3 43 E+38 

7. 1 08 E+38 
9*382 E+38 
1 . I 46E+39 
1 . 53 6 E+09 
1 . 79 6 E+09 

2. 383 E+09 

3.0 51 E+«9 
4*031 E+39 

5. 1 7 6 E+09 
7. 7 65E+09 

1 . 3 1 2 E+ 1 0 
1 .294E+13 
1 . 809 E+10 
2.266E+10 
3*31 7E+ 1 0 

4.292E+13 

6. 3 80 E+ 10 
8.282E+ 10 
1 . 1 42E+ 1 1 

1 . 579E+ 1 1 

2 • 334E+ 1 1 
2. 874E+ 1 1 

4. 1 92E+ l 1 

5. 706E+ 1 1 

7. 687E+ 1 1 

9. 867E+ l 1 
1 • 20 5E+ 1 2 
1 • 522E+ 1 2 

1 . 932 E+ 1 2 

2 • 4 1 4 E+ 1 2 


+ 5% Constant 
Mean CO Cone 
(Part, / cm 3) 

5. 1 29E+07 
7.007E+07 
9 . 465E+37 
1 • 2 63 E+ 33 
1 .673 E+3 8 

2. 1 S3 E+ 3 8 
2*833 E+3S 
3. 638E+3 8 

4. 641 E+38 

5. 931 E+0S 

7. 554E+08 
9. 631 E+3 8 
1 .224 E+09 
1 . 551 E+09 
1 • 9 82 E+09 

2. 538 E+39 
3*266 E+39 
4-293E+09 
5. 673 E+09 
7. 546E+09 

1.01 0 E+ 1 0 
1 *361 E+10 
1 -848 E+10 

2. 52 5E+13 
3 • 4 4 7 E+ 1 0 

4* 690 E+ 1 3 
6.39 9E+ 10 
3 • 7 52 E+ 1 0 
1 .20 1 E+ 1 1 
1 . 6 53 E+ 1 1 

2 . 2 7 1 E+ 1 1 

3. 1 0 6E+ 1 1 

4. 2 52 E+ 1 l 
5.S19E+1 1 
7.965E+1 1 

1 • 323 E+ 1 2 
1 • 288E+ 1 2 
1 . 60 8 E+ l 2 

1 . 9 8 6 E+ 1 2 

2 • 42 8 E+ 1 2 



TABLE. 7. 4.2 LI1&B TRANSMISSION ANALYSIS 
FOR UPPER ATMOSPHERE SINK 
CO CONCENTRATION MODEL 


No Error 


Altitude 

Shell 

Mean CO Cone. 

M 

No. 

(Part. / cm^) 

78* O'- 

40 .0 . 

4 . 8 8 5 E+ 3 6 

7 6*0 

39.0 

6. 673 E+0 6 

7,4*3 

38.0 . ... 

9.0 1 4 E+0 6 

92*0;, 

3 7.0 . 

, 1 .203E+07 

70.0 

36.3 

1.590 E+0 7 

. ,6,8 *0 

3 5*3 

2 .0 79 E+0 7 

66* 0 

34. 0 . 

2. 69 5E+0 7- 

64*0 

33*0 

. 3.4 65E+0 7 

62*0 

32.0 

4* 420 E+0 7 

63*0 

31.0 

5* 620 E+0 7 

.58*0 

33 .0 

7.1 9 5 E+0 7 

56*3 

29*3 

9. 172 E+0 7 

54*0 

28.0., ' 

■ 1.166E+38 

52*0 

27.0 

1.47 7 E+08 

53*0 

2 6.0 

1 • 887E+38 

48*0 

25-0 

2-41 7 E+0 8 

46*0 

24-3 

3-111 E+38 

44*0 

23.0 

4. 089 E+08 

42*3 

22.0, 

, 5. 799 E+08. 

40 * 0 

21.0 

8. 401 E+08 

■ 1 38*3 

20.0' 

; l . 2 40 E+39 , 

3 6*3 

19.0 

1.881 E+39 

34*0 

18.0 

. 2. 930 E+ 09, 

32*0 

1 7.0 

4. 651 E+0 9 

30*3 

1 6.0 

7. 59 8 E+0 9 

28.0 

1 5.0 

1 .293E+10 

2 6*3 

14.0 

2 • 3 1 9 E+ 1 0 

24* 0 

13-0 

4 . 1 62 E+ 1 0 

22*0 . 

12*0 

,7. 4 53 E+l 3 

. 20.0 ■ 

11.0 

r.320E+tl 

18.0 

10.3- 

2 • 1 6 3 E+ 1 1 

1 6.0 

9.0 

2. 9 58 E+l 1 

1 4.0 

8.0 • 

4. 0 49 E+ l 1 

12.0 

7.0 

'5’« 5 43 E+ 1 1 

' 13 : * 0 

6.0 

7.58 6E+ 1 1 

8*0 \ 

5-0 

9 . 7 1 .3 E+ 1 l 

6-0 

4-0 

1.22 7 E+l 2 

. 4. 3 . 

3*0 

1 *‘531 E+ 1 2 • 

2 • 0 

2.0 - 

i 1 .889 E+ 12 - 

0 • , . * 

1*3 

4 2 .3 13 E+ 1 2 


2% Random +5% Constant 

Mean CO Cone. Mean CO Cone. 

(Part. /cm 3 ) (Part, /cm 3 ) 

4.9 51 E+0 6 5. 129 E+0 6 

6.449R+06 7 .00 7 E+0 6 

9.025E+06 9.465E+.0 6 

.. 1.254E+07,. 1 * 2 6 3 E+ 0 7 

1 . 59 1 E+0 7 1 • 6 70 E+0 7 

f> • 1 2 5E+ 0 7 2.183E+3 7. 

.2 • 80 3 E+0 7 2. 833 £+0 7 

3* 2 54E+3 7 3. 638 E+0 7 

4. 62 3 E+0 7 4 • 6 4 .1 E+ 0 7 

5. 59 1 E+0 7 . 5*901 E+07 

6* 8 57 E+0 7 '7. 5 54 E+0 7 

8*7 51 E+0 7, 9* 631 E+0 7 

1 * 22 5E+ 38 1 * 2 2.4 E+ 3 8 

1* 49 6E+08. r 1*5 51 .E+0 3 

1 • 9 68 E+0 8 1*982 E+0 8 

2* 38 7 E+0 8 2* 538 E+0 8 

3 *09 8 E+0 8 3*2 668+08 

3*99 3 E+0 8 4*293 E+08 

5* 7 51 E+0 8. ", 6* 0 8.9 E+0 8 

8*01 6 E+0 8 8*821 E+08 

1*2 08 E+0 9 1 * 302 E+0 9 

1 * 8 78 E+0 9 1 • 9 7 5 E+0 9 

2*931 E+0 9 3*3 76E+09 

4. 73 4 E+ 0 9 4 * 8 84 E+0 9 

7* 71 5 E+39 7 * 9 78 E+0 9 

1 *203 E+ 10 1 * 3 5.4 E+ 10 

2* 2 63 E+ 13 2.43 5E+10 

4.227E+10 4.370E+10 

7* 223 E+ 10 7 • 822 E+ 1 0 

1 * 308E+ 11 1 * 38 6E+ 1 1 

2*1 63E+1 1 2*2 7.1 E+l 1 

2*871 E+l 1 3* 1.0 6E+1 1 

... 4* 1 73E+ 1 1 4* 2 52 E+l 1 

5*271 E+l 1 \ 5* 823 8+11 

7* 9 82 E+ 11 7*96 5E+ 1 1 

9 * 3 47 E+ 1 1 1*0 23 E+ 1 2 X 

1 • 1 72 E+ l 2 1 * 288E+ 1 2 

1 * 52 1 E+ 1 2 1 * 60 8 E+ 1 2 

1 * 891 E+l 2 1 * 9 8 4 E+ 1 2 

2 * 2 1 2 E+ 1 2 2 * 428 E+ 1 2 
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TABLE 7. 4. 3 LIMB TRANSMISSION ANALYSIS FOR LOW 

Level sink, co concentration model 




No Error 

2% Random 

+ 5% Constant 

Altitude 

Shell 

Mean CO Cone* 

Mean CO Cone. 

Mean CO Cone. 

(km) 

No. 

(Part, /cm ?) 

(Part. /cm 

(Patt. /cm^) 

78-3 

. 43.0 

4.885E+37 

4. 686E+07 

5. 129E+37 

76.3 

39.3 

6 . 6 7 3 F> 0 7 

6. 63 6E+07 

7.00 7E+0 7 

74.0 

38.3 

9.014E+07 

9.878E+07 

9. 465E+07 

72.0 

3 7.0 

1.203E+38 

1 . 10 5E+08 

1 .2 63 E+08 

73 • 3 

3 6*0 

1 . 590 E+08 

1.62 6E+38 

1 . 670 E+08 

68.0 

35. 3 

2 • 0 79 E+08 

2 . 0 1 7E+08 

2. 1 83E+38 

66.3 

34.0 

2. 695E+0B 

2* 655 E+08 

2.833 E+08 

64.0 

33.0 

3* 4 6 5E+38 

3. 423 E+08 

3* 633 E+0 8 

62 . 3 

32.0 

4. 420 E+08 

4*2 53 E+08 ' 

4*641 E+03 

60.3 

31.0 

5* 620 E+0 8 

6. 0 43 E+08 - 

5. 90 1 E+08 

58.0 

30.0 

7. 1 95E+38 

7. 108E+08 

7. 5 54 E+08 

56.0 

29.0 

9 •' 1 72 E+08 

9*382 E+08 

9. 631 E+08 

54. 0 

28-0 

1 • 1 66E+39 

1 . 1 4 6 E+09 

1 . 224E+ 39 

52.0 

27.0 

1 . 4 77 E+0 9 

1.53 6 E+09 

1 . 551 E+09 

53.0 

26.0 

1 • 8 87 E+0 9 

1 . 796£>09 r 

1 • 9 82 E+09 

48.0 

25.3 

2.417E+09 

2.383E+39 

2. 538E+39 

46*3 

24.0 

3-111 E+09 

3.0 51 E+09 

3 * 2 66 E+09 

44. 0. 

23*0 

4.089E+0? 

4.0 31 E+09 

4. 29 3 E+09 

42.3 

22.0 

5. 402 E+09 

5*1 7 6 E+09 

5. 6 73 E+0 9 . 

43 . 0 

21-0 

7. 189 E+09 

7. 765E+09 

7. 548 E+09 

38.0 

20.0 

9. 623 E+09 

1 * 0 1 2 E+ 1 3 

1 . 0 1 0 E+ 1 0 

3 6.0 

19.0 

1 .296E+10 

1 * 294E+10 

1.361E+10 

34.3 

18.0 

1 • 7 60 E+ 1 0 

1.809E+10 

1 . 848E+ 1 0 

32.0 

1 7.0 

2. 40 5E+10 

2. 266E+ 10 

2« 52 5E+ 10 

30.0 

1 6.0 

3.283E+10 

3.31 7E+10 . 

3 » 4 4 7 E+ 1 0 

28.0 

1 5.0 

4* 466E+ 1 3 

4. 2 92 E+ 10 

4 . 690 E+ 1 0 

26.0 

1 4-0 

6.09 4E+10 

6.O80E+ 1 3 

6. 399 E+ 10 

2 4. 0 

13.3 

8 • 33 5E+ 1 3 

8.282E+10 

8 • 7 52 E+ 1 0 

22.0 

12.0 

1 • 1 44E+ 1 1 

1 . 1 42 E+ 1 1 

1*231 t> 1 1 

23.0 

11.0 

1.574E+11 

1 • 579 E+l 1 

1 • 6 53E+ I 1 

1 8.3 

13.0 

2 • 1 63 E+ 1 1 

2* 334E+ 1 1 

2 • 2 7 1 E+ 1 1 

1 6.0 

9.0 

2.9 58E+1 1 

2. 874E+ 1 1 

3. 1 3 6E+ 1 1 

14.0 

8.0 

4 • 0 50 E+ 1 1 

4. 1 92 E+ 1 1 

4. 2 52 E+ 11 

12.0 

7.0 

5. 542 E+ 11 

5. 706E+ 1 1 

5. 81 9 E+ l 1 

10.0 

6-3 

7 • 58 6 E+ 1 1 

7.687E+11 . 

7.9 65E+1 1 

8.0 

5.0 

9. 71 1 E+l 1 

9. 8 67 E+ 1 1 

1 . O20E+ 12 

6*9 

4*3 

9. 839 E+ 11 

9 . 62 3 E+ l 1 

1 • 3 33 E+ l 2 

4* 0 

3.0 

9. 1 8I3E+J 1 

9.394E+ I j 

9. 639E+1 1 

2*0 

2.0 

7. 559 E+ 11 

7. 61 1 E+l 1 

7. 937E+ J 1 

3 . 

1*0 

4. 621 E+ 1 1 

5.0 51 E+ 11 

4. 852E+ l 1 
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total CO in the path is high or low by a fixed percentage. The effects of a 
systematic error of ± 10% are shown in Figure 7.4.5. The computed concen- 
tration profile deviates from that computed with no systematic error by the 
same ± 10% as is present in the instrument measurement of the total CO. 
Similarly, the computations with systematic errors of ± 2% and ± 5% yield 
calculated concentration profiles which are in error by ±2% and ± 5%, re- 
spectively. 

The computations indicate that instrumental errors of the order 
of 10% in the measurement of the total CO in the horizontal sight path through 
the atmosphere do not result in significantly larger errors in the computed 
CO concentration profile at high altitudes . While the measurement errors 
are more significant in the case where an inversion in the CO concentration 
profile occurs (i. e. the low-level sink), they do not invalidate the limb trans- 
mission experiment at the higher altitudes. The limb experiment was pri- 
marily intended for searching for an upper atmosphere sink, and in this case 
the effects of reasonable measurement errors are not of sufficient magnitude 
to permit the upper atmosphere sink profile to be mistaken for the standard 
constant mixing ratio CO profile. . , ■ 

7. 5 Multi-line Model Calculations 


7.5.1 The Program : - The first step in the theoretical determination 
of the feasability of using the correlation interferometric technique for the 
measurement of CO levels in the atmosphere is the computation of theoretical 
transmission spectra of the atmosphere in the wavelength regions of interest. 
Theoretical spectra were calculated by a fairly elaborate computer program, 
developed by GE. This program has been described in greater detail by Alyea 
and Liebling elsewhere (Bortner, Grenda, et al, 1971) and will be described 
only very briefly here. The computer program has the capability of calcula- 
ting atmospheric transmission spectra for up to three atmospheric constitu- 
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ents, with a total of up to 150 absorption lines, for three different viewing 
geometries: 

1. Ground black-body radiation including atmospheric 
absorption and emission. 

2. Solar radiation reflected from the ground including 
atmospheric absorption only. 

3. Limb transmission of solar radiation including 
absorption and emission. , 

Based on input molecular and atmospheric properties, the Voight 
profile function, absorption and emission coefficients are computed as a 
function of altitude for all lines which contribute at a given frequency value. 
Effects of overlapping lines are considered up to ± 5 cm These coeffi- 
cients are then integrated in accordance with the basic radiative transfer 
equation described in Section 4 to produce observed spectra. 

Additional subroutines, patterned after Cooley (1965) and 
Gentleman (1966) then computes the Fourier, transforms of the resultant 
spectra. These Fourier transforms which are closely related to interfero- 
grams form the data base for the theoretical feasibility study. 

This computer program was used to generate a large number 
of theoretical spectra and transforms of CO in the presence of H^O. Spectra 
were generated in the 2. 3 p spectral region (using geometry case 2) and in 
the 4. 6 p region (using geometry case 1). The regions covered were 4230 - 
4330 cm and 2100 - 2200 cm respectively. Spectra included combina- 
tions of CO models 1, 2, 3, and 6; H^O models 1, 2, and 3; temperature 
models 1, 2, 3, and 4, as described in Section 2. Ground emissivity was 
varied in the 4. 6 p spectra, and albedo was varied in the 2. 3 p spectra. In 

i i 

all, 53 theoretical spectra and transforms were computed at 2.3p, and 19 

( 

at 4. 6 p. Although there are many additional combinations which would have 
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been desirable, extremely long running time of the program (about 10 - 15 
minutes on a GE 635 or IBM 3 60 )i makes a full gamut of combinations prohib- : 
itively costly. t 

Following generation of these theoretical Fourier transforms, 
analyses were performed by use of a weighting function computer program. 
This program, patterned after a program furnished, by Barringer Research^ 
is an implementation of the weighting function concept described in detail in 
Section 6. 

The weighting function can be considered a correlation function, 
which, when multiplied by the Fourier transform of a spectrum. gives a mea- 
surement of the CP level. Inputs to the weighting function program consist 
of a basic set of several Fourier transforms of spectra generated from a. con- 
stant amount of CO and various combinations of other variables, e.g. tem- 
perature profile, H^O level, etc. One of these transforms is designated as 
the nominal case. One additional transform is input, designated the target 
case, which is identical to the nominal case with the exception of a variation 
in the CO level. The weighting function generated from these cases is then 
applied to other transforms which were, not included in the basic set. It is 
this operation of generating weighting functions from one set of data and ap- 
plying them to another set which determines the feasibility of the correlation 
operation. 

The theoretical, spectra Fourier transforms described above 
were used in various combinations as basic sets, for the weighting function 
program. Two groups of studies were made, at 2.3 p, and 4.6 p . The results 
of these studies are described below. ■ 

Calculations were made using the principles outlined in Sections 
5 and 6. The indent of these calculations is to determine the sensitivity of the 
method in the determination of CO densities, the accuracy, to be expected in 

^ The program was written by G. Levy of BRL. 
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the measurements, the- effects of various atmospheric parameters, the opti- 
mum wavelength region to be used,' and in general., to establish the feasibil- 
ity of the technique in obtaining the data needed to determine the CO sink. 

While the theoretical feasibility studies will be necessary, the 

establishment of the feasibility will only be accomplished after the instru- 

.{ 

ment is shown to be' capable: of measuring CO in the atmosphere. Thus, tests 
of the breadboard will be made' on CO and atmospheric gases to see if CO can 
be measured in amounts similar to those in the atmosphere together with 
amounts of some interferents- which, where possible/ are similar to those of 
the atmosphere:* This is hot an attempt at simulation but'an attempt to see if 
CO can be measured with the breadboard and to investigate the signal-to-noise 
ratio, under these conditions. Tests in urban polluted atmosphere will also 
be used to compare breadboard measurements with those of other methods. 

7.5.2 Calculations : - 

' 7. 5.2. 1 • Explanation of Calculations : - The results are presented in 
tables later. It must be remembered in examining the results that the object 
of this work is to find the CO sink. Thus, it is important to be able to see an 
effect produced by a rather small variation from normal CO column density, 
since, although at ground level the concentration may drop off by a large fac- 
tor, the total column density may only have small drops, as little as about 
10%, Further non-sink areas should not appear as sinks and the accuracy 
for very low GO column density may be poor, as long as it- is good enough to 
show it to be less than the norrhal atmospheric amount. It must be empha- 
sized that the most accurate measurements are needed at densities near those 
of a normal CO profile. 

Further the results should not be affected significantly by 
variations in atmospheric or" ground properties , e. g. atmospheric tempera- 
ture profile, temperature inversion layer, ground temperature, ground 
emissivity, reflectivity, and the shape of the CO profile. 
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One of the objectives of this work is to establish the types of 
conditions of the runs to be used in determining the weighting functions. 

Thus we are looking for the case giving the best results in terms of the con- 
siderations noted above. 

The calculated accuracy of the calculations are dependent on 
the conditions which were used for determining the weighting functions. The 
results which are presented in the following tables will show this. Each ta- 
ble, noted by a letter designating the case, gives the results for all interfer- 
ograms, using a specified set of interferograms to determine the weighting 
functions and then using these weighting functions to determine the CO den- 
sity in all the interferograms. The band (wavelength) is given together with 
the conditions for each run. 

Tables 7.5. 1 through 7.5.21 give the results. Each table 
gives the percentage difference of the calculated value from the actual amount 
of CO in the model. The three- right hand columns give data on this model, 
and are for all cases between the heavy lines. The eight CO model numbers 
are those of Tables 2. 1. 1 and Figure 2. 1. 1. The number in the far right 
column is the number by which the concentrations of the models (over the 
entire altitude range) are multiplied. The resulting total optical thickness 
in atmosphere-centimeters is noted. Five water models were used. Models 
W 1 , W2 , and W3 correspond to water models 1, 2, and 3 of Table 2.1.2 and 
Figure 2. 1.5; that noted as 4W 1 is that where water concentrations are four 
times those of water model 1 at all altitudes; that noted as 1. 5W2 is that 
where water concentrations are 1.5 times those of water model 2. The col- 
umn headings note the four corresponding temperature models of Table 2. 1. 3 
and Figure 2.1.6. - 

Each sheet presents the results using a specific set of runs 
to obtain the weighting functions and gives the accuracy of the calculated CO 
for each interferogram of either the 2.3 or 4. 6 y, band as noted. The sets 
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TABLE 7. 1. 15 EFFECT OF SURFACE REFLECTIVITY 
; ; ; - 2. 3 (J, BAND 


CASE 



0. 1 ' 

1.0 

A 

5.74 

5.72 

B , 

7. 99 

7.96 

c 

0. 03 

- 0. 01 

D 

- 0. 02 

0.00 

E 

4.68 

4.69 

F 

- 10. 03 

- 10 . 04 

g . 

- 18. 44 

- 18.42 

H, 

0. 02 

0. 00 

i 

1.31 

1.31 

j 

-4.38 

- 4. 39 

K 

- 19. 48 

- 19. 46 

L 

0. 00 

0. 00 

M 

-0.02 

0. 00 

N 

- 0. 04 

0. 00 
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table 7. 5 - 21 


effect 8 
parameter VAR 


SUB.EACE 

BAl^E 


parameter 


Gro 


un d Temperature 


Emi ssivity 


Gas 


Temp erature 


CASE 


P 

Q 

R 

S 

T 


288.2 

1.9 
- 5. 1 

5.2 

3.2 

_ 25.2 
0.7 


288.2 

0.0 
- 5. 1 
0.0 
o;o 
0.0 

1.0 


P 

Q 

R 

S 

T 


P 

Q 

R 

S 

T 


6.9 
.5.4 
19- *> 
11.7 
.93.8 


0.0 
.5. 1 
0.0 
0.0 
0.0 



T 


2 


0.0 

.4.1 

o.o 

0.0 

. 0.2 


0.0 
. 5. 1 
0.0 
0.0 
0.0 


0.98 


1.0 


Rel atwe 

pressure 


P 

Q 

R 

S 

T 


. 1.4 
.2.4 
.5.2 
.2.3 
18.3 


0.0 
- 5. 1 
0.0 
0.0 
0.0 
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used. in. obtaining the, weighting functions are .indicated by the rectangles j , ,, 
around, the accuracy result. . It is not, to be, expected that the results t : for 
either band will be good in all, cases and is- not -a matte, r for, concern that. . 
they are not. , It should be noted that the major ( concern is to identify a sink., , 
and to identify a source as a source. It is then possible to use weighting 
functions. derived from more appropriate conditions , to more accurately de- 
termine the CO, density. This is espepiallyj true,for , a, source region where 
the use of ,a high CO. condition. for the target condition in .establishing, weighty • 
ing functions gives the high CO concentrations much more accurately.. This, 
can be seen in Case N.- .. ,, , • ; , ,. ( ... , , . • - 

7. 5. 2. 2 The Overtone Band: - The results for the analysis are pre- 
sented in Tables 7. 5. 1 through 7. 5. 15 (cases A through N) for the overtone 
(2. 3 p,) band. Some of the results show accuracy which will fit the accuracy 
requirements detailed above while some show inaccuracies which are too 
large to be acceptable. This is as expected. 

Cases A, B, C, and D use the same CO models for the base 
and target cases used for determining the weighting functions . In cases A, 

B, and C, all runs give accuracy suitable for use with one possible exception 
(water model 4W 1 , Tl). Except for this one set of conditions, the model 1 
cases (standard CO profile), none of the runs show any significant indication 
of a sink or source. Note that the model with a temperature inversion is 
quite accurate. For the important situation of the low altitude sink, all runs 
show good accuracy, all sufficient to show a lower total CO column density 
than the standard model, that is, a sink. The models with low CO densities 
show large errors but not so large that sinks are not shown. The model with 
1. 5 times the standard CO model gives accuracies which all show significant 
increases over the standard CO. For the models with high CO densities the 
results are inaccurate, up to about 50%, but such still indicate high CO den- 
sities and thus fulfill the purpose. For case D the results are not nearly so 
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good. The conditions used to determine weighting functions do not include 
any of the low water density. models . C)ne test run with such a model shows, 
for the standard CO model, a doubling of what it should be. This would indi- 
cate a source and so is unacceptable. . In general,' the errors are too large 
to be acceptable. 

Case E shows the effect of using the very low CO density 
cases as the base runs and one standard CO model as the target run to obtain 
the weighting functions. The results are much the same as cases A and C. 

In this case all results are suitable. Cases F and G show results for similar 

y 

CO models being used in determining weighting functions but with more re- 
stricted temperature conditions. The results are poor.;. Sinks are indicated 
where there are no sinks. 

Case H is another case where the conditions of temperature 
and of water content did not cover a sufficient range. Thus large sinks are 
indicated where there is no sink. 

• • 1 • • ' : 

,Now consider cases I, J, and K, where all runs used in de- 
termining weighting functions were, with models having low CO densities (. 1 , 
.2, and . 3 of standard). For the case, I, where the ranges of temperature 
and water density employed'were wide, the, results are good, showing sinks 

' ‘ " ’ 1 f . •• • fi 

and sources where they occur. With more Restricted ranges, Cases J and K, 
the results are not good, indicating sinks where they do not exist. 

For Case L, the weighting functions were obtained using the 
standard and the low-altitude sink models , but only one water model. Re- 
sults for other water models are not accurate enough to be used. 

For case M, simU ar results are obtained to case L. No runs 
with low water densities were Used for weighting function obtention and re- 
suits for such models were inaccurate. ■- 

For Case N, where a high CO density model was used in the 

' ' '• 7 * 1 . 5 

set for obtaining weighting functions, just about all test runs were inaccurate. 

. ’ ' •' i 
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Such results are to be expected since the CO absorption curve is non-linear 
in part of the region used to determine weighting functions. 

In Table 7. 5. 15 it is seen that a change in the earth's reflec- 
tivity has no significant effect on the calculated CO density. It would also be 
small for the 4. 6 p band. 

Some general conclusions can be drawn from these results. 
They mainly concern the conditions used in obtaining weighting functions. If 
such, for CO are within the range of CO where absorption is about linear, 
that is within the range from very low CO densities to slightly above our 
standard model, the results will be good within the range of water density 
and temperature covered in the weighting function runs and probably not 
good outside such range. That is, if water density models 1 and 3 are used, 
those for in-between water densities have the desirable accuracy, but if 
models 2 and 3 are used, the results for water model 1 are inaccurate. That 
is, interpolation gives reasonable results but extrapolation does not. 

It is important to note that the variation in total water content 
does not vary over one area by a factor of more than about four during the 
year as seen in Figures 7. 5. 1 through 7. 5.4. Thus the extremes used here 
are wider than those for any one region. Hence, in treating the data for any 
one area, weighting functions based on the range of conditions for that area 
can be used rather than those which bracket the range of conditions for all 
areas . 

Conclusions about sinks and sources can also be made. 

These results show that sinks (even those at low altitude having only a 9% 
lower CO optical thickness than the standard model) are detected readily and 
regions of high CO concentration are also readily seen. It is also to be noted 
that the presence of a temperature inversion layer does not interfere with the 
results for this (2. 3 p) band. 
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Figure 7.5.1 January Mean Precipitable Water (in inches) from Tuller - 1968 
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Figure 7.5.2 April Mean Precipitable Water (in inches) from Tuller - 1968 
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Figure 7.5.3 July Mean Precipitable Water {in inches) from Tuller - 1968 









Figure 7.5.4 October Mean Precipitable Water (in inches) from Tuller - 1968 







7. 5.2.3 ‘ The Fundamental Band : - The results for the analysis for the 
fundamental (4.6^) band in Tables 7. 5. 16 through 7.5.21 (Cases P through T). 
The cases employ weighting functions obtained using CO model 1 as the base 
runs with 1. 5 x CO density of model 1 (Case S), urban CO density model 
(Case T), very low. CO density model (Case P), and the low-altitude sink 
model (Case R) for the target runs. As would be expected the best results 
are obtained for Cases P and R. However, it is to be noted that in all cases 
the temperature inversion model gives poor results with sinks not detected. 
This is an intolerable result. If CO densities for such critical conditions are 
in error by 26% (Case P) or by 71% (Case R) existing sinks will not be found. 
Cases P and R are thought to be as reasonable as any tests that could be 
made for this band and other tests would be expected to show similar results. 

Case O uses the temperature inversion layer run as a part of 
the base set in determining the weighting functions. This is not a very prac- 
tical case but it does improve the results as a function of temperature. 

There is, however, still an overlap in the CO low altitude sink and the CO 
standard model results which indicate that it could not be used to see a sink. 

In all 4. 6 p, cases, the source runs show up with large nega- 
tive errors. In most all urban model runs, negative CO densities are calcu- 
lated. These are meaningless and undoubtedly smaller sources would show 

up as sinks. This' is an intolerable situation. 

* 

Certain other tests were made to help show the effect of at- 
mospheric and earth- surface parameter variations. Results are shown in 
Table 7.5.21. 

Significant effects are found, even for the runs, P and R, 
which use better sets of runs for determining weighting functions. A few de- 
grees difference in the ground and the 0 km atmosphere temperatures cause 
a few percent difference in calculated CO. A reduction of surface emissivity 
to 0.7 causes drastic effects. Changing the temperature model of both the 
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entire atmosphere and the ground by 2° produces no significant effect. This 

o 

is as expected. However, a random variation with altitude of i 2 would pro 
duce large effects as seen in Section 7. 1.4 with the single-line model. A 
change in the total pressure significantly changes the calculated CO since it 
reduces the atmospheric emission relative to the earthshine. The effects of 
these variations for the 2. 3 y, band were shown in Section 7.1 to be small. 

7. 5. 2. 4 Summary : - The following summarizes the results for the 
two bands based on the calculations with the multi-line model. They are 
similar to those obtained with the single-line model. 

The conclusions for the two wavelengths studied may be sum 
marized on the basis of the following considerations. 

Optical Thickness: 

4. 6: For much of the range, the absorption curve is not linear. 

The higher optical thicknesses may be so high that there 
is little sensitivity in that range. 

2.3: The optical thickness is suitable for the range of concen- 

trations of interest. The only limit is the upper limit of 
the limb experiment but it will be suitable for a sink at 
any reasonable expected altitude. 

Interpretation: 

4. 6: Very difficult because of the need for accurate atmospheric 

temperature data as a function of altitude, ground tempera- 
ture, and ground emissivity at 4. 6 |j,. 

2. 3 Data directly presented in terms of CO density. 
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Sensitivity to Low Altitude Sink:* . r 

4. 6 Sensitivity to low altitude CO is very low and errors intro- 
duced are likely to be larger than the CO decrease effect 
being sought. 

,2.3: . Sensitivity at low altitude same as the similar effect at 

high altitude. Sinks of less than 10% CO density decrease 
... can be detected. 

• , - . • • • ’ ■ • i ‘ ' ‘ ' ' * . • • • 

Temperature Inversion Layer: 

,4. 6: • Drastic effects on calculated CO densities prohibit detection 
of significant sinks . . ... 

2. 3: No appreciable effect. 

Line Strengths: 

4. 6:.- Strong. 

2.3: About 1% of 4. 6 lines. This is strong enough. 

Atmospheric Emission: 

4. 6: Considerable emission in this band causes, much variation 

from straight absorption model and leads to difficulty in 
interpretation as mentioned. 

2. 3: Emission negligible. 

Night Use: . 

4.6: Possible. '• 

2. 3: Not likely to be possible, at least with acceptable sensitivity. 
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Mapping Experiment: 


4.6: Emis sion effects disastrous. 

2.3: Interferents cause problems but these are overcome by 

technique used. . . , 

Limb Experiment: 

4. 6: Atmosphere optically thick at lower altitudes of interest. 

2. 3: Sensitivity limits altitude but altitude range reasonable. 
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CONCLUSIONS 


The correlation interferometric techniques has been shown by analysis 
to be capable of measuring atmospheric trace species accurately. Specifi- 
cally it has been shown to measure atmosphere amounts of carbon monoxide 
with ah accuracy of better than 10%. Thus the technique should be capable of 
picking out CO sinks which are only 10% or less lower in CO column density 
using measurements which are made by determining the absorption in a part 
of the first overtone band of CO in the 2. 3 p, region of the spectrum. It has 
been determined that the use of this band is preferable to the use of the fun- 
damental band in the 4. 6 p, region. The latter is unsuitable because of the ef- 
fect of atmospheric parameters including the atmospheric temperature pro- 
file (and the associated atmospheric emission) (see Tables 7. 1.4 and 7. 1. 5), 
the ground temperature (see Tables 7. 1.6 and 7. 1.7), the ground emissivity 
at the wave lengths being used (see Table 7. 1. 8), atmospheric pressure and 
/ atmospheric path length. For that band, atmospheric emission causes an 
enhancement of the radiation which is a function of temperature and must be 
taken into account. If at some altitude the atmospheric temperature is greater 
than the ground temperature the effect of absorption below is minimized. 
Effectively, the instrument can not see well below an atmospheric tempera- 
ture peak. Then the absorption does not follow the CO profile at low altitudes 
but rather there is less absorption near the ground even for CO profiles when 
there is more CO in that region (see Figures 7. 1.6 and 7.1.7). These diffi- 
culties arise for any species with all spectral techniques which are primarily 
involved with radiation of wavelengths greater than about 3. 5 p. This is be- 
cause of the predominance of earthshine over reflected solar radiation at 
these wavelengths (see Figure 3.2. 1). The overtone band is not affected sig- 
nificantly by variaties in the parameters noted above and absorption in this 
band follows closely the atmospheric CO profile for any reasonable atmo- 
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spheric model. Atmospheric scintillations have little effect at either wave- 
length region,. ; *,♦ *. . 

. Spectral interferents can be overcome by the correlation interferome- 
try technique. The; chief spectral interferants in the 2. 3 p, spectral region 
are water and methane. By using the interferogram directly, carbon mon-, > 
oxide can be accurately measured in the presence of atmospheric amounts of 
these gases which would prevent accurate measurement of CO. by ordinary 
spectral methods. In order to accomplish the interferometric measurement 
it is necessary to calibrate the instrument over the entire range of density of 
interferents for which it will be used and to do so with variation of important 
atmospheric parameters over which it will be used. Thus the calibration 
must cover the range of methane of about 2 to 5 atmosphere cm. since this 
is the range that is expected to be encountered and to cover the range of about 
0,2 to 3 precipitable cm. of water since variations from dry to wet atmos- 
pheres include this range. Further, since the population of the rotational 
water levels are appreciably affected by temperature, variations correspon- 
ding to changes in atmospheric temperature profiles must be included in the 
calibration, thus necessitating calibration for conditions of a cold dry, a hot 
dry, a cold wet, and a hot wet atmosphere as well as conditions in between 
the s e . 

The calibration determines the weighting function which multiplies the 
section of the interferogram in such a way as to minimize the effect of spec- 
tral interferents and maximize the effect of the gas to be measured (see 
Chapters 5 and 6). The choice of conditions used to determine the weighting 
functions is critical. The use of a wide range of conditions and interpolation 
between these gives much better results than use of a narrow range with ex- 
trapolation. In practice under flight conditions, it may be best to use a 
weighting function derived from a wide range of conditions to obtain an ap- 
proximate measure of CO density and conditions and then use a weighting 
function derived from a narrower range of density and conditions to obtain 
a more accurate CO density measurement. 
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iThe feasibility analysis thus shows the correlation interferometry tech- 
nique to be capable of the measurement from a remote platform of carbon 
monoxide in the atmosphere over the desired range of density. The measure- 
ment can be made in the mapping mode (observing sunlight reflected by the 
earth) and the limb mode (observing direct sunlight through the earth's limb). 
It can further be stated that the technique is applicable to a variety of other 
gases present in trace amounts in the atmosphere. ! 
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APPENDIX A 


LISTING OF PROGRAM SPECTRA 
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C0CS33 480 CONTINUE 

000536 EMISCtIH) = EMSUM 

CCC540 490 ABSeCCIH = A8 SUM 



ABSORPTION AND EMISSION COEFFICIENTS COMPUTED FOR ALL CONTRIBUTING 
LINES AT ALL ALTITUDES 

NOW INTEGRATE OVER ALTITUDE USING SIMPSCN*S RULE 
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C INTERFEROMETER PLOT SCALE 
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APPENDIX B 


LISTING OF CORRELATION FUNCTION PROGRAM 
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